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Preface 


The  purpose  of  this  study  was  to  derive  a  method  of 
selecting  a  radar  waveform  which  would  satisfy  a  given  set 
of  criteria  for  target  detection  in  a  specified  environment. 
The  method  derived  in  this  study  was  to  calculate  and  graphi¬ 
cally  display  the  magnitude  of  the  radar  ambiguity  function 
for  radar  waveforms  which  met  some  general  criteria.  After 
analyzing  these  magnitude  displays,  those  which  approximated 
the  resolution  requirements  for  target  detection  were  then 
analyzed  for  their  performance  in  the  specified  environment. 

To  assist  in  this  analysis,  the  signal-to-clutter  ratio  for 
the  desired  target  in  a  model  of  the  clutter  environment  was 
calculated.  The  waveform  which  best  met  the  criteria  for 
target  detection  in  the  modeled  environment  was  then  deter¬ 
mined. 

A  limited  group  of  radar  waveforms  were  investigated  in 
this  study;  however,  all  waveform  types  generally  used  in 
contemporary  radars  were  investigated.  The  ambiguity  func¬ 
tions  were  calculated  as  a  function  of  range  delay  and  Doppler 
frequency.  Calculations  of  ambiguity  functions  have  also 
been  made  as  a  function  of  azimuth  angle  and  elevation  angle, 
but  these  are  not  significant  in  selecting  the  radar  waveform 
and  were  not  considered. 
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Abstract 

‘"/The  magnitude  of  the  radar  ambiguity  function  for  com¬ 
monly  used  radar  waveforms  is  calculated  and  displayed  in 
three-dimensional  diagrams.  The  ambiguity  functions  are  cal¬ 
culated  as  a  function  of  range  delay  and  Doppler  frequency. 
For  a  specified  clutter  environment  and  a  particular  radar 
waveform,  the  signal-to-clutter  ratio  for  a  desired  target 
is  calculated.  A  method  to  select  a  radar  waveform,  given 
the  desired  target  detection  criteria  and  clutter  environ¬ 
ment,  is  derived.  This  method  relies  on  the  analysis  of  the 
magnitude  of  the  ambiguity  function,  signal-to-clutter  ratio, 
and  other  basic  principles  of  radar  theory.  /t\ 
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I. 


Introduction 


Statement  of  the  Problem 

A  significant  problem  encountered  by  a  radar  engineer 
is  the  choice  of  a  radar  waveform  which  provides  the  desired 
detection  of  targets  in  some  given  environment.  This  problem 
is  one  of  determining  which  waveform  can  provide  the  desired 
resolution  in  boch  range  and  Doppler,  the  desired  detection 
of  targets,  minimal  ambiguities  in  the  measurement  of  target 
parameters,  and  satisfactory  clutter-rejection  properties  for 
the  specified  environment  (Ref  11:420).  The  selected  wave¬ 
form  can  only  be  implemented  if  it  can  be  transmitted  and 
received  within  the  size  and  budget  constraints  of  the  system. 
By  calculating  the  ambiguity  functions  of  radar  waveforms  and 
analyzing  the  data  from  them,  the  radar  engineer  would  have 
most  of  the  information  he  needs  to  make  such  a  determination. 
If  he  also  could  determine  the  signal-to-clutter  ratio  for 
the  desired  target  in  the  specified  environment,  he  would  have 
all  the  necessary  tools  to  solve  the  problem. 

Providing  a  fast,  simple,  and  accurate  method  for  se¬ 
lecting  a  radar  waveform  which  satisfies  the  targe-  detection 
requirements  in  a  given  environment  .is  the  purpose  of  this 
study.  This  method  relies  on  the  ability  to  calculate  the 
magnitude  of  ambiguity  functions  for  various  waveforms,  calcu¬ 
late  the  signal-to-clutter  ratio  for  the  specified  target  and 
environment,  and  analyze  these  results. 
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Background 


In  Gaussian  noise,  the  signal  ambiguity  function  plays 
an  important  part  in  obtaining  a  maximum  likelihood  estimate 
of  signal  parameters.  In  fact,  in  white  noise,  the  estimate 
is  obtained  by  constructing  the  ambiguity  function  as  a  func¬ 
tion  of  the  parameters  to  be  estimated  and  the  values  which 
maximize  the  function  are  taken  as  the  desired  estimate. 
Woodward  has  shown  how  this  is  done  for  delay  and  Doppler  and 
has  derived  formulas  for  the  variances  of  the  estimates.  For 
more  general  maximum  likelihood  estimates,  Kelly,  Reed,  and 
Root  introduced  a  more  general  version  of  the  ambiguity  func¬ 
tion,  which  in  the  case  of  white  Gaussian  noise  is  essentially 
proportional  to  the  Woodward  ambiguity  function  if  the  unknown 
parameters  are  delay  and  Doppler  (Ref  14:39). 

Since  the  concept  of  the  ambiguity  function  was  first 
introduced  in  1953,  a  multitude  of  research  has  been  performed 
in  the  field.  The  ambiguity  functions  for  a  large  number  of 
waveforms  have  been  calculated,  and  the  results  are  well  docu¬ 
mented.  Several  authors  have  provided  methods  of  selecting  a 
radar  waveform  based  on  analysis  of  ambiguity  functions.  Of 
these,  some  have  provided  only  limited  criteria  for  the  range 
and  velocity  resolution  for  a  waveform  without  much  regard  to 
the  clutter  environment  (e.g.  Sinsky)  (Ref  10).  Others  have 
attempted  to  synthesize  ambiguity  functions  to  approximate  a 
function  which  would  provide  the  optimum  performance  in  the 
clutter  environment  and  in  resolution  capabilities  (e.g. 
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Ref  13) .  The  problem  is  that  the  techniques  suggested  up  to 
now  are  either  too  simple  to  provide  accuracy  or  too  compli¬ 
cated  to  provide  the  radar  engineer  with  a  fast  method  to 
choose  a  radar  waveform. 

Approach  to  Solving  the  Problem 

The  method  to  choose  a  radar  waveform  suggested  in  this 
study  is  based  on  the  ability  to  quickly  and  accurately  calcu¬ 
late  radar  ambiguity  functions,  graphically  display  the  ambi¬ 
guity  function  magnitude,  and  calculate  the  signal-to-clutter 
ratio.  The  mathematics  required  to  perform  the  necessary 
calculations  dictated  the  use  of  a  computing  system  with  a 
high-level  language  capability,  a  main  memory  of  at  least  32K 
bytes,  access  to  a  Fast  Fourier  transform  program,  and  a 
method  of  providing  three-dimensional  diagrams.  During  this 
study,  the  Data  General  Eclipse  computer  was  used  since  it 
met  all  the  necessary  requirements.  Fortran  V  was  the  pro¬ 
gramming  language  selected  because  it  provided  the  fastest 
and  most  efficient  means  to  perform  the  calculations  on  the 
Eclipse. 

In  order  to  calculate  the  ambiguity  function,  the 
approach  in  this  study  is  to  take  cuts  of  the  ambiguity  func¬ 
tion  for  various  values  of  Doppler.  With  this  approach,  each 
cut  of  the  ambiguity  function  is  calculated  by  performing  a 
linear  convolution  of  appropriately  transformed  samples  of 
the  waveform.  The  magnitude  of  each  sample  is  then  calculated. 
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By  appropriately  spacing  these  Doppler  cuts,  an  accurate 
representation  of  the  magnitude  of  the  ambiguity  function  is 
provided.  This  approach  also  allows  simple  application  to 
three-dimensional  display  techniques.  Since  the  magnitude 
data  are  sampled  in  both  delay  and  Doppler,  the  data  can  be 
easily  formatted  for  use  by  most  available  three-dimensional 
plotting  programs. 

In  addition  to  the  magnitude  of  the  radar  ambiguity 
function,  the  calculation  of  the  signal-to-clutter  ratio  also 
has  several  other  required  inputs.  The  specified  environment 
must  be  modeled  by  a  probability  density  function  of  clutter 
as  a  function  of  delay  and  Doppler.  Additionally,  the  average 
radar  cross  section  of  the  desired  target  and  the  clutter  must 
be  determined  as  well  as  the  number  of  different  scatterers 
providing  radar  returns  in  the  environment.  With  these  in- 
puts,  an  accurate  approximation  of  the  signal-to-clutter  ratio 
is  calculated. 

The  method  of  radar  waveform  selection  in  this  study 
considers  the  range  and  velocity  resolution  requirements, 
blind  speed  limitations,  unambiguous  range  requirement,  and 
the  signal-to-clutter  ratio  requirements  for  the  radar  system. 
With  the  capability  to  calculate  and  display  the  ambiguity 
function  magnitude  and  to  calculate  the  signal-to-clutter 
ratio,  all  of  the  selection  criteria  can  be  evaluated. 


( 
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Limitations 


This  study  was  limited  in  the  types  of  radar  waveforms 
considered,  the  ambiguity  function  variables  considered,  and 
the  amount  of  target  motion  allowed.  To  begin,  only  those 
radar  waveforms  generally  used  in  contemporary  radars  were 
considered.  These  waveforms  included:  simple  pulse,  pulse 
with  linear  frequency  modulation  (FM) ,  pulse  with  V-FM,  phase- 
shift  coded  pulse,  and  repeated  pulses.  The  majority  of  wave¬ 
forms  ever  seriously  considered  for  current  use  are  contained 
in  this  group. 

Although  a  general  version  of  the  ambiguity  function  is 
available,  the  parameters  used  in  radar  ambiguity  functions 
are  range  delay,  Doppler,  azimuth  angle,  and  elevation  angle. 
Since  the  ambiguities  of  azimuth  and  elevation  angles  effect 
the  antenna  design  of  the  system,  they  must  be  considered  by 
the  radar  engineer.  Those  angles,  however,  do  not  signifi¬ 
cantly  affect  the  waveform  selection  process  in  the  narrowband 
case  and  are  not  considered  in  this  study  (Ref  15:2104). 

In  addition  to  the  above  limitations,  this  study  con¬ 
siders  only  the  case  where  there  is  negligible  distortions 
of  the  complex  envelope  of  the  waveform  by  target  motion, 
target  acceleration,  and  higher  order  derivatives.  This  case 
does,  however,  include  a  wide  range  of  practical  applications. 
When  this  is  not  the  case,  the  resolution  theory  becomes  much 
more  complex.  The  theory  has  been  extended  by  many  authors 
including  Rihaczek  (Ref  5: Chapter  11  and  Ref  8). 


Presentation 


To  provide  some  insight  into  the  method  of  selecting  a 
waveform,  this  study  begins  by  discussing  the  significant 
theory  associated  with  the  radar  ambiguity  function.  The 
concept  of  the  complex  signal  representation  of  a  radar  wave¬ 
form  is  discussed.  The  development  of  the  matched-filter 
receiver  and  the  optimum  receiver  in  white  Gaussian  noise  is 
then  presented.  Properties  of  the  ambiguity  function  are 
considered  as  well  as  the  uncertainty  relation  between  delay 
and  Doppler.  A  short  discussion  of  the  digital  signal  pro¬ 
cessing  techniques  used  to  compute  the  ambiguity  function  is 
included.  Some  characteristics  of  clutter  for  particular 
environments  and  some  effects  of  clutter  on  detection  are 
also  presented. 

This  study  continues  with  a  discussion  of  the  software 
developed  to  help  resolve  the  problem  of  waveform  selection. 
The  method  of  calculating  cuts  of  the  ambiguity  function  for 
various  Doppler  frequencies  is  explained.  Use  of  the  three- 
dimensional  graphic  technique  to  produce  diagrams  of  the 
magnitude  as  a  function  of  delay  and  Doppler  is  discussed. 

The  requirements  and  method  used  to  calculate  the  signal -to- 
clutter  ratio  ends  this  section  of  the  study. 

After  generating  the  software  products,  a  methodology 
of  analyzing  them  is  needed.  This  study  provides  a  discus¬ 
sion  of  a  methodology  and  how  it  was  developed.  Some  general 
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conclusions  on  how  to  select  a  waveform  and  some  recommenda¬ 
tions  for  further  development  of  the  software  end  this  study. 
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II .  Theoretical  Considerations 


Complex  Signal  Representation 

Radar  signals  consist  of  a  high-frequency  carrier  modu¬ 
lated  in  amplitude  or  phase  by  a  function  that  varies  much 
slower  tnan  the  carrier  frequency.  These  signals  can  be 
represented  by  the  equation: 

S(t)  =  M  { t )  cos  (u)  t  +  <t>(t))  (1) 

o 

where: 

M  (t)  is  an  amplitude  modulation  factor. 

is  the  carrier  frequency  (in  radians)  . 

<p  (t)  is  the  phase  modulation  (in  radians)  about  the 
carrier. 

This  equation  can  be  represented  in  quadrature  form  as: 

S(t)  =  X  (t)  cos  uj  t  -  Y  (t)  sin  w  t  (2) 

o  o 

where : 

X(t)  =  M  (t)  cos  <p  (t)  . 

Y  (t)  =  M  (t)  sin  $  (t)  . 

We  can  use  the  quadrature  components  to  form  a  complex  func¬ 
tion,  u  (t) ,  which  is  defined  as  the  complex  envelope  of  the 
signal  S  (t) : 

y  (t)  =  X  (t)  +  j  Y  ( t)  =  M(t)ej<1>(t)  (3) 

where: 
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M  ( t ) 


(X2  (t)  +  Y2(t))>S. 


<t>  (t)  =  tan"1  (Y(t)/X(t)  )  . 

Therefore,  since 

S(t)  =  Re{y  (t)eiu>ot}  (4) 

the  complex  envelope  eliminates  the  necessity  to  manipulate 
the  carrier  frequency  in  calculations  with  the  original  signal. 

Using  the  complex  envelope  to  represent  a  signal  yields 
several  interesting  properties.  The  first  of  these  proper¬ 
ties  is  the  relationship  between  the  energy  (E)  of  the  signal 
and  the  energy  in  the  complex  envelope.  Energy  is  defined  as: 

e  =  ;”|s(t) |2dt 

=  M2  (t)  +{*  (X2  (t)  -Y2  (t)  cos2wQt-X  (t)  Y  (t)  sin2uot)  dt 
=  /“^M2 (t)  dt 

2E  =  /“ M2 (t)  dt.  (5) 

—  00 

Now  considering  the  energy  in  the  complex  envelope,  we  write: 

[2  |y (t) |2dt  =  /“y  ( t)  y * (t)  dt 

=  />2  (t)  dt  =  2E.  (6) 

Thus,  the  complex  envelope  contains  twice  the  energy  of  the 
signal  that  it  represents. 

Another  important  property  to  consider  is  the  convolution 
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of  the  complex  envelopes  of  a  signal  and  the  impulse  response 
of  a  filter.  If  we  consider  the  frequency  content  of  the 
complex  envelope  of  a  signal,  we  note: 

S(u>)  =  /~S(t)e"-i“t  dt 

—  oo 

-  H  .Ollt)ejV  +  „*<t>e--l“ot>e-J“tat 

=  h  {M(oj-o)  )  +  M*  (-io-uo)  } .  (7) 

Therefore,  we  can  represent  the  frequency  content  of  the 
filter  complex  envelope,  A  (t) ,  as: 

H  (u3 )  =  ^{L  (w-u>  )  +  L*(-u)-a)  )}. 

o  o 

Now 

h  (t)  =  ,-5{l5TTejwot/"L(f2)ej-Qtdn+i'iTTe*ju,ot/"L*(Q)e-Jntdn} 

=  ij{e-’a)otA  (t)  +  e_-'U)otX*  (t)  } 

=  Re(A(t)  e^Uo^}.  (8) 

Hence,  given  an  input  signal,  (t)  =  Re{y  (t) eJwot} ,  an  out¬ 
put  signal,  SQ  (t)  =  Re{  8  (t)  e-^oM ,  and  the  impulse  response 
of  a  filter,  h(t)  =  Re{X  (t) eJuot} ,  we  find  that: 

S  ft)  =  ;"h(T)S,  (t-T)dT  =  ^(Oe^o1  +  8*  (t)  e"Juot}  (9) 

O  -*  ”  1 

=  ^/“x(T)u(t-T)eJ“otdT  +  !$/*X*  (x ) u*  (t-T )  e“^u)otdx 

+  V"x  (x)  y*  (t-x)  e“^uo  ^2T~t*  di  +  !^/*X*  (t)  y  (t-T) 

eJ“o(t-2T)dT. 
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Since  the  third  and  fourth  integrals  equal  zero, 

SQ(t)  =  Re{^e^a)ot/0SX  ( t )  p  (t-x)  dx}  (10) 

Therefore,  we  can  equate  values  of  SQ(t)  (eqs  9  and  10)  and 
find  that: 

B(t)  =  V“X  (x)y  (t-i)di  =  A  (t)  *u  (t )  }  (11) 

Thus,  the  complex  envelope  of  the  output  of  the  filter  is 
half  of  the  convolution  of  the  complex  envelopes  of  the  input 
signal  and  the  impulse  response  of  the  filter. 

In  addition  to  deterministic  signals,  we  want  to  con¬ 
sider  the  use  of  the  complex  envelope  with  stationary  random 
processes.  If  we  consider  a  process,  n(t),  that  is  band- 
limited  between  the  frequencies  of  ±oj  around  a  carrier  fre¬ 
quency  Q ,  we  can  represent  n(t)  as 

n(t)  =  Re{n  (t)  e-*fit}  . 

We  now  consider  the  correlation  function,  Rn(T),  of  the  pro¬ 
cess. 

R  (t)  =  E(n (t)  n* (t-T) } 
n 

=  (x)  ejP‘T  +  R~*  (t )  e~^T } 

+  !*{E(ft(t)ft(t-T)  )ej  (2nt-T)+E(n*  (t)  n*  (t-T)  )e~j  (2r2t"T)  }] 

It  can  be  shown  that  since  n(t)  is  stationary,  E  (n  (t) n ( t-x ) ) 
and  E (n* (t) n* (t-x) )  are  equal  to  zero.  Hence, 
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(12) 


Rn  (x)  =  ^Re{Rfi(T)ej'QT} 

It  follows  that  the  spectral  density  of  the  process  is 

Sn(u)  =  +  Sft*(-u-S2)}  (13) 

We  want  to  also  consider  what  happens  when  we  pass  a 
stationary  process  through  a  wideband  filter.  If  we  let 
(w)  be  the  input  spectral  density,  Sq(io)  the  output 
spectral  density,  and  H(ju>)  be  the  filter  spectral  density, 
we  know  that 

SQ  (to)  =  Si  (w)  |  H  ( jco)  |2  (14) 

It  is  possible  to  represent  j  H  < j )  |  as 

|  H  ( jui)  |  2  =  [hLp(oi-Q)+hLp*(-u)-n)  ]  [hLp*(a>-n)+hLp(-u)-n)  ] 

We  can  then  show  that 

S5(w)  =  Sj( u>)  IhLp(u)  j2  (15) 

and  correspondingly  that, 

Rq(t)  =  Re{R~  (t  )  e^T }  (16) 

We  also  want  to  consider  the  case  where  the  bandpass 
noise  process  is  Gaussian.  If  we  assume  that  the  signals 
and  noise  have  passed  through  a  filter  whose  passband  in¬ 
cludes  the  signal  spectrum,  but  is  much  wider  than  the  signal 
spectrum,  the  noise  can  be  represented  by  a  complex  envelope. 
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n(t)  =  Re{n (t) eJ^t} .  It  has  been  shown  that  with  these 
assumptions. 


Rf=  (t)  =  5  (r) 

2 


(17) 


where 


Nq/2  =  amplitude  of  the  noise  spectral  density 
6(t)  =  Dirac  delta  function. 


Determination  of  the  Matched  Filter 

In  the  infancy  of  radar,  the  signal-to-noise  ratio  was 
a  popular  measure  of  the  effectiveness  of  a  radar  receiver 
for  combating  noise.  The  idea  was  to  design  the  receiver  to 
provide  the  maximum  signal-to-noise  ratio  which  made  it  easier 
to  distinguish  a  signal  from  background  noise.  Signal-to- 
noise  ratio,  however,  is  only  indirectly  related  to  radar  per¬ 
formance  measures  such  as  the  probability  of  detecting  a  tar¬ 
get,  range  and  Doppler  estimation  accuracy,  and  so  forth.  It 
should  be  noted  that,  in  many  cases,  decision  theory  leads  to 
solutions  identical  to  those  obtained  by  maximizing  the  signal- 
to-noise  ratio  (Ref  1:143). 

To  begin,  suppose  we  have  an  input  waveform,  V(t)  = 
si(t)  +  ni(t),  where:  s^t)  =  input  signal  component 

n^(t)  =  noise  component, 

and  an  output  waveform  from  a  filter,  y(t)  =  sQ(t)  +  nQ(t). 

We  want  to  determine  the  linear  filter,  h(t),  which  maximizes 
the  output  signal-to-noise  ratio  at  time  t=tm*  The  signal- 
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to-noise  ratio  is  then  defined  as: 


S  _  s  (t  ) 
N  - 


n  ‘ ‘(t) 
o  in 


Assuming  the  noise  is  white  and  Gaussian,  we  know  using  Eq  17 


that 


Hence,  we  can  write: 


R  (t)  =  6(x)  ^o 
n  2 


S  _  sj(T)  Mtm-T)di)  |2 

N  No  /fcm  h2(t  — t ) dx 

~Y  m 


If  we  consider  the  signal-to-noise  ratio  in  the  frequency 
domain  and  apply  the  Schwarz  inequality  to  the  numerator,  we 
see  that: 


s  1  \f  H(w)Si  (u))ejwtmdu)|  2  1  J*  | S±  (u) 

5  =  »o/|H(„)|2dU  £  W° 


If  H(w)  is  chosen  so  that  equality  exists,  then 


H  ( to )  =  ^  ( w)e  •'U)tm 

o 


Thus,  we  see  that  the  optimum  filter  in  white  Gaussian  noise 


h<v  = 
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A  similar  derivation  can  be  made  using  the  complex 

envelope  representation  of  the  signals  and  the  filter  impulse 

response.  We  seek  the  linear  filter,  h(t),  that  maximizes 

the  output  signal-to-noise  ratio  at  some  time,  tm.  Let  the 

input  signal  be  x(t)  =  Re{x  (t)  e-'Ci)ot} ,  the  filter  impulse 

response  be  h(t)  =  Re{h (t) e^ Wot} ,  and  the  output  signal  be 

y(t)  =  Re{y  (t)  e^o*"} .  We  also  assume  that  we  have  white 

Gaussian  noise  that  has  passed  through  a  bandpass  filter  and 

that  R  (t)  =  Re(R~  (x )  e-® WoT )  and  R-,(t)  =  ^2.6  (t)  .  We  have 
n  n  2 

already  found  that  the  optimum  receiver  in  white  Gaussian 

noise  is  K  x* (t  -t) . 

m 

Hence  we  can  say 


y  (t) 


Vx*  (t^—t+T )  x  (x)  dx 


R~ (t-t 


m 


) 


(22) 


Using  our  definition  of  the  signal-to-noise  ratio  as  before, 
and  considering  it  in  terms  of  frequency,  we  see  that: 


S 

N 


(t=tm5 


fcl/X(f)H(f  )eJ2lTftn»df  I2 

jH(f)  |  df 

4 


(23) 


Using  the  Schwarz  inequality,  we  note  that: 

s  <  / 1 x  Cf )  l2df  ;lH(f)  |2df 

N  N  / j  H ( f )  | 2df 

o 

To  get  equality,  H(f)  =  X*  (f )  e~-,u3tm,  and  hence,  fi(t)  = 
X* (tm-t)  . 
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Uncertainty  Relationship 

In  this  section  we  want  to  develop  parameters  that 
characterize  a  signal  and  form  the  basis  of  comparisons 
between  signals  to  accurately  measure  target  parameters.  We 
develop  these  parameters  in  terms  of  the  complex  envelopes 
of  signals,  S(t)  =  Re{f(t)eJ"  }.  We  begin  with  a  series  of 
def initions : 

uT  =  mean  frequency  deviation 


/t o  IF  ( lo )  I  2du) 
/|FU)  \2du> 


j  JfMt)  [^j^-]dt 
/if (t) | 2dt 


~Z 

uJ 


mean  square  frequency  deviation 


(25) 


r 

jL. 


2 


/  F(u)  du 


/ 1  f  (t)  rdt 


(26) 


A u>  =  root  mean  square  frequency  deviation 
Au2  =  w7  -  (cr)  2  =  (u-  U)  * 
t  =  mean  duration  deviation 


=  /t]f  (t)  |2dt  =  .  /F*  (uj) 


/  f (t)  dt 


/  |  F  (u) )  |  dw 


tz  =  mean  square  duration 
_  /t2|f  jt)  J2dt  _ 


/If  (t)  pdt 


dio _ 

/  If  (cj)  I2  dco 


(27) 


(28) 


(29) 
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At  =  root  mean  square  duration 


At2  =  t2  -  (7) 2  =  (t-t)2  (30) 

St  -  -  4  -  j  /t  (31) 

/ | f (t ) |2dt 

Now  the  frequency  modulation  in  the  signal  can  be  determined 
by  calculating  the  quantity  (ait  -  ta  t)  . 

If  we  choose  ca  =  O,  t  =  0,  and  / ) f  ( t )  | 2dt  =  1  =  2E, 
then  we  can  write: 

ZJt  =  -  i  -  j  ;t  f*(t)  [^|iti]dt 

\^t\2  =  -  \  +  f/t  f*(t)  [^l^-idti2 

Applying  the  Schwarz  inequality  we  get: 

|  art  |  2  +  i  £  t2  •  u)2 

Since  cat  is  a  real  number,  we  can  define  the  uncertainty 
relation  as 

~~2  ~~2  —  2  1 

o  •  t  -  (tatr  >  j  (32) 

The  uncertainty  relation  can  also  be  defined  as 

AcjA t  >  |  (33) 

Because  of  this  relationship,  we  see  that  a  decrease  in  the 
duration  must  result  in  an  increase  in  bandwidth  and  vice-versa. 


Optimum  Receiver  in  White  Gaussian  Noise 


I 


1 


To  determine  the  optimum  receiver  in  the  presence  of 
white  Gaussian  noise,  we  consider  the  most  simple  case  - 
signal  known  exactly.  The  Bayes  detection  strategy  for  an 
exactly  known  signal  is  a  likelihood-ratio  test.  Knowledge 
of  the  form  of  the  likelihood  ratio  permits  specification  of 
a  test  statistic  which  describes  the  structure  of  the  optimum 
Bayes  detection  receiver  (Ref  1:291). 

To  derive  an  expression  for  the  likelihood  ratio,  it  is 
convenient  to  employ  the  sampl ing-theorem  waveform  represen¬ 
tation.  The  low-pass  sampling  theorem  can  be  stated  as 
follows : 

Consider  a  signal,  s(t),  bandlimited  to  (-to  ,  a)  )  .  If 

c  c 


s(t)  is  sampled  at  values  spaced 


2(j ir 


sec  apart,  then  s(t)  can 


be  completely  specified  by  those  samples  using  the  expression 


*  ,  *  ,  sin  2"£clt-2?;) 

s(t)  =  Z  s  (-^T— )  - : c— 


l2f 

k=-co  ^rC 


27T  f 

c'^  2f , 


•) 


(34) 


where  f  =  uj  /2tt  (Ref  1:53). 
c  c 

It  can  also  be  shown  that 


/ 


T 

o 


2 

s 


(t) 


2fcT 


Z 

k=l 


(35) 


where  T  is  the  signal  duration  (Ref  1:56). 

Assuming  the  signal  is  known  and  bandlimited  between 
(-fc, fc)  and  the  noise  is  also  bandlimited,  we  can  write  the 


18 


si  +  ni  w^ere 


input  waveform,  V,  where  V  =  s.  +  n,  as  = 
vi  »  V  ('Jf- )  >  si  =  s  ( 2f— )  »  and  ni  =  n  ( jf”)  ;  i  =  1,  2 _ ,  2fcT. 

C  C  C 

The  hypothesis  (signal  plus  noise)  can  be  represented  by 
the  conditional  joint  probability  density  function 


P  (V  (Hj^)  = 


<2ir)N/2(N  f  ) 
o  c 


«p  t-  ?  <Vsi)2j 

C  O  1  =  1 


(36) 


where  N  =  2fcT  and  the  power  spectral  density  of  white  noise 
is  NQ/2.  Similarly,  the  hypothesis  H0  (noise  alone)  can  be 
represented  by  the  conditional  joint  probability  density 
function 


p(V|Ho)  = 


(27r)N/2(fcNo)N 


exp  (- 


2f  N 
c  o 


N  o 

E  V  2) 
i=l 


(37) 


The  likelihood-ratio  expression  is  given  by 

1 


N 

E 


A  (V)  =  PttlHl>  =  ^P  (-  2fcNQ  (V1>si)*} 


(38) 


P(V|Hq) 


exp  {- 


2f  N  .  , 
c  o  1=1 


N  2 
Z  V.  } 

l 


Using  Eq  35,  we  can  write  Eq  38  as 

A  (V)  =  £  (V  (t) )  =  exp  {-  /*  s2  (t)  dt  +  /J'  V(t)s(t)dt) 

°  °  (39) 

Now 

In  £ (V (t) )  =  -  +  =-  /I  V (t) s (t) dt  (40) 

N0  Nq  o 

Thus,  the  likelihood  ratio  receiver  is 
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(41) 


/T  V(t)  s(t)  dt  >*  Threshold  +  j§— 

No  °  <  o 

H 

o 

The  test  statistic  (quantity  on  the  left  side  of  the  in¬ 
equality)  represents  the  cross  correlation  of  V(t)  and  s(t). 

The  matched  filter  is  used  to  realize  the  cross  correlator 
(Ref  1:292-294).  We  note  that  the  matched  filter  is  also 
the  optimum  detector  since  it  comes  from  the  likelihood  ratio 
test  using  the  Neyman-Pearson  criterion;  that  is,  for  a  fixed 
probability  of  false  alarm  (PpA=a),  the  matched  filter  is  the 
receiver  that  maximizes  the  probability  of  detection. 

Ambiguity  Function 

When  the  information  about  the  radar  environment  is 
limited,  it  appears  reasonable  to  study  the  performance  of 
matched  filter  receivers  since  they  provide  the  optimum  per¬ 
formance  in  white  Gaussian  noise.  We  now  want  to  determine 
the  resolution  capabilities  (ability  to  recognize  a  particular 
target  in  the  presence  of  others)  of  a  radar  waveform  assuming 
our  radar  uses  a  matched  filter  receiver  (Ref  5:116). 

Suppose  we  have  two  point  targets  differing  in  range.  Let 
x(t)  =  a(t)  cos  (coQt  +  <#>  ( t ) )  be  the  return  from  the  first  tar¬ 
get  and  x  ( t— t )  =  a(t-x)  cos(u)Q(t-T)  +  $>(t-x))  be  the  return 
from  the  second  target,  where  the  difference  in  range  between 
the  signals  is  Range  =  ^ ,  where  C  is  the  speed  of  light. 

From  our  discussion  of  the  complex  envelope  representation 
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increase  the  effective  bandwidth. 


Since  we  now  have  a  method  to  determine  range  resolu¬ 
tion,  we  use  the  time-frequency  duality  principle  to  formu¬ 
late  a  method  to  determine  velocity  resolution.  We  now  con¬ 
sider  two  point  targets  at  the  same  range,  but  having  different 
radial  velocities.  Again  trying  to  maximize  the  integral 
square  error  between  the  signals,  we  find 

e"^  =  /|x(t)  -  x  (t)  e^  |  2dt  (45) 


where  cuD  is  the  difference  in  the  target  Doppler-shi f t  fre¬ 
quencies  (in  radians) . 

Then 

e2  =  2/  |  x  (t)  |  2dt  -  2Re{/x(t)x*  ( t )  eJaJDtdt }  (46) 

We  want  to  minimize  /x (t) x* (t)eJWDtdt,  and  we  note  that 
/x(t)x*  (t)ej27,Vtdt  =  j'X*  (V)X(V+W)dW  =  K  (V)  where 

u>D  =  2ttV  and  X(V)  =  Fourier  Transform  of  x(t)  (47) 


Therefore,  we  want  the  spectral  correlation  function  KX(V) 
to  be  as  impulsive  as  possible  to  achieve  good  resolution. 
We  now  define  the  frequency  resolution  constant  as 


P  =  L  lu  (t)  1  4dt _ 

1  [/  i  y (t)  j  2dt ] 2 

where  X(t)  =  Re{y (t) e^Jot} . 

We  also  define  the  effective  time  duration, 


(48) 


Thus , 
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for  two  targets  at  the  same  range,  they  can  be  resolved  in 
velocity  only  if 


F 


r 


where  Awp  =  difference  in  target  Doppler  shift. 

Another  way  of  considering  target  resolution  for  a 
single  parameter  is  that  targets  must  occupy  different  bins. 

A  bin  is  the  interval  whose  width  corresponds  to  the  half¬ 
power  width  of  the  central  lobe  of  the  matched  filter  response. 
When  we  consider  this  approach,  we  assume  that  clutter  in 
other  bins  is  low  and  that  there  is  no  coupling  between  the 
two  parameters.  When  we  considered  range  resolution,  we  had 

C(t)  =  /n*(t)y  (t-r)dt  =  /|M(f)  [2ej27Tftdf 


To  define  the  bin  width,  we  say 

|  C  (To)  I  =|  =  U!M(f)  |2ej2*ftdf  |  (49) 


If  we  consider  the  Taylor  series  approximation  of  the  magni¬ 
tude  squared,  we  get 


c(t)2  =  |  =  I  C  (0) 


-C(x) 


9t‘ 


t=<5! 


(50) 


1  1  2  2 

Assuming  C(0)=1,  then  --  =  1  +  ^  (-26  )t  where  6  is  the  rms 

2  — 2 

bandwidth.  Note  that  g  =  cu  as  defined  in  Eq  26.  Then 
2  2  1  1 

6  i  =T  and  t  =  T^g.  We  must  consider  the  two-sided  nature 
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of  the  response,  and  thus. 


t'  =  2t  =  /7/e 


T  'c  c 

We  then  define  the  range  bin  as  r^  =  =  ■jpyjj-.  For  velocity 

resolution,  we  found  j  K  (V)  |  ^  =  |  /  |  u  (t)  |  2e^2lTVtdt  |  2 .  Doing 
the  Taylor  Series  approximation  about  V=0,  we  find 

2  2 

! K  (V)  j2  =  i  =  |K  (0)  I2  +  i  9  LK4V).)-.|V=0  •  V2  (52) 

y  U  sv2 


Assuming  normalized  signal  energy,  (0)  =  1,  then  we  find 

/2"  p  — 2 

Vj-,  =  — -  where  t^  is  the  rms  duration.  Note  that  t<j  =  t  as 
defined  in  Eq  29.  Therefore,  the  Doppler  bin  is  defined  as 


-  xvb  =  * 

b  2  7 ?C 


For  the  more  general  case,  we  now  consider  two  point 
targets  that  differ  in  both  range  and  velocity,  and  consider 
the  resolution  capabilities  of  a  waveform.  Again  we  want  to 
maximize  the  integral  square  error  between  the  signals,  and 
note  that 

e2  =  /IXUJe^D1  -  X  (t-x)  1 2dt 

=  2  / j  X  (t)  | 2dt  -  Re{/X(t)X*(t-x) e^wDtdt}  (54) 

We  then  want  to  minimize  the  value  /X  (t)X*  (t-T)e-'u)Dtdt.  If 
we  then  consider  the  complex  envelope  representation  of  this 
value,  we  get 
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(55) 


x(t(gjd)  =  f\i  <t)  y*  (t-r)  e-^f^dt 

which  we  will  define  as  the  ambiguity  function  of  the  wave-  * 
form.  Thus,  the  ambiguity  function  describes  the  complex 
envelope  of  the  signal  at  the  output  of  the  matched-filter 
receiver  as  a  function  of  radar  target  range  and  radial  velo¬ 
city.  This  function,  which  is  determined  solely  by  the  trans¬ 
mitted  waveform  and  matched-filter  characteristics,  has  proved 
to  be  extremely  useful  in  the  design  and  study  of  radar  sig¬ 
nals,  since  it  answers  questions  about  resolution,  ambigui¬ 
ties,  measurement  precision,  and  clutter  rejection. 

The  ambiguity  surface,  or  the  plot  of  |x(t,w)|  above 
the  t,cj  plane,  plays  a  central  part  in  the  analysis  of  reso¬ 
lution.  The  width  of  the  central  peak  of  the  ambiguity  sur¬ 
face  is  a  measure  for  close  target  resolution  in  delay  and 
Doppler,  and  subsidiary  spikes  and  low-level  parts  of  the 
surface  give  an  indication  of  the  problem  of  self-clutter  and 
target  masking  by  mutual  interference  (Ref  5:114). 

Properties  of  the  Ambiguity  Function 

The  ambiguity  function  (Eq  55)  in  the  frequency  domain 
is  given  as: 

x  ( t ,  u) )  =  /M ( f +io)  M*(f  )e"j2'ftdf  (56) 

The  maximum  value  of  the  ambiguity  function  occurs  at  the 
origin; 

|  x  (t,  w)  |  <_  x  (0, 0)  =  2E 
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Often,  the  signal  is  normalized,  and  the  maximum  value  is 
unity. 

Now  x  ( — x  ,  — co )  =  e"' WTx*  (t  ,  to)  .  Since  |  x  *  ( x ,  tu )  |  =  |x(t,u>)  |, 
we  have  |x(-x,-u>)  |  =  |x(x,gj)  |.  This  means  that  the  ambiguity 
function  is  symmetric  with  respect  to  the  origin. 

We  then  consider  the  volume  of  the  square  magnitude  of 
the  ambiguity  function: 

2 

Volume  =  // |  x  (t  ,  w)  |  dxdu> 

Now  Ft  v|x(t,V)|2  ='  //  |x(t,V)  [  2e--' 271(5  Te"-’ 2TTpVdxdV  =  |x(<S,]j)|2 

If  we  let  6=0,  y=0,  then 

// |x(T,V) |2dTdV  =  | x (0, 0) | 2  =  ( 2E) 2 

This  means  that  waveform  design  can  never  get  rid  of  the 
volume  under  |x(r,V)  |  ,*  all  proper  waveform  design  can  do  is 
to  try  to  move  the  inevitable  volume  to  regions  in  the  x,v 
plane  that  are  relatively  free  of  clutter  and  other  inter¬ 
fering  targets  (Ref  12:3-12). 

Targets  at  ranges  and  velocities  such  that  |x(t,V)|  is 
about  as  large  as  |x(0,0)|  are  indistinguishable  to  the  radar. 
The  width  of  the  peak  about  J  x  ( 0 , 0 )  |  defines  the  resolution 
of  the  waveform.  Other  peaks  away  from  |x(0,0)|  correspond 
to  ambiguities  of  the  waveform. 

Waveforms  having  x(x,V)=0  in  the  region  of  the  r,V 
plane  where  clutter  exists  generally  have  good  clutter-rejec- 


tion  properties.  More  precisely,  the  summation  over  the 
entire  t,V  plane  of  the  product  of  |x(t,V) |  with  the  clutter 
determines  the  total  interfering  clutter  signal  (Ref  12:3-10). 

Digital  Signal  Processing 

In  computing  the  ambiguity  function  of  waveforms,  a 
linear  convolution  of  samples  of  the  complex  envelope  repre¬ 
sentation  of  the  Doppler-shifted  signal  and  the  time-reversed, 
conjugate  of  the  signal  is  performed.  This  procedure  results 
in  a  cut  of  the  ambiguity  function  for  a  particular  Doppler 
frequency.  The  linear  convolution  itself  is  done  by  taking 
the  Fourier  transform  of  the  two  arguments,  multiplying  the 
transforms,  and  then  performing  the  inverse  Fourier  transform 
of  the  product.  The  Fourier  transforms  of  the  samples  are 
computed  using  a  fast  Fourier  transform  procedure.  We  shall 
now  briefly  describe  the  basic  theory  required  to  justify 
this  approach. 

To  begin,  we  must  consider  how  to  sample  the  two  argu¬ 
ments  in  order  to  accurately  reproduce  the  arguments.  This 
can  be  determined  by  the  uniform  sampling  theorem  for  lowpass 
signals  which  may  be  stated  as:  Theorem.  If  a  signal  x(t) 
contains  no  frequency  components  for  frequencies  above  f  = 

*  k 

W  Hz,  then  it  is  completely  described  by  instantaneous  sample 
values  uniformly  spaced  in  time  with  period  Tg  <  The 

signal  can  be  exactly  reconstructed  from  the  sampled  waveform 
by  passing  it  through  an  ideal  lowpass  filter  with  bandwidth 
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( 


B,  where  W<B<f  -W  with  f  =  :  The  frequency  2W  is  referred 

s  s 

to  as  the  Nyquist  frequency  (Ref  16:68).  Thus,  ,/e  know  that 
we  must  sample  the  arguments  at  intervals  less  than  or  equal 
to  the  Nyquist  rate  (^)  in  order  to  accurately  reproduce  the 
waveform.  These  samples  of  the  arguments  form  finite-duration 
sequences . 

For  a  finite  duration  sequence,  we  define  the  discrete 
Fourier  transform,  X(k),  of  that  sequence,  X(n),  as: 


N-l 


X  (k)  =  I  x  (n)  W*n,  0  <  k  <  N-l 
_  N  —  — 

n=0 

0,  otherwise 


(57) 


where 


N  is  the  number  of  samples  in  the  sequence 

=  -j  (2it/N) 
wN  -  e 


The  inverse  of  the  discrete  Fourier  transform  is  given  by: 


1  -kn 

x(n)  =  ^  Z  X  (k)  W  *n,  0  £  n  <.  N-l 
k=0 


(58) 


0,  otherwise  (Ref  4:100) 


Highly  efficient  algorithms  are  available  to  compute 
the  discrete  Fourier  transform  (DFT) ,  known  as  Fast  Fourier 
transforms,  of  a  finite-duration  sequence.  For  this  reason, 
it  is  computationally  efficient  to  consider  implementing  a 
convolution  of  two  sequences  by  computing  tneir  discrete 
Fourier  transforms,  multiplying,  and  computing  the  inverse 
discrete  Fourier  transform.  Now  consider  two  N-point 
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sequences,  (n)  and  x^n),  and  denote  their  linear  convolu¬ 
tion  as: 

N-l 

x^(n)  =  Z  x^ (m)  x2 (n-m)  (59) 

m=0 


It  is  straightforward  to  verify  that  x^  (n)  is  of  length  2N-1 . 
If  it  is  obtained  by  multiplying  the  discrete  Fourier  trans¬ 
forms  of  x^  (n)  and  x2(n),  then  each  of  these  discrete  Fourier 
transforms,  (k)  and  X^k),  must  also  have  been  computed  on 
the  basis  of  2N-1  points.  The  easiest  way  to  consider  expand 
ing  the  N-point  sequences  to  2N-1  point  sequences  is  that  N-l 
zeros  are  added  to  the  end  of  each  of  the  sequences.  Thus  if 
we  define 


x1  (k) 

i 

4 

\ 

x2  (k) 

x3  (n) 


2N-2 

Z 

n=0 

2N-2 

Z 

n=0 


_1 _ 

2N-1 


i  \  ,,nk 

X1  (n)  W2N-1 


x2(n)  w5n_x 


2N-2 

I  x. (k)X-(k)  w 
k=0  1  1 


-nk 

2N-1 


then  x^  (n)  will  be  the  linear  convolution  of  x-^  (n)  and  X2(n). 
Of  course,  a  linear  convolution  would  also  be  achieved  if  the 
discrete  Fourier  transforms  were  computed  on  the  basis  of 
more  than  2N-1  points,  but  would  not  in  general  be  achieved 
if  the  discrete  Fourier  transforms  were  computed  on  the  basis 
of  fewer  points  (Ref  4:110-111). 


There  are  several  fast  Fourier  transform  methods  avail- 


able  to  compute  the  discrete  Fourier  transform  of  a  sequence; 
however  the  one  we  shall  consider  is  called  the  decimation- 
in-time  method.  in  this  method,  the  discrete  Fourier  trans¬ 
form  is  decomposed  into  successively  smaller  DFT  computations. 
Thus,  the  original  DFT,  as  represented  by  Eq  57,  is  separated 
into  its  even  and  odd  points  such  that 

X  (k)  =  Z  x(n)  W.,nk  +  Z  x  (n)  WMnk  (60) 

n  even  n  odd 

2 

Using  the  fact  that  WN  =  WN^2'  we  can  write  the  DFT  as 

(N/2) -i  (N/2) -1 

X  (k)  =  Z  x  ( 2n)  W™  +  W*  Z  x(2n+l)  W™ 

n=0  /  n=0  ' 

=  G (k)  +  WNk  H(k)  (61) 

Then  G(k)  and  H(k)  would  be  broken  down  in  the  same  manner 

until  only  two-point  transforms  remain.  Thus,  the  number  of 

complex  multiplications  and  additions  needed  to  compute  the 

2 

DFT  would  be  reduced  from  N  to  N  log^N.  (Ref  4:290-294). 

Clutter  Models  and  Applications 

The  function  of  the  clutter  model  is  to  represent  the 
random  clutter  process  in  such  a  manner  that  system  perfor¬ 
mance  can  be  determined.  Ideally  the  clutter  model  accurately 
reflects  the  clutter  process  in  a  realistic  manner.  Practi¬ 
cally,  the  clutter  data  base  is  rarely  adequate  to  achieve 
this  objective.  The  art  of  modeling  the  clutter  involves  a 
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compromise  between  accurate  representation  and  analytic  con¬ 
venience  (Ref  9:37). 


Clutter  is  characterized  by  its  mean  or  median  equiva¬ 
lent  backscattering  cross-section,  amplitude  probability 
density  function,  and  its  power  spectral  density  or  equiva¬ 
lent  autocorrelation  function.  Angel  clutter,  radar  echoes 
primarily  caused  by  birds  and  insects,  has  velocities  that 
can  vary  from  15  to  60  knots  depending  on  wind  conditions. 

Ground  clutter  has  a  narrow  intrinsic  spectral  width  deter¬ 
mined  by  internal  clutter  motion  induced  by  the  wind.  The 
ground  clutter  spectral  density  is  usually  assumed  to  have  a 
Gaussian  shape.  Sea  clutter  power  spectral  density  is  usu¬ 
ally  assumed  to  be  Gaussian  shaped  with  an  rms  velocity  spread 
depending  upon  sea  state  (a  measure  of  the  height  of  sea  waves) . 
Sea  clutter  has  a  mean  velocity  dependent  upon  wind  speed, 
radar  polarization,  and  direction  of  the  sea.  Weather  clutter 
has  a  mean  velocity  proportional  to  the  wind  speed.  Wind 
turbulence  and  wind  shear  effect  (caused  by  the  change  in  wind 
speed  with  height)  can  cause  considerable  spectral  spreading 
(Ref  9:5-17) . 

The  spectrum  of  echoes  from  rain,  snow,  and  chaff  may 
be  deduced  from  the  motion  of  the  wind  since  these  particles 
are  good  wind  tracers.  For  the  radar  frequencies  which  we 
will  consider,  we  can  convert  scatterer  motion  to  clutter 
spectra  through  the  Doppler  equation.  The  spectrum  of  weather 
clutter  is  found  to  be  dominated  by  two  effects  which  describe 
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the  motion  of  the  wind-driven  particles:  turbulence  and  wind 
shear.  Turbulence  describes  the  fluctuations  of  the  wind 
about  its  mean  motion.  The  center  of  the  spectrum  will  be  at 
the  frequency  corresponding  to  the  mean  wind  speed  with 
respect  to  the  radar  beam  center.  Experiments  indicate  that 
the  width  of  the  turbulence  spectrum,  expressed  by  its  stan¬ 
dard  deviation,  °turb<  is  approximately  °turb  -  1.0  m/sec. 

The  wind  shear  effect  is  due  to  the  change  in  wind 
speed  with  height,  which  can  usually  be  represented  as  a 
linearly  increasing  function  over  the  elevation  extent  of  the 
beam.  The  spectrum  width  due  to  wind  shear  may  be  evaluated 
for  the  linear  wind  shear  model  by  assuming  a  Gaussian  antenna 
pattern.  As  a  result. 


a 


shear 


0.42  kR02 


(62) 


where: 

k  is  a  constant  representing  the  velocity  gradient 
component  vs  height  in  the  direction  of  the  beam 
(m/ sec/km) . 

R  is  the  slant  range  to  the  clutter  (km) . 

&2  is  the  two-way,  half-power  antenna  beamwidth  in  ele¬ 
vation  (radians) . 

Although  this  equation  seems  to  indicate  that  the  spectrum 
width  can  increase  indefinitely,  in  practice  the  value  of 

°shear  iimitec^  to  about  6  m/s ec.  For  arbitrary  radar 
azimuth,  an  average  value  of  k  is  4.0  m/sec/km.  The  _>tal 
spectrum  width,  o  ,  due  to  both  turbulence  and  wind  shear,  is: 


32 


av  =  (0turb2  +  °shear2)!s  (m/sec)  {63) 

Because  of  the  wind  shear  effect,  the  center  of  the  clutter 
spectrum  has  an  offset  value  corresponding  to  the  wind  speed 
at  the  beam  center  (Ref  2:387). 

Data  from  a  number  of  investigations  at  various  fre¬ 
quencies  show  that  the  spectrum  of  sea  clutter  in  velocity 
units  is  a  nearly  linear  function  of  wind  speed,  or  alterna¬ 
tively,  sea  state,  as  plotted  in  the  figure  below.  The  width 
of  the  Doppler  spectrum  may  be  obtained  from  this  figure 
simply  through  the  Doppler  equation.  The  literature  indi¬ 
cates  that  the  spectrum  width  is  not  very  sensitive  to  pulse 
width,  azimuth,  depression  angle,  or  polarization  (Ref  2:388). 


A 

i 


( 


Figure  1.  Spectrum  of  Sea  Clutter  (in  velocity) 
(Ref  2:388) 
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Land  clutter  extends  from  the  radar  site  out  to  the 
radar  horizon.  The  mean  velocity  of  the  land  clutter  is 
obviously  zero.  The  standard  deviation  of  the  velocity  spec¬ 
trum  of  land  clutter  is  mainly  a  function  of  the  surface  wind 
speed.  As  an  example,  the  standard  deviation  of  the  velocity 
spectrum  for  clutter  echoes  from  a  wooded  area  with  a  20  knot 
surface  wind  is  =  0.25  (m/sec)  (Ref  2:390). 

In  all  of  the  examples  given,  the  spectral  width  has 
been  expressed  in  m/sec  and  must  be  multiplied  by  2/X  to  con¬ 
vert  to  Doppler  frequency.  The  relationship  between  the 
physical  motion  of  the  clutter  scatterers  and  the  spectral 
width  has  been  found  to  conform  to  the  Doppler  equation 
throughout  most  of  the  microwave  region  of  interest  (Ref  2:390). 

When  considering  how  clutter  is  distributed  in  range, 
we  can,  in  general,  assume  that  it  is  uniformly  distributed 
in  the  specified  clutter  areas.  Therefore,  we  only  need  to 
determine  the  starting  and  ending  ranges  for  each  of  the 
clutter  areas.  The  uniform  probability  density  function  is 
given  by: 

p(T)  =  nfe y  (64) 

where 

b  is  the  ending  range  of  the  clutter  area. 

a  is  the  starting  range  of  the  clutter  area. 

c  is  the  speed  of  light. 

If  we  now  desire  to  determine  a  model  of  the  clutter 


( 


in  both  delay  and  Doppler,  we  want  a  probability  density 
function  as  a  function  of  both  variables.  In  this  paper  the 
variables  are  assumed  to  be  independent,  therefore,  we  can 
determine  the  joint  probability  function  as 


p(t,  v) 


c 

2 (b-a) 


exp 


-v-V 


(2a  2) 
v 


(65) 


where: 

a,b,c  are  as  defined  previously. 

v  is  the  mean  Doppler  frequency  of  the  clutter. 

av  is  the  standard  deviation  of  the  Doppler  frequency. 

Since  the  clutter  power  often  is  far  in  excess  of  that 
from  the  target,  it  is  necessary  to  choose  the  transmit  wave¬ 
form  on  the  basis  of  the  clutter  distributions.  The  Ambiguity 
Diagram,  a  plot  of  the  magnitude  of  the  ambiguity  function  as 
a  function  of  delay  and  Doppler,  has  been  used  as  a  tool  for 
evaluating  the  choice  of  waveform  in  specific  environments. 

It  has  proven  quite  successful  for  suggesting  specific  wave¬ 
forms  for  resolving  multiple  targets  in  range  or  range  rate 
and  in  certain  highly  specific  clutter  environments.  In  the 
extreme  case,  all  signal  waveforms  are  equally  good  (or  bad) 
as  long  as  they  are  not  compared  against  a  specific  radar 
environment.  Fortunately,  there  are  inevitably  some  bounds 
on  the  location,  velocity,  and  intensity  of  clutter  (Ref  2: 
389)  . 
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Assuming  that  we  can  now  model  the  clutter  for  the  en¬ 
vironment,  we  want  to  use  the  model  to  assist  in  selecting 
the  radar  waveform.  If  we  assume  the  transmitted  signal  to 
be 

S  (t)  =  Ay  (t)  e^u)ot 

the  returns  from  the  scatterers  will  be 

N 

S  (t)  =  Z  A  u  (t-T  )ej  (Ljo-SV  (t-Tn)  (66) 

c  n=l  n  n 

where 

An  is  the  amplitude  of  the  nth  return. 

N  is  the  number  of  scatterers. 

t n  is  the  delay  of  scatterers  from  matched  filter. 

ft  is  the  Doppler  of  scatterers  from  matched  filter, 
n 

Assuming  a  matched-filter  receiver,  the  total  output  response 
will  be: 

1  N  N  .  ,  0  ,  ,  . 

1  1  r  r  »  rt  /  ,  J  (U>  -ft  )  (Y~T  )  i 

y(t)  =  —  Z  Z  A  flu  (y-  t  )eJ  o  n  n  J 

2  m=l  n=l  n  n 

[p*(t  -t+y-t  )  e~^u)o-^m^To~t+^o-Tm^  ] dY 
o  '  m 

.  N  N  ,  ,  _  . 

1  .  .  _  —  1  ( u)  —  ft  )  T  f  .  \  /£T\ 

=  —  Z  Z  A  e  J  o  m  n  x(t  ,  <p  )  (67) 

2  .  .  n  n  n 

m=l  n=l 

Now  we  consider  the  term  e--*  Tn.  If  the  scatterers  are 

modeled  as  randomly  distributed  in  range,  then  6  =  (u>  -ft  )  t 

n  o  m  n 

is  a  random  angle.  Assuming  independence,  then  6^  is  modeled 
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by  a  uniform  probability  density  function  (pdf),  and  we  get: 


y  (t)  =  ~  Z  A_  e  ^9n  x(x  ,  4>  ) 


2  n=l  n 


n  n 


(68) 


One  measure  of  performance  is  the  peak  signal  energy  to  clut¬ 
ter  power  ratio: 


= 


at  o 

X.  1  x..(  Q ,  QJ  _l__ 


E{^iT  A  e  j8n  x(t  ,  )  |2> 

n=l  n  n  n 


By  appropriate  manipulation,  we  find: 


E{%  =  AT2  |  x  (0,0)  l2 
C 

Z  A  2  |  x  (t  ,  )  |2 

n=l  n  n  n 


(69) 


(70) 


When  the  clutter  is  dense,  it  is  useful  to  replace  the  dis¬ 
crete  scatterers  by  a  continuous  density  function  in  delay 
and  Doppler  space.  Thus,  we  can  write: 

2 

y  _  _ qT  |  x  iu ,  u;  | 

Lc 


E(#>  = 


Na//f  ( t ,  a))  !x(t,oj)  |  ^didw 


(71) 


where 


f  (  t  ,  ui )  is  the  clutter  pdf. 

aT  is  the  target  radar-cross  section. 

N  is  the  total  number  of  scatterers. 
a  is  the  average  radar-cross  section  (E(An  }=a) . 

Thus,  with  an  appropriate  probability  density  function 
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modeling  the  clutter  and  the  ambiguity  function,  the  signal- 
to-clutter  ratio  can  be  estimated  (Ref  12:3-39). 

In  order  to  approximate  the  average  radar-cross  section 
of  clutter  returns,  a  brief  discussion  of  clutter  amplitude 
models  will  be  given.  The  Gaussian  or  Rayleigh  envelope 
clutter  model  is  applicable  to  distributed  clutter  that  is 
returned  from  a  spatially  continuous  distribution  of  scatterers 
with  no  subset  of  scatterers  predominating.  This  type  of 
clutter  is  usually  associated  with  weather  clutter,  chaff, 
sea  clutter  observed  with  a  low  resolution  (pulse  width  > 

0.5  ysec)  or  with  a  high-resolution  radar  at  high  grazing 
angles  (4>>5°),  and  land  clutter  observed  from  high  grazing 
angles  over  undeveloped  terrain.  The  clutter  cross-sectional 
area  is  proportional  to  the  power  envelope  resulting  in  an 
exponential  probability  density  function  for  the  clutter 
cross-section  (C) : 

f  (C)  =  Jr  exp  [-C/C] 

where  C  is  the  mean  clutter  cross-section  and  is  given  by 

C  =  a°C'  (72) 

where  o°  is  the  mean  backscattering  coefficient. 

C'  is  the  area  intercepted  by  the  radar  cell  (Ref  9:38-39). 

The  log-normal  clutter  model  has  been  used  to  model 
high-resolution  sea-clutter  data,  where  the  sea  clutter  is 


observed  at  grazing  angles  less  than  5  degrees.  In  addition, 
ground  clutter  observed  at  low  grazing  angles  has  also  been 
modeled  using  the  log-normal  model.  The  mean  value  of  the 
clutter  cross-section  is  given  by: 

C  =  PcCm  (73) 

where 

Cm  is  the  median  value  of  the  clutter  cross-section. 

Pc  is  the  mean-to-median  ratio,  where 
Pc  =  exp  (op2/ 2) 

where 

a  =  twice  the  standard  deviation  of  the  underlying 
p  normal  distribution. 

The  distribution  function  for  the  cross-sectional  area  is 
given  by 

F  (C)  =  \  [1  +  erf  (7^-  In  §  )]  (74) 

*  pm 

The  mean  or  median  cross-sectional  area  is  determined  from 
a  knowledge  of  the  clutter's  backscattering  coefficient. 

Table  I  gives  some  reported  values  of  ap  for  various  types 
of  clutter  (Ref  9:45-47).  Other  useful  clutter  models  can  be 
found  in  Scheleher. 
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Table  I.  Log-Normal  Clutter  Parameters 
(Ref  9:47) 


Terrain/Sea  State 

Frequency 

<P 

°p 

Sea  State  2-3 

X 

4.7° 

1.382 

Sea  State  3 

KU 

1-5° 

1.440-.960 

Sea  State  4 

X 

0.24 

1.548 

Sea  State  5* 

KU 

0.50° 

1.634 

Land  Clutter 

S 

low 

3.916 

(Discrete) 

Land  Clutter 

S 

low 

1.380 

(Distributed) 

Land  Clutter* 

P-Ka 

10-70 

.728-2.584 

Rain  Clutter 

X  and  95  GHz 

- 

.680 

Radar  Angels 

- 

- 

1.352-1.620 

*Extrapolated  Values 

$  is  the  grazing  angle 

0 


Ill .  Software  Developed  to  Resolve  the  Problem 


Approach 

The  method  of  waveform  selection  discussed  in  this 
paper  requires  three  primary  products  associated  with  soft¬ 
ware  development: 

(1)  a  method  to  calculate  the  magnitude  of  a  Doppler 
cut  of  the  ambiguity  function  for  typical  radar  waveforms. 

(2)  a  method  to  produce  a  three-dimensional  display  of 
the  ambiguity  function  magnitude  as  a  function  of  delay  and 
Doppler. 

(3)  a  method  to  calculate  the  signal-to-clutter  ratio 
for  a  specified  target,  waveform,  and  clutter  environment. 

All  the  software  developed  to  perform  these  requirements  was 
done  by  the  author  with  the  exception  of  the  fast  Fourier 
transform  program  and  the  general  three-dimensional  graphics 
program.  The  software  developed  by  the  author  was  programmed 
using  the  Fortran  V  language  on  the  Data  General  Eclipse  com¬ 
puter  . 

In  this  section  of  the  paper,  the  algorithms  developed 
by  the  author  will  be  explained.  Each  of  the  three  products 
mentioned  above  will  be  discussed  in  detail.  Appendix  B  con¬ 
tains  a  system  flowchart  for  the  software.  Appendix  D  con¬ 
tains  the  program  listings  of  the  software  developed  by  the 
author . 
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Calculating  the  Magnitude  of  the  Ambiguity  Function 

As  we  have  stated  before,  the  ambiguity  function  de¬ 
scribes  the  complex  envelope  of  the  signal  at  the  output  of 
the  radar  receiver  as  a  function  of  radar  target  range  and 
radial  velocity  (Ref  12:3-3).  It  was  defined  as: 

x(t,oj)  =  f\s  (t)  u*  (t-x)  e-’a)tdt 

The  ambiguity  function  (which  is  a  correlation  function)  can 
also  be  written  as  a  convolution: 

x  ( t  ,  u) )  =  u  (t)  e^WT*u*  (-t)  (75) 

If  we  defined  the  arguments  of  the  convolution  as  finite 
sequences,  we  could  calculate  the  ambiguity  function  by  per¬ 
forming  a  linear  convolution.  To  get  the  arguments  into 
finite  sequence  representation,  we  must  appropriately  sample 
the  complex  envelope  representation  to  form  an  N-point  se¬ 
quence,  and  then  add  at  least  N-l,  zero-valued  samples  to 
the  end  of  the  sequence  to  get  it  into  a  form  to  perform  a 
linear  convolution.  The  linear  convolution  itself  can  be 
implemented  by  calculating  the  discrete  Fourier  transform  of 
the  sequences,  multiplying,  and  calculating  the  inverse  dis¬ 
crete  Fourier  transform  of  their  product.  The  discrete 
Fourier  transform  and  its  inverse  can  be  quickly  calculated 
using  the  decimation-in-time  method  of  the  fast  Fourier 
transform. 
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To  appropriately  sample  the  waveform  to  produce  the 
N-point  sequences,  the  Nyquist  rate  of  the  waveform  must  be 
determined.  The  time-bandwidth  product,  signal  duration  mul¬ 
tiplied  by  the  signal  bandwidth,  of  the  complex  envelope 
representation  provides  the  Nyquist  rate  (Ref  3:254).  If  the 
waveform  is  sampled  at  intervals  of  1/time-bandwidth  or  less, 
then  N-point  finite  sequences  are  developed.  By  adding  at 
least  N-l  zero-valued  samples  to  tt  se  sequences,  the  linear 
convolution  can  be  performed. 

The  Eclipse  computer  has  a  library  program  which  calcu¬ 
lates  the  fast  Fourier  transform  and  its  inverse  using  the 
decimation-in-time  method  in  which  the  transform  is  decomposed 
into  successively  smaller  increments.  The  sequences  which 
will  be  transformed  must  have  their  total  number  of  samples 
equal  to  some  power  of  two,  up  to  a  maximum  of  1024,  for  the 
program  to  run  correctly  (Ref  17:1-1,3-4).  Because  of  this 
requirement,  we  must  sample  the  waveform  at  intervals  of  1/N 
or  less,  where  N  is  the  smallest  value  of  two  raised  to  a 
power  which  is  greater  than  the  time-bandwidth  product.  Also, 
we  must  add  N  zero-valued  samples  to  the  sequence  to  have  the 
proper  form  for  linear  convolution.  If  we  want  good  resolu¬ 
tion  of  the  ambiguity  function,  we  may  have  to  increase  the 
number  of  samples  of  the  waveform  to  get  the  desired  results. 

As  an  example,  assume  we  want  to  calculate  the  ambiguity 
function  of  a  rectangular  pulse  which  has  a  pulse  width  of  1 
sec  and  an  amplitude  of  1.  The  time-bandwidth  product  of  the 

4  3 
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complex  envelope  would  be  1,  and  therefore,  two  samples  of 
the  waveform  should  allow  us  to  calculate  the  ambiguity  func¬ 
tion.  Let  us  assume  that  to  get  the  resolution  desired  for 
the  ambiguity  function,  we  make  eight  samples  of  the  wave¬ 
form.  The  samples  would  be  made  at  1/8  sec  intervals,  start¬ 
ing  at  1/8  sec.  We  would  also  have  add  eight  zero-valued 
samples  to  the  sequences  to  perform  the  linear  convolution. 
Hence,  the  finite  sequence  representing  would  be: 


e j u)  (1/8)  eju)(l/4)  eju(3/8) 


ejw(l/2)  e j w  ( 5/8 ) 


ejo(3/4) 


ejw(7/8) 


0, 0,  0, 0, 0,  0, 0,  0. 


The  finite  sequence  representing  u*(-t)  would  be: 

0,0, 0,0, 0,0, 0,0,  1,1, 1,1, 1,1, 1,1. 


Before  ending  this  discussion  of  forming  appropriate 
sequences,  we  must  consider  how  to  sample  waveforms  with 
repeated  pulses  or  with  a  phase-shift  coded  pulse.  The  pro¬ 
blem  to  be  considered  here  is  to  equally  sample  each  element 
of  the  waveform  to  eliminate  any  possible  bias.  Suppose,  for 
example,  that  the  waveform  has  two  pulses,  each  1  sec  in  dur¬ 
ation,  and  a  repetition  time  of  3  secs.  If  we  choose  to 
make  16  samples  of  the  waveform,  then  the  samples  will  be 
separated  by  intervals  of  (3+3)/16  sec  or  3/8  sec.  Then 
there  will  be  two  samples  made  of  the  first  pulse  and  three 
samples  made  of  the  second  pulse.  If  the  calculations  are 
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done  now,  the  linear  involution  and  ambiguity  function 
results  would  be  incorrect.  To  correct  this  situation,  the 
samples  are  "weighted"  to  insure  each  elemental  part  of  the 
waveform  is  sampled  properly  and  equally.  In  our  example, 
we  know  the  duration  of  the  waveform  and  the  number  of  sam¬ 
ples  we  wish  to  make.  By  truncating  the  value  of  the  expres¬ 
sion,  (number  of  samples/waveform  duration)  x  (pulse  width) 
to  the  integer  value,  we  now  have  the  correct  number  of  sam¬ 
ples  to  make  of  each  pulse.  In  this  example,  the  waveform 
will  be  specified  by  the  first  twelve  samples?  thus,  we 
simply  zero-fill  the  remaining  four  samples  and  continue  the 
process  as  discussed  previously. 

After  forming  the  appropriate  finite  sequences  for 
p(x)eJU,T  and  p*(-x),  the  linear  convolution  must  be  performed. 
The  library  program  to  compute  the  fast  Fourier  transform  of 
both  sequences  is  called.  Due  to  an  error  in  this  program, 
the  results  must  be  conjugated  to  provide  the  correct  results. 
This  program  also  requires  the  results  to  be  divided  by  the 
total  number  of  samples  (2N)  to  be  completely  accurate.  The 
transformed  sequences  are  then  multiplied  together  on  a  sample- 
by-sample  basis.  The  resulting  product  is  then  inverse  trans¬ 
formed  using  the  library  program.  The  results  must  again  be 
conjugated  and  divided  by  the  number  of  samples  to  get  the 
correct  ambiguity  function  values.  These  results  are  samples 
in  delay  of  a  Doppler  cut  of  the  ambiguity  function.  The 
magnitude  of  these  samples  is  then  calculated,  and  we  have 


the  desired  results. 


In  the  discussion  so  far,  we  have  not  considered  how  to 

handle  waveforms  with  Gaussian-shaped  pulses  nor  those  having 

some  form  of  frequency  modulation.  In  order  to  normalize  the 

magnitude  of  the  ambiguity  function  when  Gaussian-shaped 

pulses  are  used,  we  must  determine  a  weighting  factor.  This 

weighting  factor  was  determined  on  the  basis  that  the  mean  of 

a  Gaussian  pulse  was  ^ (pulse  width)  and  the  standard  deviation 

was  1/6 (pulse  width).  Using  +  3  standard  deviations  from  the 

mean,  insures  that  our  approximation  is  within  2%  of  the 

actual  result.  Thus  we  want  to  determine  the  weighting  fac- 

2 

tor,  A  such  that  A/s  (t)dt=l.  Now 

s  (t)  =  '^2  exp  {- (t-t)  2/2a2}  (76) 

where 

o=~  (pulse  width) . 
o 

t  =  ^  (pulse  width) . 

Thus  we  want 

A/P  r  -~2  exp  { - (t-t) 2/o2} dt  =  1  (77) 

O  ZTTO 

Substituting  x  =  (t-t) /a  into  Eq  77,  we  get: 


A 

2^0 


/ 


A  / 


3 

-3 


(PW-t) /o 
-t/o 

.  2 
exp (-x 


exp 

)  dx 


(-x2)dx  =  1 

o  (PW) 

=  3 
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Using  a  standard  normal  distribution  table,  we  can  determine 
3  2 

the  value  of  exp  (-x  )  dx  =  0.9974/tT.  Therefore, 

A  =  [  7T  (PW)  ]/ [  3  (0.9974)  /rf]  =  0.59236  (PW)  (78) 

Thus,  to  normalize  the  magnitude  of  the  ambiguity  function 
samples,  simply  multiply  each  sample  by  the  weighting  factor 
A. 

When  the  waveform  has  been  frequency  modulated,  the 
effects  on  the  ambiguity  function  are  significant.  We  recall 
that  the  complex  envelope  representation,  p (t)  was  also  de¬ 
fined  as  y  (t)  =  |  u  (t)  J  e-1  ^  ^  ,  where  $  (t)  will  contain  the 
frequency  modulation.  In  this  study  only  linear  frequency 
modulation  (FM)  and  V-FM  techniques  were  considered.  For 
linear  FM,  the  complex  envelope  has  the  form:  p ( t )  =  | p  ( t )  | 

•j  7T]tt2 

e  .  The  instantaneous  frequency  of  the  complex  envelope 

is  fi  =  =  kt  with  k  determining  the  slope  of  the  FM. 

If  the  instantaneous  frequency  is  swept  over  the  band,  B, 
during  the  signal  duration,  T,  the  absolute  value  of  k  is 
given  as  jk|  =  B/T.  The  sign  of  k  determines  the  direction 
of  the  sweep,  with  k  positive  for  a  frequency  upsweep  and 
negative  for  a  downsweep  (Ref  5:169-170). 

Meanwhile,  if  V-FM  is  chosen,  the  waveform  is  split 
into  two  parts,  one  for  the  downsweep  and  the  other  for  the 
upsweep.  The  complex  envelope  for  V-FM  has  the  form: 


< 
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(t)  =  (u(t)-u(t-PW/2)]  |p(t)  |e"j7Tkt< 


+  [u (t-PW/2)  -  u(t-PW) ] | y (t)  | e 


j:rkt 


(79) 


where  u(t)  is  the  unit  step  function.  Again,  assuming  the 
instantaneous  frequency  is  swept  over  the  band  B  during  the 
signal  duration  T,  then  the  absolute  value  of  k  is  given  as 
| k |  =  B/T  (Ref  5:176) . 

As  an  example,  suppose  we  again  choose  our  waveform  to 
have  a  pulse  width  of  1  sec,  an  amplitude  of  1,  and  been 
sampled  every  1/8  sec  (i.e.  .8  samples  are  made  of  the  wave¬ 
form)  .  Suppose  we  also  want  to  FM  the  waveform  and  sweep 
over  a  band  of  25  Hz.  If  we  choose  linear  FM,  the  sequence 
representing  u(T)eJWT  will  be: 

j  (w/8+25it/64)  J  (u)/4+25n  (4)/64)  j  (w+25ir) 

0,  0,  0,  0,  0,  0,0,0. 


The  sequence  representing  u*(-t)  will  be: 


0.0.0.0.0.0.0.0,  e-J(25’/6'1>.  e-J(25’"4>/64) 


»  •  »  •  » 


-j  (25tt) 


When  a  waveform  has  phase-shift  coding,  the  amplitude 
of  the  complex  envelope  is  modulated.  For  example,  suppose 
we  choose  our  waveform  to  have  a  pulse  width  of  1  sec  and  an 
amplitude  of  1,  and  to  have  8  samples  made  of  it.  If  we 
assume  it  is  phase-shift  coded  as  a  5-element  Barker  code. 
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the  phase  for  each  element  would  be  0,0,0,Tr,0  (i.e.  the  ampli 

tude  of  the  waveform  for  each  element  would  be  1,1,1, -1,1). 

,  .  Jut  . , ,  , 

Then  the  sequence  representing  U(T)e  will  be: 

e  (1/5)  ej^>(2/5)<  e  juJ  (3/5 )  _ej  ^  (4/5 )  eJ<^ 

0,  0,  0,  0,  0,  0,  0, 0,  0, 0, 0 

since  we  also  had  to  weight  our  sampling  of  the  waveform. 

The  sequence  representing  u*(-t)  will  be:  0, 0, 0, 0, 0, 0, 0, 0, 0, 

0,0, 1,1,1, -1,1. 

Thus,  we  have  described  the  general  principles  used  to 
calculate  the  magnitude  of  the  ambiguity  function  for  the 
waveforms  considered.  The  software  developed  to  calculate 
the  ambiguity  function  magnitudes  implements  these  principles 
Magnitudes  calculated  using  this  software  were  consistently 
within  one  percent  of  the  actual  results  as  long  as  adequate 
sampling  was  made  of  the  waveform. 

Producing  Three-Dimensional  Plots  of  the  Magnitude 

One  of  the  most  useful  tools  in  selecting  the  waveform 
is  the  plot  of  the  ambiguity  function  magnitude.  We  have 
already  discussed  the  method  used  to  calculate  the  magnitude 
of  samples  made  in  delay  for  a  Doppler  cut  of  the  ambiguity 
function.  A  method  is  needed  to  select  and  properly  format 
these  magnitude  samples  (made  for  various  Doppler  cuts  of 
the  ambiguity  function)  for  use  by  the  three-dimensional 
graphics  program.  This  graphics  program,  called  PLTTRNS, 


( 
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was  written  by  R.W.  Schafer,  July  3,  1978,  at  the  Georgia 
Institute  of  Technology.  The  PLTTRNS  program  requires  the 
input  to  be  a  file  of  a  two-dimensional  array  of  complex 
numbers  in  binary  format.  The  array  is  assumed  square  and 
of  minimum  size  64  x  64  and  of  maximum  size  512  x  512.  The 
program  provides  three  possible  display  modes  on  the  Textronics 
storage  scope:  contour,  perspective,  and  single-line.  The 
contour  plot  and  three-dimensional  perspective  plot  can  only 
display  an  array  of  data  of  si2e  65  x  65.  if  the  input  array 
is  larger  than  64  x  64,  an  appropriate  number  of  rows  and 
columns  are  decimated  prior  to  display.  Therefore,  the  data 
we  want  displayed  must  be  of  size  64  x  64. 

The  problem  we  now  face  is  how  to  select  the  data  we 
want  displayed  and  how  to  properly  format  those  data  for  the 
PLTTRNS  program.  Since  our  two-dimensional  array  contains 
64  delay  samples  for  each  of  64  different  Doppler  cuts  of  the 
ambiguity  function,  we  must  decide  which  delay  samples  and 
Doppler  frequencies  to  choose.  The  software  developed  in 
this  study  allows  the  user  to  choose  which  delay  sample  to 
begin  the  display  and  the  option  to  skip  samples.  For  example, 
if  128  delay  samples  are  available  over  the  delay  times  of  -1 
sec  to  1  sec,  the  user  could  choose  to  start  at  sample  64 
(corresponding  to  time  0  sec)  and  observe  every  sample  out  to 
sample  127  (corresponding  to  time  63/64  sec).  The  user  could 
also  choose  to  observe  every  other  sample,  starting  at  sample 
1  and  ending  at  sample  127,  to  have  a  display  covering  the 
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entire  delay  interval.  This  choice  is  simply  a  matter  of  the 
resolution  desired  by  the  user. 

Although  the  choice  of  the  Doppler  frequencies  is  again 
a  matter  of  the  resolution  desired  by  the  user,  greater  care 
should  be  taken  in  the  selection  process.  Prior  knowledge  of 
the  resolution  required  for  similar  waveforms  can  assist  in 
this  selection.  With  no  prior  knowledge,  the  user  should 
begin  by  observing  closely  spaced  (interval  of  1  Hz)  cuts  and 
then  increasing  the  spacing  until  the  resolution  meets  the 
user's  desires.  As  with  the  delay  samples,  the  user  chooses 
the  starting  value  of  the  Doppler  cuts.  The  user  also  speci¬ 
fies  the  interval  between  the  Doppler  cuts,  and  the  64  cuts 
are  then  calculated  as  explained  previously. 

After  the  magnitudes  of  the  samples  of  the  ambiguity 
function  cut  have  been  calculated,  the  magnitudes  are  then 
quantized  to  a  value  between  0  and  200.  The  quantization 
occurs  because  the  PLTTRNS  program  must  have  the  data  in 
integer  form.  The  200  quantization  levels  provide  the  reso¬ 
lution  necessary  for  smoothing  of  the  curves  of  the  three- 
dimensional  graphs.  After  quantizing  the  magnitude  samples, 
they  are  transformed  into  complex  data  with  the  imaginary 
part  set  to  zero,  and  then  stored  into  a  two-dimensional  com¬ 
plex  array.  When  all  64  Doppler  cuts  have  been  calculated 
and  stored  in  the  complex  array,  the  data  are  written  onto  a 
file  called  "Data".  The  PLTTRNS  program  is  then  called,  and 
the  three-dimensional  plot  of  the  magnitude,  as  well  as  a 


contour  plot,  can  then  be  made. 

Calculating  the  Signal-to-Clutter  Ratio 

As  stated  previously,  one  of  the  constraints  in  select¬ 
ing  a  waveform  is  how  well  the  waveform  is  able  to  detect 
desirable  targets  in  the  environment.  The  signal-to-clutter 
ratio  is  the  value  which  will  be  used  to  indicate  the  perfor¬ 
mance  of  the  waveform  in  the  environment.  As  derived  pre¬ 
viously,  the  signal-to-clutter  ratio  is  given. by  Eq  71.  in 
the  theoretical  discussion  of  this  paper,  guidance  was  given 
on  how  to  establish  the  probability  density  function  of  clut¬ 
ter  returns  as  a  function  of  delay  and  Doppler  and  on  how  to 
establish  the  average  radar  cross-section  of  the  clutter 
returns.  The  user  must  establish  the  radar  cross  section  of 
the  desired  target  and  estimate  the  number  of  scatterers 

reflecting  energy.  Since  the  magnitudes  of  the  ambiguity 

2 

function  are  always  normalized  in  this  study,  |x(0,0)  |  =1. 

The  significant  task  left  for  the  software  is  to  calcu¬ 
late  the  volume  of  the  product  of  the  magnitude  squared  of 
the  ambiguity  function  and  the  clutter  probability  density 
function.  In  this  procedure  the  boundaries  of  tne  clutter 
areas  are  defined  in  terms  of  delay  and  Doppler,  as  well  as 
the  boundaries  of  the  ambiguity  function.  The  interval 
between  the  Doppler  cuts  of  the  ambiguity  function  must  be 
determined  and  specified.  It  is  important  to  note  that  no 
limit  exists  in  this  software  for  the  number  of  Doppler  cuts 


that  can  be  made.  The  interval  should  be  chosen  to  insure 
that  all  significant  changes  in  the  ambiguity  function  magni¬ 
tudes  are  considered.  The  necessary  resolution  in  delay 
should  also  dictate  the  number  of  samples  to  be  made  of  the 
waveform.  The  three-dimensional  plots  of  the  ambiguity  func¬ 
tion  magnitude  can  greatly  assist  in  these  parameter  selections. 

After  the  necessary  parameters  are  entered,  the  software 
begins  its  work.  The  algorithm  begins  by  calculating  the 
magnitudes  of  the  samples  for  the  first  two  Doppler  cuts  and 
storing  those  data  in  an  array.  Then  the  first  two  magnitude 
samples  of  both  Doppler  cuts  are  considered.  The  value  of  the 
clutter  probability  density  function  is  then  calculated  for 
each  of  the  four  points.  If  the  user  wishes  to  input  his  own 
clutter  distribution,  he  will  be  prompted  to  do  so  if  any  of 
the  points  is  included  in  a  clutter  area.  The  product  of  the 
magnitude  squared  and  the  clutter  probability  density  function 
is  then  calculated  for  each  of  the  four  points.  The  average 
value  of  the  four  products  is  then  calculated,  and  the  minimum 
and  maximum  products  are  identified.  These  three  products 
are  each  then  multiplied  by  the  incremental  area  of  the  t,u 
plane  to  produce  an  incremental  value  of  the  minimum,  maximum, 
and  average  volume.  The  incremental  area  of  the  t,oj  plane  is 
simply  the  product  of  the  interval  between  the  samples  in 
delay  and  the  interval  between  the  Doppler  cuts.  Also,  the 
incremental  volumes  are  added  to  a  cumulative  sum  for  the 
respective  minimum,  maximum,  and  average  volumes. 
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i  After  this  process  is  completed,  the  second  and  third 

magnitude  samples  for  both  Doppler  cuts  are  considered. 

Again,  a  product  of  the  magnitude  squared  and  the  value  of 
the  clutter  probability  density  function  for  each  of  these 
four  points  is  calculated.  Incremental  volumes  for  the  mini¬ 
mum,  maximum,  and  average  products  are  calculated  and  added 
to  the  respective  cumulative  sums. 

This  process  of  considering  the  next  magnitude  samples 
for  the  first  two  Doppler  cuts  is  continued  until  all  the 
samples  of  those  cuts  are  considered.  The  final  incremental 
values  will  be  those  for  the  last  and  next-to-last  samples. 
Then  we  simply  consider  the  second  and  third  Doppler  cuts 
and  repeat  the  entire  process  again.  The  entire  procedure 
ends  with  the  last  and  next-to-last  Doppler  cuts  repeating 
the  process. 

As  an  elementary  example,  suppose  we  have  three  magni¬ 
tude  samples  for  each  of  three  Doppler  cuts  of  an  ambiguity 
function.  Suppose  the  Doppler  cuts  are  made  at  frequencies 
of  1,0,  and  -1  Hz  and  the  magnitude  samples  occur  at  delay 
values  of  1,0,  and  -1  sec.  Also,  suppose  the  magnitude  of 
the  samples  and  the  values  of  the  probability  density  func¬ 
tion  for  the  specified  points  are  given  in  the  following 
table : 
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Table  II.  Example  Values  for  Calculating 
the  Signal-to-Clutter  Ratio 


Doppler  (Hz) 

Delay  (sec) 

Magnitude  of 
Amb  Fen 

Value  of 
Clutter  PDF 

-1 

-1 

0.1 

0.1 

-1 

0 

0.2 

0.1 

-1 

1 

0.1 

0.1 

0 

-1 

0.2 

0.1 

0 

0 

1.0 

0.2 

0 

1 

0.2 

0.1 

1 

-1 

0.1 

0.1 

1 

0 

0.2 

0.1 

1 

-1 

0.1 

0.1 

Considering  the  first  two  Doppler  cuts  (-1  and  0  Hz)  and  the 
first  two  delay  samples  (-1  and  0  sec)  for  each,  we  find  the 
products  of  the  magnitude  squared  and  clutter  pdf  to  be: 
0.001,  0.004,  0.004,  and  0.200.  The  average  value  would  be 
0.05225,  the  minimum  would  be  0.001,  and  the  maximum  would 
be  0.20 0.  The  incremental  volumes  would  be  the  same  since 
the  incremental  area  is  1.  Considering  the  next  samples  in 
delay  (0  and  1  sec) ,  we  find  the  incremental  volumes  to  be 
0.001,  0.200,  and  0.05225  for  the  minimum,  maximum,  and 
average  volumes,  respectively.  Then  we  consider  the  next 
Doppler  cuts  (0  and  1  Hz)  and  find  the  incremental  volumes 
for  the  first  two  delay  samples  (-1  and  0  sec)  to  be:  0.001, 
0.200,  and  0.05225  again.  Similarly,  we  find  the  same  in¬ 
cremental  volumes  for  the  iast  two  delay  samples  (0  and  1 
sec).  The  cumulative  volumes  would  then  be  0.004,  0.800, 
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and  0.209  for  the  mini;;. am ,  maximum,  and  average  volumes, 
respectively. 

After  the  volumes  have  been  calculated,  a  comparison 
can  be  made  to  determine  if  the  volumes  have  a  wide  spread. 

As  long  as  the  spread  between  the  minimum  and  maximum  volumes 
is  less  than  20  percent  of  the  maximum  volume,  we  can  assume 
the  average  volume  to  be  an  adequate  approximation  of  the 
true  volume.  If  the  spread  is  greater  than  20  percent,  the 
program  should  be  run  again  with  smaller  intervals  between 
the  samples  and  the  Doppler  cuts. 

Once  we  have  an  acceptable  approximation  of  the  volume 
of  the  magnitude  squared  of  the  ambiguity  function  and  the 
clutter  probability  density  function,  we  then  make  the  calcu¬ 
lation  of  the  signal-to-clutter  ratio.  The  cumulative  average 
volume  is  the  value  used  in  this  calculation.  Several  signal- 
tm-clutter  ratios  may  need  to  be  calculated  for  differing  tar¬ 
gets  and  types  of  clutter. 

We  note  that  a  software  program  was  also  developed  which 
generates  a  three  dimensional  plot  of  the  probability  density 
function  of  the  clutter  model  as  a  function  of  delay  and 
Doppler.  This  plot,  used  with  a  corresponding  plot  of  the 
magnitude  of  the  ambiguity  function,  can  give  the  user  a  rough 
approximation  of  the  signal-to-clutter  ratio  that  exists  for 
the  particular  waveforms  being  considered. 
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IV.  Applications  of  Software  in  Waveform  Selection 


General  Approach 

As  stated  previously,  the  purpose  of  this  study  was  to 
develop  a  method  of  selecting  a  waveform  given  the  desired 
radar  coverage  and  operational  environment.  Software  programs 
have  been  developed  in  this  study  to  provide  products  to  assist 
the  user  in  the  selection  process.  In  this  section  of  the 
study,  an  explanation  will  be  given  of  how  to  use  these  soft¬ 
ware  products  in  selecting  a  waveform.  This  section  begins 
with  a  discussion  of  how  to  analyze  the  given  requirements  of 
radar  coverage  and  the  environment.  The  next  topic  discussed 
is  a  method  of  determining  the  capabilities  of  a  waveform, 
given  the  coverage  requirements  and  environment.  A  discussion 
of  how  to  make  tradeoffs  between  the  desired  requirements  and 
actual  capabilities  available  is  then  given.  The  section  con¬ 
cludes  with  the  analysis  of  an  example  problem. 

Analyzing  Requirements 

In  this  section  we  will  briefly  explain  how  to  take  a 
standard  set  of  radar  coverage  requirements  and  environmental 
location  and  determine  what  factors  to  use  with  the  software 
products  available.  We  assume  the  standard  requirements  to 
be  given  as  the  unambiguous  range,  range  resolution,  accept¬ 
able  blind  speeds,  and  velocity  resolution  of  desired  targets. 
The  general  type  of  location,  in  which  the  radar  is  to  be 
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operating,  will  also  be  given. 


I 


The  unambiguous  range,  which  is  the  range  between  range 
bins,  can  be  used  to  determine  the  minimum  delay  required 
between  range  bins  by  the  equation: 

At  =  (80) 

where : 

UR  is  the  unambiguous  range. 

At  is  the  delay  between  delay  bins. 

C  is  the  velocity  of  light. 

The  range  resolution  requirement  specifies  the  maximum  range 
bin  width  and  can  be  used  to  determine  the  minimum  effective 
signal  bandwidth  by  the  equation: 


where: 


(81) 


8  is  the  minimum  effective  signal  bandwidth  (see  Eq  51) . 
RBW  is  the  range  bin  width. 


The  blind  speed  requirement,  which  is  the  velocity 
between  Doppler  bins,  can  be  used  to  determine  the  Doppler 
interval  required  between  bins  by  the  equation: 


Af  =  (82) 

where: 

Af  is  the  Doppler  interval. 

V  is  the  blind  speed  velocity. 

X  is  the  wavelength. 
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The  velocity  resolution  requirement  specifies  the  maximum 
Doppler  bin  width  and  can  be  used  to  determine  the  minimum 
effective  signal  duration  by  the  equation: 

fcd  =  7?  VBW  (83) 

where: 

t^  is  the  minimum  effective  signal  duration  (see  Eq  53) . 
X  is  the  wavelength. 

VBW  is  the  velocity  bin  width. 

With  the  minimum  effective  signal  bandwidth  and  minimum  effec¬ 
tive  signal  duration,  the  minimum  effective  time-bandwidth 
product  of  the  signal  can  be  determined. 

Since  the  general  type  of  environment  is  also  given, 
the  types  of  associated  clutter  can  be  deduced.  With  an  idea 
of  the  clutter,  an  analysis  can  be  performed,  as  outlined  in 
the  theoretical  discussion,  to  estimate  the  average  radar 
cross-section  of  the  clutter  and  to  develop  a  probability 
density  function  of  the  clutter  as  a  function  of  delay  and 
Doppler.  These  data  will  be  required  to  calculate  the  signal- 
to-clutter  ratio. 

Determining  Waveform  Capabilities 

By  analyzing  the  given  requirements,  the  user  can  begin 
to  use  the  software  products  in  the  selection  process.  At 
this  point  in  the  process,  the  user  knows  the  minimum  time- 
bandwidth  product,  the  minimum  spacing  between  delay  bins 
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and  between  Doppler  bins,  and  the  maximum  delay  bin  width 
and  Doppler  bin  width.  Also,  a  minimum  signal-to-clutter 
ratio  value  of  30,  which  is  an  arbitrary  estimate,  must  be 
achieved  to  detect  a  target. 

Radar  waveforms  can  be  classified  into  groups  according 
to  the  general  properties  of  their  ambiguity  function  magni¬ 
tudes.  We  can  assume  the  waveforms  fall  into  one  of  four 
categories:  constant-carrier  pulse,  irregular  or  noiselike 
pulse  (e.g.,  phase-shift  coded  pulse),  linear  frequency  modu¬ 
lated  pulse,  and  repetitive  pulses  (Ref  6:1079). 

The  basic  radar  waveform  is  the  constant-carrier  pulse, 
which  has  a  time-bandwidth  product  of  unity.  The  magnitude 
of  the  ambiguity  function,  thus,  is  concentrated  in  a  unity 
size  cell  in  delay-Doppler  space.  Depending  on  the  exact 
pulse  shape,  there  may  be  sidelobes  in  delay  and  Doppler, 
but  these  will  contain  only  an  insignificant  part  of  the 
total  energy  response.  This  type  of  waveform  will  be  desig¬ 
nated  as  a  Class  A  waveform  (Ref  6:1079). 

The  irregular  or  noiselike  pulse,  like  a  phase-shift 
coded  pulse,  has  an  ambiguity  function  magnitude  resembling 
an  inverted  thumbtack.  The  dimensions  of  the  central  spike 
of  the  thumbtack  function  are  (1/bandwidth)  x  (1/signal  dura¬ 
tion)  ,  which  define  the  size  of  the  resolution  cell.  This 
type  of  waveform  trades  a  small  resolution  cell  against  a 
potentially  severe  self-clutter  problem.  The  resolution  per¬ 
formance  will  be  poor  in  a  dense-target  environment  or  when 
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a  small  target  is  to  be  detected  in  the  presence  of  much 
stronger  reflectors.  This  waveform  type  will  be  designated 
as  a  Class  B1  waveform  (Ref  6:1079). 

When  the  instantaneous  frequency  or  the  carrier  fre¬ 
quency  is  linearly  varied  or  stepped,  the  waveform  is  classi¬ 
fied  as  a  linear  FM  signal.  The  magnitude  of  the  ambiguity 
function  of  this  signal  can  be  viewed  as  a  sheared  version  of 
the  ambiguity  function  magnitude  of  the  constant-carrier 
pulse.  It  follows  that  the  size  of  the  resolution  cell  must 
also  be  unity,  and  that  little  sidelobe  energy  need  be  out¬ 
side  the  main  response  ridge.  This  implies  an  excellent 
resolution  potential  when  targets,  in  delay-Doppler  space, 
do  not  fall  along  the  ridge,  but  essentially  no  resolution 
performance  for  targets  that  do  fall  along  the  ridge.  The 
amount  of  shearing  of  the  ridge  is  determined  by  the  bandwidth 
swept  by  the  FM;  hence,  the  larger  the  bandwidth  swept,  the 
greater  the  shearing  of  the  ridge.  Signals  of  this  type  are 
designated  as  Class  B2  waveform  (Ref  6:1080). 

If  the  waveform  is  a  periodic  repetition  of  a  pulse  whose 
bandwidth  is  large  compared  to  the  repetition  frequency,  the 
signal  is  called  a  repetitive  pulse  signal  and  is  designated 
as  a  Class  C  waveform.  The  magnitude  of  the  ambiguity  func¬ 
tion  consists  of  a  regular  array  of  spikes  in  the  delay- 
Doppler  domain,  with  the  delay  spacing  of  the  spikes  equal  to 
the  repetition  period  and  the  Doppler  spacing  equal  to  the 
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repetition  frequency.  The  dimensions  of  the  ambiguous  spikes 
are  (1/bandwidth)  x  (1/pulse  width)  as  in  the  case  of  the 
thumbtack  function.  However,  instead  of  the  uniform  sidelobe 
pedestal,  the  response  energy  is  concentrated  in  isolated 
spikes.  This  waveform  type  offers  excellent  resolution  per¬ 
formance  in  a  dense-target  environment  at  the  cost  of  a  pos¬ 
sible  ambiguity  problem  (Ref  6:1080).  The  following  table 
summarizes  the  characteristics  of  the  magnitude  of  the  ambi¬ 
guity  functions  for  the  waveform  classifications. 


Table  III.  Waveform  Classification  and 
Ambiguity  Functions 


Class  A 

Class  B1 

Class  B2 

Class  C 

Time-bandwidth 

product 

unity 

large  compared  with  unity 

ridge 

thumbtack 

sheared  ridge 

bed-of-nails 

Resolution 
cell  size 

unity 

1/TB 

unity 

1/TB 

Ambiguities 

no 

no 

range-Doppl er 
coupling 

spikes 

Sidelobes 

low 

high 

low 

low 

(Ref  6:1080) 

A  three-dimensional  diagram  and  a  contour  diagram  of 
each  waveform  classification  is  contained  in  Appendix  C. 
These  will  clearly  demonstrate  the  characteristics  which  we 


have  discussed. 


With  the  information  about  the  waveform  requirements 
which  is  now  available,  the  user  should  be  able  to  eliminate 
certain  waveform  classifications  from  further  consideration. 
Other  waveform  types  should  also  be  recognized  as  the  most 
likely  candidates  to  be  selected.  Because  high-performance 
radars  are  very  expensive,  it  is  important  that  the  user 
choose  the  optimum  waveform  in  the  sense  that  it  will  do  the 
job  in  the  least  complicated  manner.  Thus,  the  user  will 
have  a  prioritized  list  of  the  remaining  waveform  candidates 
to  analyze  further. 

Applying  Software  Products 

With  the  information  now  available,  the  user  can  fur¬ 
ther  analyze  the  waveform  by  using  the  software  products 
developed  in  this  study.  The  user  should  begin  this  analysis 
phase  with  the  least  complicated  waveform  candidate.  Since 
the  analysis  of  the  delay  and  Doppler  bin  widths  and  bin 
separations  were  only  rough  approximations  in  the  preceding 
sections,  the  user  will  want  to  obtain  a  three-dimensional 
diagram  and  a  contour  plot  (indicating  the  half-power  points) 
of  the  magnitude  of  the  ambiguity  function.  With  this  infor¬ 
mation,  the  user  can  more  precisely  establish  the  resolution 
capabilities  and  ambiguities  associated  with  the  waveform. 

If  these  resolution  capabilities  and  ambiguities  do  not  meet 
the  initial  requirements,  the  user  should  attempt  to  modify 
the  waveform  parameters  until  the  requirements  are  met  or 
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until  the  user  concludes  that,  in  fact,  the  waveform  can  not 
meet  the  requirements. 

Assuming  that  the  waveform  now  meets  the  resolution 
and  ambiguity  requirements,  the  user  must  insure  that  desired 
targets  can  be  detected  in  the  modeled  clutter  environment. 

To  determine  this,  the  user  should  use  the  software  to  calcu¬ 
late  the  signal- to-clutter  ratio.  If  the  signal-to-clutter 
ratio  is  30  or  above,  then  the  user  may  assume  the  target 
will  be  detected.  If  the  result  is  less  than  30,  the  user 
should  eliminate  the  waveform  from  further  consider. 

This  process  of  establishing  the  resolution  capabili¬ 
ties,  ambiguities,  and  signal-to-clutter  ratio  should  be  con¬ 
tinued  for  all  candidate  waveforms.  If  more  than  one  waveform 
fulfills  the  requirements,  the  user  may  want  to  consider 
tradeoffs  between  signal  complexity  to  gain  improved  waveform 
performance.  If  none  of  the  candidate  waveforms  can  fulfill 
the  original  requirements,  the  user  must  make  tradeoffs  in 
the  requirements  to  obtain  a  waveform  which  would  be  the  most 
acceptable. 

Tradeoffs 

For  the  case  of  several  waveforms  fulfilling  the  re¬ 
quirements,  the  user  would  normally  use  the  least  complex 
waveform.  If  a  more  complex  waveform  provides  significantly 
improved  performance,  the  user  might  accept  the  more  complex 
waveform,  especially  if  the  cost  to  implement  the  more  complex 
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waveform  is  within  budget  constraints.  The  user  must  deter¬ 
mine  whether  the  improved  performance  is  great  enough  to 
outweigh  the  increased  cost. 

On  the  other  hand,  if  none  of  the  candidate  waveforms 
fulfill  the  specified  requirements,  the  user  must  determine 
which  requirements  can  be  changed  and  by  how  much.  Perhaps 
the  least  flexible  requirement  is  the  signal-to-clutter  ratio 
required  to  detect  desired  targets.  The  tradeoffs  in  resolu¬ 
tions  and  ambiguities  must  be  determined  by  the  user.  By 
prioritizing  the  requirements,  the  user  can  revise  the  re¬ 
quirements  until  a  waveform  is  found  that  fulfills  them.  If 
the  waveform  can  be  implemented  within  the  cost  constraints, 
it  should  then  be  selected.  If  the  implementation  will  cost 
too  much,  the  requirements  must  continue  to  be  revised.  Thus, 
the  user's  judgement  is  the  most  important  factor  in  deter¬ 
mining  tradeoffs. 

Example  Problem 

To  make  the  preceding  discussion  of  waveform  selection 
more  clear,  an  example  problem  of  waveform  selection  will  now 
be  presented.  We  shall  assume  the  following  requirements 
have  been  given  to  us  by  the  contractor: 

Minimum  unambiguous  range  -  25  km 
Range  resolution  -  . 5  km 

Minimum  first  blind  speed  -  200  km/hr  =  556  m/s 
Velocity  resolution  -  20  km/hr  =  55.6  m/s 
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Operating  environment  -  coastal  radar 

2 

Minimum  desired  target  size  -  4m 
Cost  constraints  are  inflexible. 

Operating  frequency  of  radar  -  3  GHz 

We  begin  with  a  cursory  analysis  to  determine  the  mini¬ 
mum  time-bandwidth  product  needed  by  the  waveform.  The  mini¬ 
mum  effective  bandwidth  is  firs!  estimated  by  Eq  81: 

6  =  3  -X  ,.10,.8_m/s  —  =  4  24  x  1q5  Hz 
Si  (.5  x  10  m) 


The  minimum  effective  signal  duration  is  estimated  by  Eq  83: 


fcd  - 


3  x  108  m/s 


/7f(VBW)  /J  (3  x  109Hz)  (55.6  m/s) 


=  1.27  x  10“3s 


Hence,  the  minimum  effective  time-bandwidth  product  must  be 
approximately  538.  Therefore,  we  know  that  the  simple  pulse. 
Class  A,  can  not  be  used. 

We  now  can  approximate  the  delay  and  Doppler  bin  widths. 
The  delay  bin  width  is  approximated  as: 

Bin  width  =  2-iRan,ae-r_esolution)  =  2j.50^m,L  ,  3>33  ysec 

C  3x10  m/s 

The  Doppler  bin  width  is  approximated  as: 
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Doppler  bin  width  =  ^Velocity  Resolution) 

A 

=  2(55,6  in/sliMo!  Hz)  =  x  ^  R  Hz 
3  x  10®  m/s 

The  delay  bins  must  be  separated  by  approximately  (using 
Eq  80)  : 

At  =  — 25xl-9®m .  =  4a?  x  1q-5  sec  =  41>?  ysec 
2(3x10  m/s) 

The  Doppler  bins  must  be  separated  by  approximately  (using 
Eq  82) : 

Af  .  2<5'56  m/s)  <3X109bZ)  ,  na  k  Hz 
3  x  10°  m/s 

Now  since  we  normally  would  use  phase-shift  coding  as 
the  Class  B1  waveform,  we  would  find  it  impossible  to  achieve 
the  required  time-bandwidth  product  using  a  Class  B1  waveform. 
More  sophisticated  noise-like  coding  could  produce  the  re¬ 
quired  time-bandwidth  product,  but  the  cost  of  implementing 
it  would  be  prohibitive.  Therefore,  we  want  to  consider  the 
two  remaining  waveform  types,  linear  FM  and  repeated  pulses. 

If  repetitive  pulses  are  used,  the  pulse  repetition 
time  must  be  at  least  41.7  ysec  to  achieve  the  spacing  in 
delay.  This  would  correspond  to  a  pulse  repetition  frequency 
of  at  most  24  kHz  which  would  exceed  the  Doppler  spacing 
requirement  of  11.1  kHz.  If  we  selected  a  PRT  =  50  ysec,  a 


pulse  width  of  5  nsec,  and  a  train  of  10  pulses,  we  would 
achieve  the  desired  time-bandwidth  product. 

If  linear  FM  is  used,  the  bandwidth  to  be  swept  could 
be  25  MHz  and  the  pulse  width  20  psec.  This  would  provide 
the  required  time-bandwidth  product.  The  resolution  proper¬ 
ties  can  be  met  except  along  the  delay-Doppler  coupled  ridge. 
We  can  investigate  the  V-FM  waveform  which  will  reduce  the 
coupling  with  a  small  increase  in  complexity.  Therefore,  we 
will  now  run  the  programs  to  produce  the  three-dimensional, 
ambiguity  function  magnitude  diagrams  and  contour  plots  for 
the  repetitive  pulse,  linear  FM,  and  V-FM  waveforms. 

First  we  consider  the  repetitive  pulse  waveform.  As 
shown  in  Figure  2,  a  plot  of  the  magnitude  of  the  central 
portion  of  the  ambiguity  function,  the  Doppler  spacing  be¬ 
tween  the  3-dB  areas  is  approximately  17.9  KHz.  Observing 
Figures  3  and  4,  diagrams  of  the  peak  portion  of  the  ambi¬ 
guity  function  magnitude  and  3-dB  contour  of  magnitude 
squared,  respectively,  the  Doppler  resolution  is  approxi¬ 
mately  995  Hz,  and  the  delay  resolution  is  approximately  2.0 
psec.  Figure  5  indicates  that  the  delay  spacing  between  the 
3-dB  areas  is  approximately  45.0  usee.  As  shown  in  Figure  6, 
the  small  subsidiary  peaks  between  the  main  Doppler  areas 
are  at  most  9dB  below  the  main  regions.  Hence,  the  effects 
of  ambiguities  from  these-subsidiary  levels  is  negligible. 
Thus,  the  repetitive  pulse  waveform  exceeds  all  resolution 
and  ambiguity  requirements. 
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Next  we  consider  the  linear  FM  waveform.  As  shown  in 
Figures  7  and  8  (a  plot  of  the  ambiguity  function  magnitude 
and  a  plot  of  the  -3dB  and  -6dB  contours  of  the  magnitude 
squared,  respectively)  the  delay  resolution  is  approximately 
0.234  usee  while  the  Doppler  resolution  is  approximately  31 
kHz.  The  delay  resolution  is  well  within  limits;  however, 
the  Doppler  resolution  is  poorer  than  required.  Ambiguities 
exist  only  along  the  delay-Doppler  ridge;  however,  it  will  be 
virtually  impossible  to  resolve  targets  which  correspond  to 
points  along  the  ridge.  These  may  be  acceptable  to  the  user, 
especially  if  the  signal-to-clutter  ratio  is  high;  therefore 
we  will  continue  to  consider  this  waveform. 

As  shown  in  Figures  9  and  10  (a  plot  of  the  ambiguity 
function  magnitude  and  a  plot  of  the  -6dB  contour  of  the  mag¬ 
nitude  squared,  respectively) ,  the  V-FM  waveform  provides 
better  delay  and  Doppler  resolution  than  the  linear  FM  wave¬ 
form.  The  delay  and  Doppler  resolutions  for  V-FM  at  the 
-6dB  contour  correspond  to  those  for  linear  FM  at  the  -3d3 
contour.  With  V-FM,  however,  two  ridges  of  delay-Doppler 
ambiguities  exist  compared  to  the  one  ridge  for  linear  FM. 

In  any  case,  we  will  continue  to  consider  the  V-FM  waveform 
as  a  candidate  waveform. 

We  next  must  develop  a  clutter  model  for  the  operational 
environment.  By  modeling  the  environment,  we  can  then  deter¬ 
mine  the  signal-to-clutter  ratio  for  the  candidate  waveforms. 
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Since  we  assume  a  coastal  environment,  the  radar  will  receive 
sea  clutter,  ground  clutter,  and  occassionally  weather  clut¬ 
ter.  We  will  want  to  have  a  reasonably  accurate  model  for 
ail  three  types  of  clutter. 

The  first  clutter  model  we  will  develop  will  be  the  one 
for  the  sea  clutter.  From  the  discussion  of  clutter  pre¬ 
viously  provided  in  this  study,  we  know  that  the  spectrum  of 
sea  clutter  in  velocity  is  approximately  a  linear  function 
of  wind  speed.  The  user  must  determine  the  average  wind  speed 
for  the  general  location  of  the  radar.  For  our  purposes,  we 
will  assume  an  average  wind  speed  of  4  m/sec  (corresponding 
to  80  Hz) .  From  Figure  1  we  find  the  width  of  the  Doppler 
spectrum  to  be: 


2(2.3m/sec)  (3xl09sec~* 1 
8 

(3x10  m/sec) 


=  46  Hz 


We  also  assume  the  sea  clutter  to  be  distributed  as  Gaussian 
with  zero  mean  and  standard  deviation  of  15  Hz.  We  also 
assume  the  sea  clutter  to  be  distributed  uniformly  in  range 
from  the  radar  origin  oat  to  the  radar  horizon  (assume  it  to 
be  6  km).  On  the  delay  axis,  this  would  correspond  to  clutter 
from  0  to  10  ysec.  Hence  the  joint  probability  density  func¬ 
tion  would  be 


P(t'V)  ~  20usec 


7==-  — ?  exp  -  ( - j ) 

'2it  (15;  2(15; 


( 
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We  now  want  to  determine  the  average  cross  section  of 
the  sea  clutter  and  estimate  the  number  of  scatterers  from 
sea  clutter  in  the  radar  beam  at  any  time.  Since  the  wave¬ 
forms  under  consideration  have  pulse  widths  of  at  least  5 
ysec,  the  sea  clutter  amplitude  can  be  modeled  by  a  Gaussian 
or  Rayleigh  envelope  model.  The  mean  cross-sectional  area 
can  then  be  determined  by  using  Eq  72.  From  Figure  11  and 
assuming  a  grazing  angle  of  10°,  aQ  can  b;;  approximated  to 
-42dB  for  horizontal  polarization  and  -35dB  for  vertical 
polarization.  Since  we  will  use  horizontal  polarization,  we 
assume  oq  =  -42dB.  We  can  approximate  Ac  by  the  equation 

Ac  =  R6g (Cx/2)  sec  <p  (84) 

where 

R  =  range  of  clutter 

C  =  velocity  of  propagation 

0q  =  azimuth  beamwidth 

<j>  =  grazing  angle 

t  =  pulse  width  (Ref  11:472). 

Assuming  an  azimuth  beamwidth  of  3°  and  an  average  range  of 

5  2 

3  km,  we  approximate  Ac  as  1.2  x  10  m  for  the  5  ysec  pulse 
and  4.8  x  105m2  for  the  20  ysec  pulse.  Hence  the  average 

2 

cross  section  of  the  sea  clutter  will  be  approximately  7.6m 

2 

for  the  5  ysec  pulse  and  30.3m  for  the  20  ysec  pulse.  We 
estimate  that  we  would  receive  10  returns  from  sea  clutter 
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GRAZING  ANGLE  (fleg-ees) 

Fig  11.  Composite  of  o°  data  for  average  conditions 
with  wind  speeds  ranging  from  10  to  20  knots 


Gronnq  angle 


Fig  12.  L  Band,  Horizontal  Polarization  Values 
of  a  for  Various  Environments 
(Ref  11:475,493) 


at  the  same  time  we  would  want  to  detect  the  target. 

We  next  want  to  consider  the  ground  clutter.  Assuming 
the  area  will  be  flat  and  uncultivated,  we  can  use  the  graphs 
Figure  12  to  estimate  a°  =  -42dB  (Ref  11:495).  Hence  our 
average  clutter  returns  over  land  will  be  similar  to  those 
for  sea  clutter. 

To  determine  the  distribution  of  ground  clutter,  we 
recall  that  the  Doppler  spectrum  will  be  centered  at  zero. 

For  a  wind  speed  of  10  knots,  we  estimate  the  standard  devia¬ 
tion  of  the  Doppler  to  be  0.125  m/sec  (which  corresponds  to 
2.5  Hz).  Again  we  assume  the  clutter  to  be  uniformly  dis¬ 
tributed  in  range  from  the  radar  origin  out  to  the  radar 
horizon.  Hence  the  probability  density  function  for  the 
ground  clutter  in  terms  of  delay  and  Doppler  is: 


P  (t  ,  V)  = 


20ysec  /W.5) 


exp  (- 


2(2.5) 


r)> 


Finally,  we  want  to  consider  the  effect  of  weather 
clutter  on  the  detection  performance.  To  model  the  weather 
clutter,  the  user  will  need  weather  statistics  from  the 
general  operational  area.  If  we  assume  rain  showers  with  a 
rainfall  rate  of  4  mm/hr  occur  5  percent  of  the  time,  the 
clutter  amplitude  returns  can  be  modeled  by  Gaussian  (or 
Rayleigh  envelope)  probability  density  function.  The  average 
cross  sectional  area  can  be  estimated  at  approximately 
10  cm  /m  (cross  section  per  unit  volume)  (Ref  11:501). 
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Since  the  range  of  weather  clutter  can  be  from  the  radar 
origin  out  to  the  maximum  range  of  50  km,  very  large  cumula¬ 
tive  clutter  returns  can  be  expected.  We  assume  an  average 
2 

return  of’  0.04m  with  1000  returns  in  a  beamwidth. 

Assuming  a  mean  wind  speed  of  6  m/sec  (approximately 
18  mph) ,  the  mean  of  the  weather  clutter  velocity  will  be 
similar.  From  our  previous  discussion  of  weather  clutter, 
we  know  oturb  -  1.0  m/sec  and,  assuming  an  average  range  of 
20  km  and  a  velocity  gradient  component  of  4.0  m/sec/km, 
°shear  ~  m/sec.  Therefore,  the  total  spectrum  width,  c  , 

will  be  about  1.0  m/sec.  Again  we  assume  a  uniform  distri¬ 
bution  in  range  from  the  radar  origin  out  to  40  km  (corres¬ 
ponding  to  67  ysec) .  Hence  the  clutter  probability  density 
function  can  be  modeled  as: 


P  ( T  /  V) 


1 

134ysec 


exp 


(- (V-30) 2/ ( 5) 2) } 


By  using  the  clutter  plot  program,  we  now  can  plot  the 
clutter  distributions  as  functions  of  delay  and  Doppler.  A 
plot  of  the  sea  clutter  distribution  is  shown  in  Figure  13. 
Plots  of  the  ground  clutter  and  weather  clutter  are  shown  in 
Figures  14  and  15,  respectively.  Also  included  are  plots 
(near  the  clutter  regions)  of  the  magnitude  of  the  ambiguity 
functions  the  repeated  pulse,  LFM,  and  VFM  waveforms  in  Figures 
16,  17,  and  18,  respectively.  A  cursory  examination  of  the 
ambiguity  function  magnitudes  indicates  that  the  LFM  waveform 
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should  have  the  best  signal-to-clutter  ratio  since  its  delay 
width  is  much  narrower  than  the  repetitive  pulse,  and  since 
the  VFM  waveform  has  the  pedestal-like  magnitude  away  from 
the  ridge. 

Now  we  use  the  signal-to-clutter  program  to  determine 
quantitatively  how  each  waveform  will  perform  in  the  environ¬ 
ment.  We  assume  the  smallest  desired  target  has  a  radar  cross 
2 

section  of  4m  .  When  sea  clutter  was  considered,  the  signal- 
to-clutter  ratio  was  calculated  to  be  0.5  for  the  repetitive 
pulse  waveform,  2.7  for  the  LFM  waveform,  and  1.9  for  the 
VFM  waveform.  The  signal-to-clutter  ratio  over  the  ground 
clutter  was  calculated  to  be  0.5  for  the  repeated  pulse  wave¬ 
form,  2.7  for  the  LFM  waveform,  and  1.9  for  the  VFM  waveform. 
Finally,  the  signal-to-clutter  ratio  over  the  weather  clutter 
was  calculated  to  be  2.4  for  the  repeated  pulse  waveform,  139 
for  the  LFM  waveform,  and  97.9  for  the  VFM  waveform. 

Thus,  we  find  .at  none  of  the  candidate  waveforms  pro¬ 
vide  the  necessary  signal-to-clutter  ratio  to  detect  a  4m^ 
target  over  the  sea  and  ground  clutter  areas.  The  LFM  and 
VFM  waveforms  would  allow  detection  of  a  4mz  target  in  the 
weather  clutter,  while  the  repeated  pulse  waveform  would  not. 
The  user  must  now  consider  making  tradeoffs  between  the  wave¬ 
forms  and  requirements. 

If  the  user  is  willing  to  sacrifice  coverage  over  the 
clutter  areas,  then  the  repeated  pulse  waveform  should  be 
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chosen  since  it  exceeds  both  the  delay  and  Doppler  resolu¬ 
tion  requirements.  If  coverage  over  the  clutter  areas  is 
more  critical  than  resolution  requirements,  then  the  user 
will  choose  the  LFM  waveform.  A  much  better  signal-to-clutter 
ratio  can  be  obtained  by  using  an  LFM  waveform  with  a  5  usee 
pulse  width  and  a  frequency  modulation  band  of  100  MHz  (see 
Figure  19) .  With  this  waveform,  the  signal-to-clutter  ratio 
would  increase  to  43.8,  43.8,  and  55.7  over  the  sea,  ground, 
and  weather  clutter,  respectively.  The  delay-Doppler  ridge 
would  be  narrowed  in  delay,  but  would  be  extended  in  Doppler. 
Thus,  tha  user  must  determine  what  requirements  are  the  most 
important  and  then  make  tradeoffs  in  the  requirements  to 
determine  a  waveform  for  the  radar  system. 


( 
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Delay  (usee) 


Delay  (usee) 


( 


Fig  18.  Magnitude  of  Central 
of  VFM  Ambiguity  Fui 
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Magnitude  of  Central  Portion  of 
Revised  LFM  Ambiguity  Function 


Discussion 


As  stated  at  the  beginning  of  this  paper,  this  study 
was  to  develop  a  fast,  simple,  and  accurate  method  of  select¬ 
ing  a  radar  waveform  when  provided  the  target  detection 
requirements  and  the  expected  operating  environment.  The 
majority  of  work  for  this  study  was  involved  in  developing 
software  to  calculate  the  magnitude  of  radar  ambiguity  func¬ 
tions,  to  depict  these  magnitudes  in  three-dimensional  dia¬ 
grams,  and  to  calculate  the  signal-to-clutter  ratio  of  a 
desired  target  in  the  modeled  environment.  This  work  was 
limited  to  considering  only  the  most  common  waveforms  used 
in  radar  systems;  however,  these  waveforms  provided  the  basis 
for  more  complex  waveforms. 

Conclusions 

From  the  work  done  and  presented  in  this  study,  the 
method  proposed  for  waveform  selection  provides  the  fast, 
simple,  and  accurate  technique  desired.  By  calculating  and 
analyzing  the  magnitude  of  ambiguity  functions  and  the  signal- 
to-clutter  ratio,  the  waveform  selection  process  becomes 
simpler.  The  discussion  of  modeling  the  clutter  environment 
describes  how  to  reasonably  model  the  projected  operating 
environment  for  the  waveform  selection  process.  Also,  the 
discussion  of  basic  resolution  properties  of  waveforms  allows 
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the  user  to  select  only  those  waveforms  with  a  reasonable 
probability  of  success  for  detailed  analysis. 

Once  a  list  of  candidate  waveforms  is  selected,  the 
three-dimensional  diagram  of  the  ambiguity  function  magnitude 
(as  a  function  of  delay  and  Doppler)  for  each  waveform  allows 
the  user  to  insure  the  resolution  and  ambiguity  requirements 
are  fulfilled.  With  proper  sampling,  these  magnitudes  are 
within  one  percent  of  the  actual  results.  For  those  wave¬ 
forms  fulfilling  these  requirements,  the  signal-to-clutter 
ratio  for  desired  targets  is  calculated  for  each  waveform. 

The  discussion  on  tradeoffs  in  this  paper  describes  how  the 
user  should  continue  the  selection  process  when  several  wave¬ 
forms  fulfill  the  requirements  or  when  none  of  the  waveforms 
fulfill  the  requirements. 

To  use  this  method  of  waveform  selection,  the  user 
must  have  access  to  a  computing  system  with  a  high-level 
language  capability,  at  least  32K  bytes  of  memory,  a  fast 
Fourier  transform  program,  and  a  three-dimensional  graphics 
program.  All  software  developed  in  this  study  will  provide 
accurate  results  as  long  as  the  proper  parameters  are  input 
to  the  program.  The  user  must  have  a  basic  understanding  of 
radar  theory  and  digital  signal  processing  techniques  to  suc¬ 
cessfully  use  the  software.  A  brief  discussion  of  both  topics 
is  included  in  this  paper. 


Recommendations 


As  stated  previously,  this  study  was  limited  to  wave¬ 
forms  generally  used  in  contemporary  radars.  Ambiguity  func¬ 
tions  using  azimuth  angle  and  elevation  angle  were  also  not 
considered.  The  study  also  only  considered  the  case  where 
there  was  negligible  distortions  of  the  complex  envelope  of 
the  waveform  by  target  motion,  target  acceleration,  and 
higher  order  derivatives. 

Many  other  radar  waveforms  exist  other  than  those  con¬ 
sidered  in  this  study.  As  technology  improves,  the  cost  to 
implement  these  waveforms  in  radar  systems  will  decrease  and 
make  their  use  more  appealing.  Thus,  further  study  of  addi¬ 
tional  waveforms  is  suggested. 

The  ambiguity  functions  calculated  using  azimuth  angle 
and  elevation  angle  can  significantly  contribute  to  the  de¬ 
sign  and  selection  of  antennas  for  radar  systems.  All  re¬ 
quired  theory  to  calculate  these  ambiguity  functions  exists, 
and  therefore,  a  method  for  selecting  an  antenna  for  the 
radar  can  be  developed  similar  to  the  waveform  selection 
method  of  this  study. 

Although  there  is  a  wide  range  of  practical  applications 
in  which  target  motion  does  not  distort  the  waveform,  situa¬ 
tions  exist  in  which  a  radar  must  detect  high-speed  targets. 
Extension  of  the  theory  for  such  situations  exist,  but  the 
theory  becomes  much  more  complex.  Consideration  should  be 


^  given  to  developing  a  method  of  waveform  selection  when  high¬ 

speed  targets  must  be  detected. 

An  improved  three-dimensional  graphics  program  would 
improve  the  waveform  selection  process.  Because  the  graphics 
program  used  in  this  study  was  very  limited  in  the  number  of 
samples  it  could  display,  the  resolution  and  range  of  the 
three-dimensional  diagrams  was  only  adequate.  With  a  pr  ■'gram 
which  could  plot  256  x  256  arrays,  the  use  of  these  die  ams 
would  increase  and  fewer  diagrams  would  be  required  in 
analysis. 

Finally,  the  performance  of  receivers  which  do  not  use 
a  matched  filter  can  still  be  evaluated.  The  cross-ambiguity 
function  can  be  determined  and  used  in  the  waveform  selection 
process.  This  study  made  no  provision  to  calculate  the  cross¬ 
ambiguity  function;  however,  the  development  of  software  to 
consider  cross-ambiguity  functions  could  follow  directly  from 
that  for  the  ambiguity  function.  Hence,  the  cross-ambiguity 
function  software  should  be  developed. 


( 
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APPENDIX  A 


Operation  of  Software  on  the  Eclipse  Computer 

The  basic  principles  of  operating  the  Eclipse  computer 
are  explained  so  that  the  following  items  can  be  produced: 
the  magnitude  of  cuts  of  ambiguity  functions,  three-dimensional 
diagrams  of  the  magnitude  of  the  ambiguity  function,  and  the 
signal-to-clutter  ratio.  The  Eclipse  uses  the  Command  Line 
Interpreter  (CLI)  as  the  basic  operating  command  set.  When 
the  Eclipse  can  accept  a  command,  an  'R'  prompt  appears  on 
the  terminal.  The  user  should  enter  the  command,  DIR  DP4 : 

REED  < CR> ,  which  puts  the  user  in  the  directory  containing 
the  software  programs. 

The  user  should  enter  the  command,  MAIN  <CR>,  to  run 
the  program  which  simply  provides  the  magnitudes  of  the  am¬ 
biguity  function  for  a  particular  Doppler  cut.  To  calculate 
the  signal-to-clutter  ratio,  the  user  should  enter  the  com¬ 
mand,  IMAIN  <CR> .  When  a  three-dimensional  plot  of  the  magni¬ 
tude  of  the  ambiguity  function  is  desired,  the  user  should 
enter  the  command,  GMAIN  <CR>.  For  a  three-dimensional  plot 
of  the  clutter  probability  density  function,  the  user  should 
enter  the  command,  CLUTR  <CR>.  After  GMAIN  or  CLUTR  runs, 
the  user  will  then  have  to  use  the  Textronix  storage  scope 
to  observe  the  three-dimensional  plot.  Once  the  control  box 
for  the  foreground  is  set  to  the  Textronix  scope,  the  user 
should  enter  the  command,  PLTTRNS  DATA/I  Num/N  <CR>  for  the 
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ambiguity  function  magnitude  plot  or  PLTTRNS  CLUTTER/I  64 /N 
for  the  clutter  plot.  The  value  of  Num  will  be  specified  at 
the  end  of  the  GMAIN  program,  and  the  user  must  note  it  to 
enter  when  PLTTRNS  is  run. 

When  MAIN,  GMAIN,  IMAIN,  or  CLUTR  is  run,  the  user  will 
be  required  to  make  entries  during  the  program.  These  entries 
will  specify  the  data  to  enter,  and  all  data  should  be  entered 
in  decimal  notation.  A  carriage  return,  <CR>,  is  required 
after  all  entries.  If  the  user  makes  a  mistake  in  the  entry 
while  using  the  Background  terminal,  the  DEL  key  may  be  used 
to  erase  characters.  The  '/'  key  may  be  used  to  delete  the 
entire  entry  on  both  the  Foreground  and  Background  terminals. 
The  user  may  enter  ' ^A'  at  any  time  to  terminate  the  program 
being  run.  When  a  program  has  ended,  the  user  will  see  an 
1 R'  prompt  on  the  terminal. 

When  PLTTRNS  is  run,  the  user  must  use  the  Textronix 
scope.  This  program  will  also  require  the  user  to  make  en¬ 
tries.  If  the  entry  is  an  alphabetic  character,  no  <CR>  is 
required.  After  any  numerical  entry,  the  user  must  use  a 
<CR>  to  indicate  the  end  of  the  number.  If  a  copy  of  the 
plots  on  the  scope  is  desired,  the  user  simply  depresses  the 
COPY  key  on  the  Textronix  keyboard.  The  user  should  insure 
the  copier  has  been  on  for  at  least  ten  minutes  before  copies 
are  made.  The  diagrams  produced  by  PLTTRNS  do  not  have  the 
axes  identified  or  labeled.  The  user  can  obtain  the  delay 
and  Doppler  limiting  values  from  the  GMAIN  or  CLUTR  program. 
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If  copies  of  the  program  listings  are  desired,  the  user 
should  enter  the  command,  'PRINT  filename. FR,  f ilename .FR, . . . 
filename. FR' .  The  program  listings  indicate  all  the  subrou¬ 
tines  called,  and  these  are  the  filenames  to  enter  along  with 
the  program  name. 

After  the  user  has  completed  the  use  of  all  the  programs, 
the  computer  should  be  returned  to  the  DP0  directory.  This 
is  done  by  entering  the  command,  DIR  DP0  <CR>. 


( 
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APPENDIX  B 


System  Flowchart 


Flowchart  for  Calculating  Magnitude  of  Ambiguity  Function 


Run  Program  Main 


END 
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Flowchart  for  Three-Dimensional  Diagram  of 
Ambiguity  Function  Magnitude 


Run  Program  'Main' 


End  of  Gmain 


Run  graphics  program  ' PLTTRNS '  to  pro¬ 
duce  the  three-dimensional  plots  and 
contour  plots _ _ 
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Flowchart  for  Calculating  Signal -to-Cl utter  Ratio 


Run  Program  'Imain' 


Continued 


Continuation 


END 
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APPENDIX  C 

Plots  of  the  Magnitude  of  Ambiguity  Function 
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-3dB  and  -6dB  Contour  of  5-F.lomenl  Barker  Code 


Repented  Gaussian  Pulses 


Delay  (sec) 


Delay  (sec) 


n  o  r»  o  n 


*  PROGRAM  MAIN 

«  WRITTEN  BY.  CAPT  JOHN  REED 

•  IS  APRIL  1982 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


MAIN  IS  A  FORTRAN  V  PROGRAM  WHICH  CALCULATES  THE  MAGNITUDE  OF  A  CUT  OF 
THE  AMB1GUTITY  FCN  FOR  A  SPECIFIED  DOPPLER  FREQUENCY.  THE  OUTPUT  IS  THE 
MAGNITUDE  AT  VARIOUS  DELAY  TIMES  (TAU)  FOR  THE  SPECIFIED  DOPPLER 
FREQUENCY  (W). 

THE  FOLLOWING  ARE  USER  SELECTED  PARAMETERS  WHICH  THE  PROGRAM  (OR 
THE  SUBROUTINES  WHICH  IT  CALLS)  WILL  REQUEST: 

1.  PULSE  WIDTH. 

2.  WHETHER  OR  NOT  THE  PULSE  IS  REPEATED. 

A.  WHETHER  PULSE  REPETITION  IS  STAGGERED  OR  THE  SAME. 

B.  PULSE  REPETITION  TIME(S)  (PRT) 

C.  NUMBER  OF  REPETITIONS 

3  NUMBER  OF  SAMPLES  TO  BE  MADE  OF  THE  WAVEFORM. 

4.  WHETHER  OR  NOT  THE  PULSE(S)  IS  PHASE-SHIFT  CODED. 

A.  NUMBER  OF  ELELMENTS  IN  THE  PULSE  CODE. 

B.  THE  PHASE  OF  EACH  ELEMENT. 

3  WHETHER  OR  NOT  THE  PULSE(S)  IS  FREQUENCY  MODULATED. 

A.  BANDWIDTH  TO  BE  SWEPT  BY  THE  FM. 

B.  IF  LINEAR  FM  OR  VFM. 

6.  WHETHER  THE  PULSE(S)  IS  RECTANGULAR  OR  GAUSSIAN  SHAPED. 

7.  DOPPLER  FREQUENCY.  U  (RADIANS),  FOR  WHICH  THE  MAGNITUDE  DF 
THE  AMBIGUITY  FUNCTION  IS  TO  BE  CALCULATED. 

MAIN  MAY  MAKE  CALLS  TO  THE  FOLLOWING  SUBROUTINES: 

INPUT,  PHASE.  FM,  CPULS,  RECPL,  SETUP.  OUT. 

IT  ALSO  CALLS  THE  FAST  FOURIER  TRANSFORM  SUBROUTINE  (EDFT  LB >  FROM 
THE  COMPUTER  LIBRARY. 

MAIN  IS  COMPILED  AND  LOADED  ON  THE  ECLIPSE  USING  THE  FOLLOWING  COMMANDS. 
FORTRAN  MAIN 

RLDR/N  MAIN  INPUT  PHASE  FM  GPULS  RECPL  SETUP  OUT  EDFT  LB  «FLIB« 

THIS  PROG°AM  IS  STORED  IN  FILE  'MAIN. FR '  IN  THE  DIRECTORY  DP4  REED. 

THE  FOLLOWING  VARIABLES  ARE  USED  IN  'MAIN': 

PW  -  PULSE  WIDTH 

SAME  -  PRT  IF  PULSES  ARE  REPEATED  AT  SAME  INTERVAL 

NPRT  -  •  OF  REPETITIONS  OF  PULSE  TO  BE  MADE 

PRT  -  ARRAY  CONTAINING  THE  PRT  FOR  EACH  REPETITION 

TOTAL  -  TOTAL  TIME  DURATION  OF  RADAR  WAVEFORM 

NUM  -  t  OF  SAMPLES  MADE  OF  WAVEFORM 

SAMPLE-  •  OF  SAMPLES  OF  WAVEFORM  PLUS  7ER0-F ILL INO  SAMPLES  NEEDED 
TO  PERFORM  LINEAR  CONVOLUTION  ,»2*NUM> 

PERIOD-  ARRAY  CONTAINING  #  OF  SAMPLES  MADE  DURING  EACH  PRT 
N  -  *  OF  SAMPLES  MADE  OF  PULSE  AFTER  WEIGHTING 

M  -  #  OF  NON-ZERO-FILL  SAMPLES  MADE  Or  WAVEFORM  AFTER  WEIGHTING 

Z  -  ARRAY  CONTAINING  THE  AMPLITUDE  OF  EACH  SAMPLE 
GFACTOR-SCALINO  VALUE  OF  MAGNITUDE 

FREQ  -  ARRAY  CONTAINING  FREQUENCY  MODULATION  FOR  EACH  SAMPLE 
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C  SCALE 

C  VALUE 

C  XVALS 

c 
c 

C  YVAL3 

C 

C  MAG  - 

C 

COMPLEX  XVALS ( 1024),  YVALS( 1024) ,  Z ( 1024 ! ,  VALUE 
INTEGER  SAMPLE, PERIOD* 10) 

REAL  MAG  < 1024).  PRT< 10)  PREQ<1024) 

TYPE  '  WELCOME  TO  MAIN' 

TYPE  'THIS  PROGRAM  IS  DESIGNED  TO  COMPUTE  THE  MAGNITUDE  OF  THE' 

TYPE  'AMBIGUITY  FUNCTION  OF  A  RADAR  WAVEFORM  YOU  WILL  BE  ASKED' 

TYPE  'TO  ENTER  PARAMETERS  TO  DEFINE  THE  WAVEFORM  AT  VARIOUS' 

TYPE  'TIMES  DURING  THE  PROGRAM  WHEN  THE  PROGRAM  REQUESTS  AN' 

TYPE  'ENTRY,  SIMPLY  TYPE  YOUR  REPLY  AND  , IF  CORRECT,  HIT  CR  KEY  ' 

TYPE  'IF  YOU  MAKE  A  MISTAKE  ON  AN  ENTRY,  USE  THE  <DEL>  KEY  TO' 

TYPE  'ERASE  THE  CHARACTER  OR  THE  <\>  KEY  TO  DELETE  THE  ENTIRE' 

TYPE  ENTRY  YOU  MAY  STOP  THE  PROGRAM  AT  ANY  TIME  BY  ENTERING' 

TYPE  ' < " A>  THE  OUTPUT  WILL  BE  THE  MAGNITUDE  OF  THE  AMBIGUITY' 

TYPE  'FUNCTION  AT  VARIOUS  DELAY  TIMES  ( TAU  FOR  THE  DOPPLER  FREQUENCY' 
TYPE  'SPECIFIED  EARLIER  TO  STOP  THE  SCREEN  DISPLAY  OUTPUT,  ENTER' 
TYPE  'OSIX;  TO  RESUME  THE  DISPLAY  OUTPUT,  ENTER  <"Q>.  ' 

PAUSE  'ENTER  ANY  CHARACTER  TO  BEGIN' 

call  I NPUT  <  PER  I OD.  PRT,  PW,  NPRT,  TOTAL,  SAMPLE, NUM, M, N) 

ACCEPT  ENTER  0  IF  PULSE<S)  IS  PHASE-SHIFT  CODED .  1  IF  NOT  ',J 

IF  <U  EQ.  1  )  GO  TO  5 

CALL  PHA5E(FREQ.  PRT,  PERIOD,  Z, M. N,  NPRT,  GF ACTOR,  PW,  SAMPLE,  TOTAL) 

GO  TO  20 

5  TYPE  PULSE <S>  IS  NOT  PHASE-SHIFT  CODED' 

ACCEPT  ENTER  0  IF  PULSE(S)  IS  FREOENC Y  MODULATED,  1  IF  NOT  ',U 
CALL  FM ( FREQ,  PERIOD,  NPRT,  PW,  N,  J,  SAMPLE) 

10  ACCEPT  'ENTER  0  FOR  RECTANGULAR-SHAPED  PULSECS),  1  FOR  GAUSSIAN  ',J 

IF  ( J  EQ.  0)  GO  TO  15 

CALL  GPULSIPERIOO,  Z,  PW,  GFACTCR, NPRT,  N) 

GO  TO  20 

15  CALL  RECPUPERIOD,  Z,  GFACTOR,  NPRT,  N) 

20  SCALE=GFACTOR/<N*NPRT) 

ACCEPT  'ENTER  VALUE  OF  W  TO  BE  CONSIDERED  ',W 
25  VALUE=CMPLX <0  0. W> 

WRITE  <10,300)  W 

CALL  SETUP*  Z,  FREQ,  XVALS,  YVALS,  VALUE,  TOTAL,  SAMPLE,  M) 

CALL  DFT5< XVALS. SAMPLE, 0) 

CALL  DFT5(YVALS, SAMPLE. 0) 

DO  30  1*1. SAMPLE 

XVALS < I )=CONJG<  XVALS < I )  ) 

YVALS  < I >  *CONJG  <  YVALS  < I ) ) 

XVALS* I )  =  XVALS  *  I )*YVAL5( I  >/SAMPLE 
30  CONTINUE 

CALL  DCT5<  XVALS. SAMPLE,  1 ) 


-  SCALING  FACTOR  OF  MAGNITUDE 

-  COMPLEX  EXPONENTIAL  OF  DOPPLER  FREQUENCY  <EXP(J*W>> 

-  ARRAY  CONTAINING  THE  TIME-REVERSED.  CONJUGATE  SAMPLES  OF 
THE  COMPLEX  ENVELOPE  OF  THE  WAVEFORM;  ALSO  CONTAINS  THE 
AMBIGUITY  FUNCTION  VALUES  LATER  IN  THE  PROGRAM 

-  array  CONTAINING  ThF.  DOPPLER  AND  RANGE-DELAYED  SAMPLES 
OF  THE  COMPLEX  ENVELOPE  OF  WAVEFORM 

ARRAY  CONTAINING  THE  MAGNITUDE  OF  THE  AMBIGUITY  FCN  CUT 
FOR  EACH  SAMPLE 
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DO  35  1=1. SAMPLE 
MAG  < I >=CABS( XVALSi I ) ) »SCALE 
35  CONTINUE 

CALL  OUT (MAG. NUM.  PW,  TOTAL,  m 

ACCEPT  'ENTER  999  TO  STOP  OR  NEXT  U  TO  CONTINUE  ',W 
IF  (W  £Q.  999)  GO  TO  500 
GO  TO  25 

300  FORMAT <  'W='. F 12  2) 

50 0  '  END 
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c*  subroutine  input 

C*  WRITTEN  BV  CAPT  JOHN  REED 

C*  15  APRIL  1982 

. . . . . 
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INPUT  IS  A  FORTRAN  V  SUBROUTINE  WHICH  REQUESTS  INFORMATION  FROM 
THE  USER  TO  DEFINE  THE  BASIC  RADAR  WAVEFORM 

THE  INPUT  SUBROUTINE  IS  CALLED  BV  THE  FOLLOWING  PROGRAMS 
MAIN,  GMAIN,  IHAIN 

THIS  SUBROUTINE  RECEIVES  NO  VARIABLES  FROM  THE  CALLING  PROGRAM 
IT  RETURNS  THE  FOLLOWING  VARIABLES  TO  THE  CALLING  PROGRAM 
PERIOD  PRT,  PW.  NPRT,  TOTAL.  SAMPLE.  NUM,  M,  N 

INPUT  IS  COMPILED  ON  THE  ECLIPSE  USING  THE  COMMAND 
FORTRAN  INPUT 

THIS  SUBROUTINE  IS  STORED  IN  FILE  'INPUT.  FR '  IN  THE  DIRECTORY  DP 4  REED 

THE  following  VARIABLES  are  used  in  this  subroutine. 

PW  -  PULSE  WIDTH 

SAME  -  PRT  IF  PULSES  ARE  REPEATED  AT  THE  SAME  INTERVAL 
NPRT  -  #  OF  REPETITIONS  OF  PULSE  TO  BE  MADE 

PRT  -  ARRAY  CONTAINING  THE  PRT  FOR  EACH  REPETITION 
TOTAL  -  TOTAL  DURATION  OF  WAVEFORM  (IN  TIME) 

NUM  -  (•  OF  SAMPLES  TO  BE  MADE  OF  WAVEFORM 

SAMPLE-  #  OF  SAMPLES  OF  WAVEFORM  PLUS  THE  ZERO-FILL  SAMPLES 

NEEDED  TO  PERFORM  THE  LINEAR  CONVOLUTION  (=2»NUM> 

N  -  #  OF  SAMPLES  MADE  OF  EACH  PULSE  AFTER  PROPER  WEIGHTING 

M  -  #  OF  NON-ZERO-FILL  SAMPLES  MADE  OF  WAVEFORM  AFTER  WEIGHTING 

PERIOD-  ARRAY  CONTAINING  «  OF  SAMPLES  MADE  OF  EACH  PRT 


SUBROUTINE  INPUT  <  PER  I OD*  PRT,  PW,  NPRT,  TOTAL,  SAMPLE,  NUM,  M,  N) 


INTEGER  SAMPLE, PERIOD! 10) 
REAL  PRT < 10) 


ACCEPT  ENTER  PULSE  WIDTH  '.PW 

ACCEPT  'ENTER  1  IF  WAVEFORM  HAS  REPEATED  PULSES,  0  IF  NOT  ',  I 

IF  ( I  £0  0)  GO  TO  25 

TYPE  WAVEFORM  HAS  REPEATED  PULSES' 

ACCEPT  ENTER  1  IF  STAGGERED  PRT,  0  IF  PRT  SAME  U 
IF  ( J  EQ  1 )  CO  TO  10 

ACCEPT  'ENTER  PULSE  REPITITION  TIME  '.SAME 
ACCEPT  'ENTER  *  OF  REPIT1TICNS  TO  BE  MADE  ' ,  NPRT 
DO  5  1*1.  NPRT 
PRT ( I ) *SAME 
5  CONTINUE 

total*same»nprt 

WRITE!  10.  300)  SAME,  NPRT 
GO  TO  30 

10  TYPE  "STAGGERED  PRT'S  TO  BE  USED" 

ACCEPT  ENTER  »  OF  DIFFERENT  PRT  TO  BE  USED  '.NPRT 


ACCEPT  'ENTER  FIRST  PRT  ',PRT.'1> 

DO  15  1=2, NPRT 

ACCEPT  'ENTER  NEXT  PRT  ',PRT(I> 

15  CONTINUE 

T0TAL=0  0 
DO  20  1=1. NPRT 
TOTAL-TOTAL+PRT! I ) 

WRITE! 10. 305)  l.PRT(l) 

20  CONTINUE 

CO  TO  30 

25  TYPE  'WAVEFORM  HAS  ONLY  ONE  PULSE' 

TOTAL=PW 
NPRT=1 
PRT ( 1 ) =PW 

ACCEPT  'ENTER  #  OF  SAMPLES  TO  BE  MADE  OF  WAVEFORM  ' , NUM 

SAMPLE=NUM*2 

PER  I OD  < 1 )=NUH 

M=NUM 

N=NUM 

TYPE  'N=',N,  '  M='.M,  '  PERIOD!  1  )=',  PERIOD(  1  ) 

CO  TO  33 

30  TYPE  'ENTER  #  OF  SAMPLES  TO  BE  MADE  OF  WAVEFORM;  YOU  MAY' 

ACCEPT  '  CHOOSE  16, 32. 64, 12S. 256, OR  512  ',NUM 
SAMPLE=NUM*2 

IF  <PW  GT  O  0001)  CO  TO  32 
CPW=PW* 10000  0 
CTOTAL=TOTAL* 10000  0 
N=NUM/C70TAL*CPW 
M=N*CTOTAL/CPW 
DO  31  1=1. NPRT 
CPRT< I)=PRT! 11*10000.  0 
PERIOD! I )=N*CPRT!I)/CPW 
TYPE  ' N= ' ,  N.  '  M= ' ,  M 

TYPE  'I=',I.  '  PERIOD! I)  = '. PERIOD! I ) 

31  CONTINUE 
GO  TO  33 

32  N—NUM/ TOTAL *PW 
M=N*TOTAL/PW 
TYPE  'M=',M 
TYPE  'N=',N 

DO  35  1=1, NPRT 

PERIOD! I )=N*PRT! I )/PW 

TYPE  'I=',I. '  PERIOD! I )='. PERIOD! I) 

35  CONTINUE 

33  RETURN 

300  FORMAT! "PRT 'S  ARE  THE  SAME  PRT=" . F 1 2  10,  5X,  '#  OF  REP  I T I T 10NS= ' ,  12 ) 

305  FORMAT!  'FOR  I = ' .  1 3. 5X,  ' PRT- ' , F 1 2  10) 

END 
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*  SUBROUTINE  PHASE  * 

*  WRITTEN  BY;  CAPT  JOHN  REED  * 

*  15  APRIL  1982  • 
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PHASE  IS  A  FORTRAN  V  SUBROUTINE  WHICH  CALCULATES  THE  AMPLITUDE  OF 
THE  SAMPLES  OF  A  PULSE  WHICH  IS  PHASE-SHIFT  CODED.  IT  ALSO  SETS  THE 
FREQUENCY  OF  EACH  SAMPLE  TO  ZERO  SINCE  FM  IS  NOT  AN  ALLOWED  OPTION. 

THE  PHASE  SUBROUTINE  MAY  BE  CALLED  BY  THE  FOLLOWING  PROGRAMS; 

MAIN,  CHAIN,  1MA1N. 

THIS  SUBROUTINE  RECEIVES  THE  FOLLOWING  VARIABLES  FROM  THE  CALLING  PROGRAM: 

PRT,  PERIOD.  M,  N,  NPRT.  PW.  SAMPLE,  TOTAL. 

IT  RETURNS  THE  FOLLOWING  VARIABLES  TO  THE  CALLING  PROGRAM: 

FREQ,  PRT,  PERIOD,  Z,  M,  N,  GFACTOR.  TOTAL. 

'PHASE'  IS  COMPILED  ON  THE  ECLIPSE  USING  THE  COMMAND: 

FORTRAN  PHASE 

'PHASE'  IS  STORED  IN  FILE  'PHASE.  FR '  IN  THE  DIRECTORY  DP4:  REED. 

THE  FOLLOWING  VARIABLES  ARE  USED  IN  SUBROUTINE  PHASE: 

NELEM  -  #  OF  ELEMENTS  IN  THE  PULSE  CODE 

IT I ME  -  #  OF  SAMPLES  TO  BE  MADE  OF  EACH  ELEMENT 

SUM  -  TOTAL  #  OF  SAMPLES  MADE  OF  EACH  PULSE  AFTER  WEIGHTING 

PERIOD-  ARRAY  CONTAINING  »  OF  SAMPLES  MADE  OF  EACH  PRT  AFTER  WEIGHTING 

RATIO  -  ARRAY  CONTAINING  THE  PHASE  OF  EACH  ELEMENT  I /PI) 

MEAN  -  1/2  LENGTH  OF  AN  ELEMENT 

STDEV  -  1/3  OF  MEAN)  3  STANDARD  DEVIATIONS  FROM  MEAN 
VAR  -  VARIANCE  ( STDEV*STDEV> 

GAUSS  -  SCALING  VALUE  OF  GAUSSIAN  ELEMENTS  C 1 /SORT ( 2*P I*VAR > > 

GFACTOR -SCALING  VALUE  OF  MAGNITUDE  OF  AMBIGUITY  FCN 
TIME  -  VALUE  OF  DELAY  (TAU) 

POWER  -  (TIME-MEAN) 

WT  -  VALUE  OF  GAUSSIAN  PULSE  (DISTRIBUTION)  FOR  A  PARTICULAR  TIME 
AMT2  -  IMAGINARY  VALUE  -  RATIO*PI 

Z  -  ARRAY  CONTAINING  AMPLITUDE  OF  EACH  SAMPLE  OF  WAVEFORM 

FREQ  -  ARRAY  CONTAINING  THE  FM  OF  EACH  SAMPLE 

M  -  TOTAL  (*  OF  NON-ZERO-FILL  SAMPLES  OF  WAVEFORM  AFTER  WEIGHTING 
N  -  »  OF  SAMPLES  OF  EACH  ELEMENT  AFTER  WEIGHTING 


SUBROUTINE  PHASE < FREQ,  PRT,  PERIOD, Z. M, N, NPRT, GFACTOR, PW, SAMPLE, TOTAL) 


COMPLEX  Z< 1024) . AMT2 

INTEGER  PERIODC 10),  COUNT,  SUM,  TALLY.  SAMPLE 
REAL  PRT ( 10). RATI0(20>. FREO< 1024), MEAN 

TYPE  'PULSES  ARE  PHASE-SHIFT  CODED.  .  NO  OTHER  OPTIONS  AVAILABLE’ 

ACCEPT  'ENTER  •  ELEMENTS  IN  THE  PULSE  CODE  '.NELEM 

ITIHE-N/NELEM 

SUM«ITIME»NELEM 

TYPE  'SUM-'.  SUM 

IF  (NPRT.  NE  1)  GO  TO  1 
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PERIOD* 1 >=SUM 

TVPE  PERIOD! 1 >  =  ',  PERIOD!!  1 
50  TO  6 

1  IF  (PW  CT  0  0001)  GO  TO  3 
DO  2  1*1. NPRT 

CPRT  < I ) =PRT ( I )*10000  0 

CPW=PU*10000  0 

PERIOD! I )=SUM»CPRT! I > /CPW 

TVPE  I.  '  PERIOD! I  1= PERIOD! I ) 

2  CONTINUE 
CO  TO  6 

3  DO  5  1=1, NPRT 

PERIOD! I >=SUM*PRT< I > /PW 

TVPE  'I=',I,'  PER  I OD ! I ) = ' . PER  I OD ( I ) 

5  CONTINUE 

6  DO  10  1=1 , NEl EM 

ACCEPT  ENTER  PHASE  AS  A*PI/Bi  ENTER  A/B  RATIO  ' 
WRITE! 10.  300)  I- RATIO! I) 

10  CONTINUE 

ACCEPT  'ENTER  0  FOR  RECTANGULAR-SHAPED  ELEMENTS, 

IF  <  J  EG. 0)  GO  TO  35 

TYPE  'ELEMENTS  ARE  GAUSSIAN  SHAPED' 

MEAN=PW/!2  0*NELEM> 

STDEV=MEAN/3.  0 

UAR=STDEV*STDEV 

GAUSS* 1  O/SGRT < 2  0*VAR*3.  1416) 

GF AC  TOR =0.  5923581  * ( PW/NELEM ) * ! PW/NELEM ) 

WRITE! 10.  305)  MEAN,  STDEV.  VAR,  GF ACTOR 
TALLY=0 

DO  30  IA=1 , NPRT 

C0UNT=TALLY*1 

DO  20  1*1,  NELEM 

DO  15  J=COUNT , COUNT  + 1  TIME- 1 

TIME-IPW/ITIME)*! J-COUNT+ 1 ) / NELEM 

POUER=TIME-MEAN 

WT=GAUSS*EXP!-P0WER*POWER/!2  0*VAR ) ) 

AMT2=CMPLX!0  0, RATIO! I ) *3. 1416) 
Z<U>=WT*CEXP<AMT2) 

15  CONTINUE 

C  OUNT =C OUNT  + I T I ME 
20  CONTINUE 

DO  25  J=COUNT . COUNT *PER I OD! I A ) -SUM- 1 
Z<  J>»(0  0.  0.  0) 

25  CONTINUE 

TALLV=TALLV*PERIOD< IA) 

30  CONTINUE 

GO  TO  60 

35  TYPE  'ELEMENTS  ARE  RECTANGULAR  SHAPED' 

GFACTOR*!  0 
TALLY=0 

DO  55  I A=1 . NPRT 

COUNT.TALLY+1 

DO  45  1*1. NELEM 

DO  40  U=COUNT- C0UNT*ITIME-1 


( 


.  RATIO! I ) 

1  FDR  GAUSSIAN 
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AMT2=CMPLX<0.  0.  RATIO<  I)*3.  1416) 

Z(U>=1  0*CEXP  <  AMT2) 

40  CONTINUE 

C  OUNT =C  OUNT* I T I ME 
45  CONTINUE 

DO  50  J=COUNT ,  COUNT ♦PER  1 0D< I A ) -SUM-1 
Z(  J)  =  <0.  0,  0  0) 

50  CONTINUE 

TALLY=TALLY+PERIOD< IA) 

55  CONTINUE 

60  DO  65  1=1, SAMPLE 

FREQ<  1  )  =0.  0 
65  CONTINUE 

M=SUM*TOTAL/PW 

N-SUM 

RETURN 

300  FORMAT <  'FOR  I  =  ' ,  12, 5X,  'PHASE* ' ,  F5.  2,  '*P I  ' ) 

305  FORMAT (  'M=',  F12  10,  3X,  'SD=',  FI 2.  10,  3X,  'VAR=',  F12  10,  3X,  'GFAC=  ' ,  FI 2.  10) 

END 
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c ****************************** ********** 

C*  SUBROUTINE  FM  * 

C*  WRITTEN  BY:  CAPT  JOHN  REED  • 

C*  IS  APRIL  1982  * 

C**************************************** 

' FM '  IS  A  FORTRAN  V  SUBROUTINE  WHICH  DETERMINES  THE  FREQUENCY 
MODULATION  OF  EACH  SAMPLE  OF  THE  RADAR  WAVEFORM  (EXCEPT  FOR 
PHASE-SHIFT  CODED  SIGNALS). 

THE  FM  SUBROUTINE  MAY  BE  CALLED  BY  THE  FOLLOWING  PROGRAMS: 

MAIN.  GMAIN.  I  MAIN. 

THIS  SUBROUTINE  RECEIVES  THE  FOLLOWING  VARIABLES  FROM  THE  CALLING  PROGRAM 
PERIOD.  NPRT,  PW,  N.  J,  SAMPLE. 

IT  RETURNS  THE  FOLLOWING  VARIABLES  TO  THE  CALLING  PROGRAM: 

FREG. 

'FM'  IS  COMPILED  ON  THE  ECLIPSE  USING  THE  COMMAND: 

FORTRAN  FM 

THIS  SUBROUTINE  IS  STORED  IN  FILE  'FM.  FR '  IN  THE  DIRECTORY  DP4  REED 

THE  FOLLOWING  VARIABLES  ARE  USED  IN  'FM': 

J  -  SPECIFIES  IF  THE  WAVEFORM  IS  FREQUENCY  MODULATED  OR  NOT 
BAND  -  BANDWIDTH  OF  THE  FM 
TIME  -  VALUE  OF  TAU  (DELAY) 

FREQ  -  ARRAY  CONTAINING  THE  FM  OF  EACH  SAMPLE 
PERIOD-  ARRAY  CONTAINING  THE  #  OF  SAMPLES  OF  EACH  PRT 
PW  -  PULSE  WIDTH 

SAMPLE-  #  OF  SAMPLES  OF  WAVEFORM  PLUS  THE  ZERO-FILL  SAMPLES  NEEDED 
TO  PERFORM  LINEAR  CONVOLUTION 
NPRT  -  •  OF  REPETITIONS  OF  PULSE  MADE  IN  WAVEFORM 

N  -  #  OF  SAMPLES  MADE  OF  PULSE  AFTER  WEIGHTING 

• 

SUBROUTINE  FM ( FREQ.  PERIOD.  NPRT,  PW,  N,  J, SAMPLE) 

INTEGER  PERIOD! 10).  SUM,  SAMPLE,  COUNT, REMAIN 
REAL  FREQU024) 

IF  <  J  EG.  1  >  GO  TO  40 

TYPE  'WAVEFORM  IS  FREQUENCY  MODULATED' 

ACCEPT  'ENTER  BANDWIDTH  TO  BE  SWEPT  '.BAND 
ACCEPT  ENTER  0  IF  LFM.  1  IF  VFM  ', I 
IF  (I  EQ.  1)  GO  TO  19 
WRITE( 10. 300)  BAND 
SUM-0 

DO  10  1-1.  NPRT 
COUNT-SUM*! 

REMAIN-COUNT+N-1 
DO  5  V-COUNT,  REMAIN 
TIME— (PW/N) *<  U-SUH) 

FREQ ( V ) =BAND*TIME*T IME*3  I416/PW 
5  CONTINUE 

SUM-SUM+PER I 00 ( I ) 
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10  CONTINUE 

GO  TO  50 

15  WRITE! 10, 305)  BAND 

SUM=0 

DO  30  1=1. NPRT 

C0UNT=SUM+1 

REMA1N=C0UNT»N-1 

DO  20  J=COUNT, REMAIN/2 

TIME=PW/N*(N-J*2-SUM> 

FREQ ( J ) =-0  5*8AND*T I ME*T 1  ME *3  1416/PW 
20  CONTINUE 

DO  25  J=REMA I N/2+1,  REMAIN 
TIME=PW/N*< J*2-N-SUM) 

FREQ ( J ) =0-  S*BAND*T!ME*TIM£*3  1416/PW 
25  CONTINUE 

SUM=SUM-<-PER10D( I ) 

30  CONTINUE 

GO  TO  50 

40  TYPE  'WAVEFORM  IS  NOT  FREQUENCY  MODULATED' 

DO  45  1  =  1,  SAMPLE 
FREQ ( I >=0  O 
45  CONTINUE 

50  RETURN 

300  FORMAT ( 'LFM  USED.  BANDWIDTH  TO  BE  SWEPT= ' , F 1 2. 3 ) 

305  FORMAT! 'VFM  USED,  BANDWIDTH  TO  BE  SWEPT® ' , F 1 2  3 > 

END 


I 


r>  n  r>  r>  o 


♦**#******##**#***#**»****##************ 

*  subroutine  gpuls  * 

»  WRITTEN  BY.  CART  JOHN  REED  * 

*  AFIT  BOX  4027  * 

»»**»**#•»»*#«*****»**#*»*•**#************ 

C  GPULS  IS  A  FORT0 AN  V  SUBROUTINE  WHICH  DETERMINES  THE  AMPLITUDE 

C  OF  EACH  SAMPLE  OF  THE  RADAR  WAVEFORM  IF  THE  PULSEIS)  IS  GAUSSIAN  SHAPED 

c  ‘Except  for  phase-shift  coded  signals). 
c 

C  THE  GPULS  SUBROUTINE  MAY  BE  CALLED  BY  THE  FOLLOWING  PROGRAMS 
C  MAIN,  GMAIN,  IMAIN 

C  GPULS'  RECEIVES  THE  FOLLOWING  VARIABLES  FROM  THE  CALLING  PROGRAM: 

C  PERIOD,  PW,  NPRT,  N. 

C  IT  RETURNS  THE  FOLLOWING  VARIABLES  TO  THE  CALLING  PROGRAM: 

C  Z ,  GFACTOR. 

C 

C  GPULS'  IS  COMPILED  ON  THE  ECLIPSE  USING  THE  COMMAND 

C  FORTRAN  GPULS 

C 

C  THIS  SUBROUTINE  IS  STORED  IN  FILE  GPULS  FR'  IN  DIRECTORY  DP 4  REED 
C 

C  THE  FOLLOWING  VARIABLES  ARE  USED  IN  THIS  SUBROUTINE: 

C  MEAN  -  1/2  LENGTH  OF  PW 

C  STDEV  -  1/3  OF  MEANi  USE  3  STANDARD  DEVIATIONS  FROM  MEAN 

C  VAR  -  VARIANCE.  STDEV*STDEV 

C  GAUSS  -  SCALING  VALUE  OF  GAUSSIAN  PULSE 

C  GFAC TOR-SCALING  VALUE  for  the  magnitude  of  the  ambicutiy  fcn 

C  NPRT  -  *  OF  PRTS  IN  WAVEFORM 

C  PW  -  PULSE  WIDTH 

C  N  -  #  OF  SAMPLES  MADE  OF  PULSE  AFTER  WEIGHTING 

C  Z  -  ARRAY  CONTAINING  THE  AMPLITUDE  OF  EACH  SAMPLE  OF  WAVEFORM 

C  PERIOD-  ARRAY  CONTAINING  «  OF  SAMPLES  OF  EACH  PRT  AFTER  WEIGHTING 

SUBROUTINE  GPULSIPERIOD.  Z,  PW.  GFACTOR,  NPRT. N > 

COMPLEX  Z ( 1 024 ) 

INTEGER  SUM.  COUNT,  REMAIN, PERIOD! 10) 

REAL  MEAN 

TYPE  ' PULSE ( S >  IS  GAUSSIAN  SHAPED' 

MEAN-0  5*PW 

STDEV=MEAN/3. 0 

VAR=STOEV«STDEV 

GAUSS-1  0/SQRTI2  0*VAR*3  1416) 

GFACTOR-O.  5R235Bi*PW*PW 

WRITE! 10, 300)  MEAN,  STDEV,  VAR,  GFACTOR 

SUM=0 

DO  15  I-l, NPRT 
COUNT  =  SUM+ 1 
REMAIN* COUNT  +  N- 1 
DO  5  J-COUNT, REMAIN 
T I  ME* ' Pw/N ) *  <  J-SUM ) 

P0WER*T1ME-MEAN 
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C ****************************** ********** 

C*  SUBROUTINE  RECPL  * 

C*  WRITTEN  BY  CART  JOHN  REED  • 

C*  15  APRIL  t  R82  * 

C ***************** *********************** 

C  RECPL  IS  A  FORTRAN  V  SUBROUTINE  WHICH  DETERMINES  THE  AMPLITUDE 
C  OF  EACH  SAMPLE  OF  THE  RADAR  WAVEFORM  IF  THE  PULSE(S)  IS  RECTANGULAR 
C  SHAPED  (EXCEPT  FOR  PHASE-SHIFT  CODED  SIGNALS) 

C 

C  THE  RECPL  SUBROUTINE  MAY  BE  CALLED  BY  THE  FOLLOWING  PROGRAMS 
C  MAIN, GMAIN, IMAIN 

C  RECPL'  RECEIVES  THE  FOLLOWING  VARIABLES  FROM  THE  CALLING  PROGRAM: 

C  PERIOD.  NPRT,  N 

C  IT  RETURNS  THE  FOLLOWING  VARIABLES  TO  THE  CALLING  PROGRAM: 

C  Z,  GFACTOR. 

C 

C  RECPL'  IS  COMPILED  ON  THE  ECLIPSE  USING  THE  COMMAND: 

C  FORTRAN  RECPL 

C 

C  this  SUBROUTINE  IS  STORED  IN  FILE  RECPL  FR '  IN  DIRECTORY  DP 4  REED. 

C 

C  the  FOLLOWING  VARIABLES  ARE  USED  IN  THIS  SUBROUTINE 
C  GFACTOR -SC AUNG  VALUE  FOR  THE  MAGNITUDE  OF  THE  AMBIGUITY  FCN 

C  NPRT  *  OF  PRTS  IN  THE  WAVEFORM 

C  Z  -  ARRAY  CONTAINING  THE  AMPLITUDE  OF  EACH  SAMPLE  OF  THE  WAVEFORM 

C  PERIOD-  ARRAY  CONTAINING  THE  «  QF  SAMPLES  OF  EACH  PRT  AFTER  WEIGHTING 

C  N  -  *  OF  SAMPLES  OF  PULSE  AFTER  WEIGHTING 

SUBROUTINE  RECPL ( PERIOD, Z, GFACTOR. NPRT, N) 

COMPLEX  Z ( 1024) 

INTEGER  SUM, COUNT, REMAIN, PERIOD ( 10 ) 

SUM=0 

GFACT0R=1  0 

TYPE  ' PULSE < S I  IS  RECTANGULAR  SHAPED' 

DO  15  1=1, NPRT 
C0UNT=SUM*1 
REMAIN=C0UNT*N-1 
DO  5  J=COUNT, REMAIN 
Z ( J)=( 1  0, 0  0) 

5  CONTINUE 

DO  10  J=REMAIN* 1 . SUM+PER I OD < I ) 

Z<  J)»<C  0.  0  0) 

10  CONTINUE 

SUM«SUM*PERIOD(  I ) 

15  CONTINUE 

RETURN 
END 
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Z<J)=CMPLX(GAUSS*£XP(-POWER*POWER/<2  0»VAR)),0  0) 

5  CONTINUE 

DO  10  J=REMAIN-M  <  SUM+PER 1 0D< I ) 

U J>  =  <0  0, 0  0) 

10  CONTINUE 

SUM=SUM*PERIODCI) 

15  CONTINUE 

RETURN 

300  FORMAT!  'M= F15  1  0,  3X,  'SD* F12  1 0,  3X,  '  VAR-  FI 2  1 0,  3X,  ' GFAC=  FI 2  10 > 

END 
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c  »»*****»«*»**«* •  «#*********»***■»  ******** 


c* 

SUBROUTINE  SETUP 

* 

c* 

WRITTEN  BY  CAPT  JOHN  REED 

# 

c* 

1SAPR1L  1RS2 

« 

C **************************** ************ 

C  'SETUP'  IS  A  FORTRAN  7  SUBROUTINE  WHICH  SETS  UP  ARRAYS  WITH  THE 
C  TIME-REVERSED  CONJUGATE  OF  THE  COMPLEX  ENVELOPE  (XVALS)  AND  THE  DELAYED, 

C  DOPPLER-SHIFTED  VALUES  OF  THE  COMPLEX  ENVELOPE ( YVALS >  THIS  WILL 
C  THEN  ALLOW  US  TO  PERFORM  THE  LINEAR  CONVOLUTION  LATER  ON 
C 

C  THIS  SUBROUTINE  IS  CALLED  BY  THE  'MAIN'  PROGRAM  AND  THE  PHOTO' 

C  AND  'INTEG'  SUBROUTINES. 

C  this  SUBROUTINE  RECEIVES  THE  FOLLOWING  VARIABLES  FROM  THE  CALLING  ROUTINE: 
C  Z,  FREQ- VALUE,  TOTAL.  SAMPLE,  M. 

c  this  subroutine  returns  the  following  variables  to  the  calling  routine: 

C  XVALS,  YVALS. 

C  SETUP'  IS  COMPILED  ON  THE  ECLIPSE  USING  THE  COMMAND. 

C  FORTRAN  SETUP 

C 

C  'SETUP'  IS  STORED  IN  FILE  'SETUP. FR'  IN  DIRECTORY  DPR  REED 

C 

C  THE  FOLLOWING  VARIABLES  ARE  USED  IN  THIS  SUBROUTINE: 

C  M  -  #  OF  NON-ZERO-FILL  SAMPLES  OF  WAVEFORM  AFTER  WEIGHTING 

C  SAMPLE-  *  OF  SAMPLES  OF  WAVEFORM  WITH  ALL  ZERO-FILL  SAMPLES 

c  YVALS  -  ARRAY  CONTAINING  THE  DELAYED, DOPPLER-SHI FTED  SAMPLES  OF  THE 

C  COMPLEX  ENVELOPE  OF  THE  WAVEFORM 

C  XVALS  -  ARRAY  CONTAINING  THE  TIME-REVERSED,  CONJUGATE  SAMPLES  OF  THE 

C  COMPLEX  ENVELOPE  OF  THE  WAVEFORM 

C  FREQ  -  ARRAY  CONTAINING  THE  FM  OF  EACH  SAMPLE 

C  TIME  -  VALUE  OF  TAU  (DELAY) 

C  '  TOTAL  -  TOTAL  TIME  DURATION  OF  THE  WAVEFORM 

C  VALUE  -  IMAGINARY  VALUE  »  DOPPLER  FREQUENCY ( W ) 

C  Z  -  ARRAY  CONTAINING  THE  AMPLITUDE  DF  EACH  SAMPLE  OF  WAVEFORM 

SUBROUTINE  SETUP (Z, FREQ, XVALS. YVALS, VALUE,  TOTAL.  SAMPLE,  M) 

COMPLEX  XVALS ( 1024). YVALS < 1024),  Z ( 1024), VALUE, AMT  1 , AMT2 
INTEGER  SAMPLE 
REAL  FREQ(1024) 

DO  5  I =M+1. SAMPLE 
YVALS < I )*(0  0,  0  0) 

5  CONTINUE 

DO  10  1=1. SAMPLE-M 
XVALS ( I  )*<0  0,  0  0) 

10  continue 

DO  15  1  =  1.  M 

AMT2«CMPLX(0  0,  FREQ< I ) ) 
time=total/m#i 

YVALS ( I )»Z< I )*CEXP< VALUE*TIME*AMT2) 

15  CONTINUE 

DO  20  I-SAMPLE-M+1. SAMPLE 
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AMT  1=2  <  I  -SAMPLE-1-!"!) 

AMT2=CMPLX (0.  0,  -FREQ  <  I -SAMPLE-"-M  > ) 
XVALS< I ) =AMT 1 *CEXP ( AMT2 ) 

CONTINUE 

RETURN 

END 


1 


C«*ft»***«*»*****»*#***««**«»«*»««*«W«4«** 

C*  SUBROUTINE  OUT  * 

C*  WRITTEN  BY'  CAPT  JOHN  REED  » 

C*  15  APRIL  1982  * 

C**#********#**#*************#***#-*#***** 

C  'OUT'  IS  A  FORTRAN  V  SUBROUTINE  WHICH  DISPLAYS  THE  MAGNITUDE 
C  OF  THE  AMBIGUITY  FUNCTION  AT  VARIOUS  DELAY  TIMES  (VALUES  OF  TAU). 
C 

C  OUT'  IS  CALLED  BY  THE  'MAIN'  PROGRAM. 

C 

C  OUT'  RECEIVES  THE  FOLLOWING  VARIABLES  FROM  'MAIN' 

C  MAG,  NUM,  PW,  TOTAL.  M. 

C  OUT'  DOES  NOT  RETURN  ANY  VARIABLE  TO  MAIN'. 

C 

C  OUT'  IS  COMPILED  ON  THE  ECLIPSE  USING  THE  COMMAND: 

C  FORTRAN  OUT 

C 

C  OUT'  IS  STORED  IN  FILE  OUT  FR '  IN  DIRECTORY  DP4  REED 

C 

C  THE  FOLLOWING  VARIABLES  ARE  USED  IN  'OUT': 

C  PW  -  PULSE  WIDTH 

C  NUM  -  *  OF  SAMPLES  MADE  OF  WAVEFORM 

C  M  -  *  OF  SAMPLES  MADE  OF  WAVEFORM  AFTER  WEIGHTING 

C  TOTAL  -  TOTAL  TIME  DURATION  OF  WAVEFORM 

C  MAG  -  ARRAY  CONTAINING  THE  MAGNITUDE  OF  EACH  SAMPLE  OF  THE 

C  AMBIGUITY  FCN  CUT 

SUBROUTINE  OUTIMAG. NUM, PW, TOTAL, M) 

REAL  MAG  < 1024) 

WRITE! 10.  300)  PW, NUM 
WRITE! 10,  305) 

DO  5  1  =  1,  M 

T I  ME- <  TOTAL/M*! 1-1 ) )- TOTAL 
WRITE! 10.  310)  TIME  . MAG ( I ) 

5  CONTINUE 

DO  10  I=M+1,M*2 
T I ME=TOTAL/M* ( I-M-l ) 

WRITE! 10.  310)  T I  ME ,  MAG  < I ) 

10  CONTINUE 

RETURN 

300  FORMAT!  PW-'. F 12  10.  5X.  '*  OF  SAMPLES  MADE  OF  WAVEFORM- 1 5 ) 

305  FORMAT!  'TIME(DELAY)  5X,  MAGNITUDE  OF  AMBIGUITY  FUNCTION') 

310  FORMAT <  F  1 2  10,  7X.F10  6) 

END 


< 


c  **********  **«•*  *****#**♦**■***##■«■***  ***♦*# 


c* 

PRGRAM  GMAIN 

* 

c# 

WRITTEN  BY.  CAPT  JOHN  REED 

* 

c» 

15  APRIL  1982 

* 

c  *•****♦■*♦*  ***********  ********  ************ 


C  GMAIN '  IS  A  FORTRAN  V  PROGRAM  WHICH  CALCULATES  THE  MAGNITUDE  OF  THE 

C  AMBIGUITY  FUNCTION  FOR  A  SPECIFIED  RADAR  WAVEFORM.  THE  MAGNITUDE  IS 
C  CALCULATED  FOR  64  DIFFERENT  DOPPLER  FREQUENCIES,  BEGINNING  AND  ENDING 
C  WITH  FREQUENCIES  THE  USER  SPECIFIES  AND  AT  EQUAL  INTERVALS  IN  BETWEEN. 
C  THESE  MAGNITUDES  ARE  QUANTIZED  INTO  VALUES  RANGING  FROM  0  TO  200. 

C  THEN,  UP  TO  64  SAMPLES  OF  EACH  OF  THESE  QUANTIZED  MAGNITUDES  APE 
C  WRITTEN  INTO  A  FILE  'DATA'  THESE  DATA  ARE  THEN  USED  BY  ANOTHER 
C  PROGRAM  PLTTRNS'  TO  PROVIDE  A  THREE-DIMENSIONAL  PLOT  OF  THE  MAGNITUDE 
C  OF  THE  AMBIGUITY  FUNCTION 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c- 

c 

c 

c 

c 

c 

c 

c 


THE  FOLLOWING  ARE  USER  SELECTED  PARAMETERS  WHICH  THE  PROGRAM  (OR 
THE  SUBROUTINES  WHICH  IT  CALLS)  WILL  REQUEST: 

1  PULSE  WIDTH 

2  WHETHER  OR  NOT  THE  PULSE  IS  REPEATED. 

A  WHETHER  PULSE  REPETITION  IS  STAGGERED  OR  THE  SAME. 

B.  PULSE  REPETITION  TIME(S)  (PRT). 

C  NUMBER  OF  REPETITIONS 

3  NUMBER  OF  SAMPLES  TO  BE  MADE  OF  WAVEFORM. 

4  UHETHER  OR  NOT  THE  PULSE(S)  IS  PHASE-SHIFT  CODED 

A  NUMBER  OF  ELEMENTS  IN  THE  PULSE  CODE 
B  PHASE  OF  EACH  ELEMENT 

5  WHETHER  OR  NOT  THE  PULSE < S )  IS  FREQUENCY  MODULATED 

A.  BANDWIDTH  TO  BE  SWEPT  BY  THE  FM 
B  IF  LINEAR  FM  OR  VFM 

6  WHETHER  THE  PULSE ( S )  IS  RECTANGULAR  OR  GAUSSIAN  5HAPED. 

7  STARTING  VALUE  OF  DELAY  <TAU>  FOR  THE  SAMPLES  OF  THE  MAGNITUDE 
TO  BE  WRITTEN  INTO  THE  DATA'  FILE. 

8  WHETHER  OR  NOT  TO  SKIP  SAMPLES,  AND  IF  SO,  HOW  MANY. 

9  STARTING  t<  FINAL  VALUES  OF  THE  DOPPLER  FREQUENCIES  FOR  WHICH 
MAGNITUDES  ARE  CALCULATED 

GMAIN'  MAY  MAKE  CALLS  TO  THE  FOLLOWING  SUBROUTINES. 

INPUT,  PHASE,  FM  GPULS,  RECPL,  PHQTOtWHICH  CALLS  THE  FAST  FOURIER 
TRANSFORM  SUBROUTINE  < EDFT  LB)  FROM  THE  LIBRARY). 

GMAIN'  IS  COMPILED  AND  LOADED  ON  THE  ECLIPSE  USING  THE  FOLLOWING 
COMMANDS: 

FORTRAN  GMAIN 

RLDR/N  GMAIN  INPUT  PHASE  FM  GPULS  RECPL  PHOTO  EDFT  LB  SFLIBS 

THIS  PROGRAM  IS  STORED  IN  FILE  'GMAIN  FR '  IN  DIRECTORY  DP4  REED. 

THE  FOLLOWING  VARIABLES  ARE  USED  IN  'GMAIN' 

PW  -  PULSE  WIDTH 

SAME  -  PRT  IF  PULSES  ARE  REPEATED  AT  SAME  INTERVAL 

NPRT  -  *  OF  REPETITIONS  OF  PULSE  TO  BE  MADE 

pRT  -  ARRAY  CONTAINING  THE  PRT  FDR  EACH  PULSE  REPETITION 

TOTAL  -  TOTAL  TIME  DURATION  OF  WAVEFORM 
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C  MUM  -  #  OF  SAMPLES  MADE  OF  WAVEFORM 

SAMPLE-  #  OF  SAMPLES  OF  WAVEFORM  PLUS  ZERO-FILLING  SAMPLES  NEEDED 
C  TO  PERFORM  LINEAR  CONVOLUTION  <=2*NUM> 

C  PERIOD-  ARRAY  CONTAINING  #  OF  SAMPLES  MADE  OF  EACH  PRT 

C  N  -  #  OF  SAMPLES  MADE  OF  rULSE  AFTER  WEIGHTING 

C  M  -  #  OF  NON-ZERO-FILL  SAMPLES  MADE  OF  WAVEFORM  AFTER  WEIGHTING 

C  Z  -  ARRAY  CONTAINING  THE  AMPLITUDE  OF  EACH  SAMPLE 

C  GFACTOR-SCALING  VALUE  OF  magnitude 

C  FREQ  -  ARRAY  CONTAINING  THE  FM  FOR  EACH  SAMPLE 

C  SCALE  -  SCALING  FACTOR  OF  MAGNITUDE 


COMPLEX  Z ( 1 024 ) 

INTEGER  SAMPLE, PERIOD< 10) 

REAL  PRT< 10),  FREG( 1024) 

TYPE  '  WELCOME  TO  GMAIN' 

TYPE  THIS  PROGRAM  IS  DESIGNED  TO  COMPUTE  THE  MAGNITUDE  OF  64' 

TYPE  'DIFFERENT  DOPPLER  CUTS  OF  THE  AMBIGUTIY  FUNCTION,  ' 

TYPE  'QUANTIZE  THE  MAGNITUDES,  AND  OUTPUT  UP  TO  64  SAMPLES  OF' 

TYPE  'EACH  CUT  TO  FILE  DATA.  YOU  WILL  BE  REQUIRED  TO  ENTER' 

TYPE  PARAMETERS  AT  VARIOUS  TIMES  DURING  THE  PROGRAM.  WHEN  THE' 

TYPE  PROGRAM  REQUESTS  AN  ENTRY,  SIMPLY  TYPE  YOUR  REPLY  AND, ' 

TYPE  'IF  CORRECT,  HIT  CR  KEY  USE  THE  <DEL>  KEY  TO  ERASE  CHARACTERS' 
TYPE  AND  THE  C\>  TO  ERASE  THE  ENTIRE  ENTRY  YOU  MAY  STDP  THE' 

TYPE  PROGRAM  AT  ANY  TIME  BY  ENTERING  <'A>.  YOU  MAY  STDP  THE  ' 

TYPE  SCREEN  OUTPUT  USING  OS>  AND  RESUME  USING  <~Q>  ' 

PAUSE  ENTER  ANY  CHARACTER  TO  BEGIN' 

CALL  I NPUT (PERIOD.  PRT,  PW,  NPRT ,  TOTAL-  SAMPLE,  NUM,  M,  N) 

ACCEPT  ENTER  0  IF  RULSE<S)  IS  PHASE-SHIFT  CODED.  1  IF  NOT  ',J 
IF  (J  EO  II  GO  TO  5 

CALL  PHASE ( FREQ,  PRT,  PERIOD,  Z,  M,  N, NPRT,  GF ACTOR,  PW,  SAMPLE,  TOTAL ) 

GO  TO  20 

5  TYPE  PULSE  <  S )  IS  NOT  PHASE-SHIFT  CODED' 

ACCEPT  ENTER  0  IF  PULSE<S)  IS  FREQUENCY  MODULATED,  1  IF  NOT  ',J 
CALL  FM < FREQ.  PERIOD,  NPRT,  PW.  N,  U,  SAMPLE) 

10  ACCEPT  'ENTER  0  FOR  RECTANGULAR-SHAPED  PUL5E(S).  1  FOR  GAUSSIAN  '.U 

IF  IJ  EQ.OI  GO  TO  15 

CALL  GPULS ( PER  I OD,  Z. PW- GF ACTOR, NPRT, N) 

GO  TO  20 

15  CALL  RECPL (PERIOD,  Z, GF ACTOR,  NPRT,  N) 

20  SCALE=GFACTOR / ( N»NPRT ) 

CALL  PHOTO <  Z . FREQ,  TOTAL, M, SAMPLE- SCALE) 

TYPE  "DATA  HAVE  BEEN  TRANSFERRED  TO  FILE  'DATA'" 

TYPE  'AFTER  THIS  PROGRAM  ENDS,  ENTER  THE  FOLLOWING  COMMAND  ON' 

TYPE  'THE  TEKTRONIX  SCOPE  PLTTRNS  DATA/I  (ARRAY  SIZE)/N' 

END 


1 

4 


(;#***-**♦*.****##**#**#*****#***#*■##*****## 
c*  subroutine  photo  * 

C*  WRITTEN  BY  CAPT  JOHN  REED  * 

C*  15  APRIL  1982  * 

C**************************************** 


C  PHOTO'  IS  A  FORTRAN  V  SUBROUTINE  WHICH  TAKES  A  DESIRED  RADAR 

C  WAVEFORM,  PERFORMS  THE  CALCULATIONS  TO  PRODUCE  THE  AMBIGUITY  FCN 
C  MAGNITUDES  AT  VARIOUS  CUTS  OF  DOPPLER  FREQUENCIES,  QUANTIZES  THE 
C  MAGNITUDES  TO  VALUES  BETWEEN  0  AND  200.  SELECTS  THE  SAMPLES  OF  THE 
C  MAGNITUDES  TO  BE  OBSERVED.  AND  WRITES  THESE  DATA  INTO  FILE  'DATA'. 

C 

C  PHOTO'  IS  CALLED  BY  THE  PROGRAM  ' GMA IN'. 

C  PHOTO'  MAKES  CALLS  TO  THE  'SETUP'  SUBROUTINE  AND  TO 

C  THE  FFT  LIBRARY  SUBROUTINE  EDFT  LB  ( DFT5 ) . 

C  -PHOTO'  RECEIVES  THE  FOLLOWING  VARIABLES  WHEN  CALLED  BY  'GMAIN': 

C  Z,  FREQ,  TOTAL.  M.  SAMPLE.  SCALE. 

C  PHOTO'  DOES  NOT  RETURN  ANY  VARIABLES  TO  'GMAIN'. 

C 

C  -PHOTO'  IS  COMPILED  ON  THE  ECLIPSE  USING  THE  COMMAND: 

C  FORTRAN  PHOTO 

C  PHOTO'  IS  STORED  IN  FILE  PHOTO  FR '  IN  DIRECTORY  DP4  REED 

C 

C  THE  FOLLOWING  VARIABLES  ARE  USED  IN  'PHOTO': 

C  INTERVAL-TIME  (TAU)  BETWEEN  EACH  SAMPLE  OF  WAVEFORM 

C  M  -  *  OF  NON-ZERO-FILL  SAMPLES  MADE  OF  WAVEFORM  AFTER  WEIGHTING 

C  JL  -  *  OF  SAMPLES  OF  MAGNITUDE  OF  AMB  FCN  FOR  EACH  DOPPLER 

C  TOTAL  -  TOTAL  TIME  DURATION  OF  WAVEFORM 

C  A  -  STARTING  VALUE  OF  TIME  FOR  MAGNITUDE  OF  AMB  FCN  FOR  EACH  CUT 

C  N  -  HOW  TO  SAMPLE  AVAILABLE  SAMPLES  OF  MAGNITUDE 

c  start  -  starting  SAMPLE  TO  BE  OBSERVED 

C  STOP  -  FINAL  SAMPLE  TO  BE  OBSERVED 

C  WS  -  STARTINC  DOPPLER  FREQUENCY  TO  BE  CONSIDERED 

C  WF  -  FINAL  DOPPLER  FREQUENCY  TO  BE  CONSIDERED 

C  WIN  -  INTERVAL  BETWEEN  DOPPLER  FREQUENCIES  TO  BE  CONSIDERED 

c  W  -  ARRAY  CONTAINING  ALL  DOPPLER  FREQUENCIES  TO  BE  CONSIDERED 

C  IER  -  ERROR  LABEL  FOR  A  FILE  OPERATION  (OKAY  IF  »  1) 

C  Z  -  ARRAY  CONTAINING  THE  AMPLITUDE  FOR  ALL  SAMPLES  OF  WAVEFORM 

C  FREQ  -  ARRAY  CONTAINING  THE  FMOF  ALL  WAVEFORM  SAMPLES 

C  X VALS  -  ARRAY  CONTAINING  THE  T I ME-REVERSED, CONJUGATE  SAMPLES  OF 

C  THE  COMPLEX  ENVELOPE  OF  THE  WAVEFORM.  ALSO  CONTAINS  THE 

C  AMBIGUITY  FCN  VALUES  LATER  IN  THE  ROUTINE 

C  YVALS  -  ARRAY  CONTAINING  THE  DOPPLER  AND  RANGE -DEL AYED  SAMPLES  OF 

c  THE  COMPLEX  ENVELOPE  OF  WAVEFORM 

C  VALUE  -  IMAGINARY  VALUE  *  W 

C  SAMPLE-  *  OF  SAMPLES  OF  WAVEFORM  PLUS  ZERO-FILLINO  SAMPLES  NEEDED 

C  TO  PERFORM  LINEAR  CONVOLUTION 

C  MAG  -  ARRAY  CONTAINING  THE  SAMPLES  OF  THE  MAGNITUDE  OF  THE 

C  AMBIGUITY  FCN  FOR  A  DOPPLER  CUT 

C  G  -  AN  UNQUANT  I  ZED  VALUE  OF  MAGNITUDE 

c  K  -  A  QUANTIZED  VALUE  OF  MAGNITUDE  (BETWEEN  C  AND  200) 

SUBROUTINE  PHOTO(Z.  FREQ.  TOTAL, M,  SAMPLE,  SCALE) 

COMPLEX  XVALS( 1C24> . YVALS ( 1024) ,  Z ( 1024) ,  VALUE, STORE (64) 


( 
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INTEGER  DIFF, A,  START.  STOP.  GRAPH<64,  64).  SAMPLE 
REAL  INTERVAL. W<64 ) . FREQi 1024) ,  MAG < 1024) 

interval=total/m 

MR  I T£ ( 10. 300)  -TOTAL. TOTAL.  INTERVAL 

ACCEPT  'CENTER  VALUE  OF  TAU  IS  A*  ( TOT AL/ M) .  ENTER  A  ' ,  A 

ACCEPT  'TO  GRAPH  EACH  SAMPLE,  ENTER  1,  EVERY  OTHER  SAMPLE. 2.  ETC  N 

WRITE  (10. 305)  A* INTERVAL 

NUM=SAMPLE/2 

oiff=num-m 

IF  (  (A+NUf1-32*N)  GT  0)  GO  TO  1 
START=A- ( 32  *N ) *NUM 
STOP=A* ( 32*N) *NUH 
GO  TO  2 

1  START=A-(32*N)+M 
ST0P=A*(32*N)+M 

2  IF  (START  LT  0)  START=0 

IF  (STOP  GT  SAMPLE)  STOP=SAMPLE 
WRITE! 10. 310)  <START-M)«INTERVAL 
WRITE( 10. 315)  <STQP-NUM+DIFF)*INTERVAL 
UL= (STOP-START) /N 

IF  (JL  NE.  64)  TYPE  'ERROR  IN  SAMPLING  SELECTION' 

ACCEPT  ENTER  STARTING  VALUE  OF  W  ' , WS 
ACCEPT  'ENTER  FINAL  VALUE  OF  U  WF 
WIN=!WF-US)/64  O 
WRITE! 10. 320)  (WF*WS)/20 
DO  5  1  =  1,  64 
W(  I  >=WS-M4IN*I 
5  CONTINUE 

CALL  OFILW! "DATA", IER) 

TYPE  'FOR  DFILW.  IER='.IER 
CALL  CFILU< "DATA", 2, IER) 

TYPE  'FOR  CFILW  ,  IER= IER 
CALL  OPEN  <6.  "DATA", 2,  2ER> 

TYPE  'FOR  OPEN,  IER='.1ER 
DO  25  1  =  1.  64 
WRITE! 10.  325)  W(J) 

VALUE=CMPLX!0  0.  W( I ) ) 

CALL  SETUP 1Z,  FREO,  XVALS, YVALS. VALUE, TOTAL, SAMPLE,  M) 

CALL  DFT5< XVALS. SAMPLE, 0) 

CALL  DFT5(YVALS.  SAMPLE, 0) 

DO  10  U=l. SAMPLE 
XVALS! U)=CONUG<  XVALS! J> > 

YVALS! U)=CONVG( YVALS! U) > 

XVALS ! U ) =X VALS1 J ) *YVALS ( J  > /SAMPLE 
10  CONTINUE 

CALL  DFT5! XVALS,  SAMPLE,  J ) 
do  is  j= i , sample 

MAG! J)=CABS( XVALS! J) >*SCALE 
15  CONTINUE 

WRITE! 10. 330)  W( I ) 

DO  20  U=START*N, STOP. N 
G- MAG  (  -J  ) 

F=0 


( 
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DO  21  L= 1 . 2C0 
R=0  OC5*L 
IF  <G  GE  R)  K=L 
21  CONTINUE 

GRAPH! IJ-5TART) /N, I  )  =K 
20  CONTINUE 

25  CONTINUE 

DO  35  1=1. 64 
DO  30  J=l- JL 
STORE ( J ) =GRAPH  <  J.  I ) 

30  CONTINUE 

CALL  WRBLH  <6.  I . STORE.  1 .  I ER ) 

35  CONTINUE 

CALL  CLOSE ! 6.  IER 1 
TYPE  FOR  CLOSE,  I£R=',IER 
WRITE!  10.  310)  (  ST  ART-M  )  *  INTERVAL 
WRITE! 10. 315)  !STOP-NUM+DIFF)*INTERVAL 
‘TYPE  'ARRAY  5IZE='.UL 
RETURN 

300  FORMAT (  ' T AU  GOES  FROM  '.F12  10.  'TO  '.F12  10.'  IN  INTERVALS  OF  '.F12  10) 

305  FORMAT! 'TAU  IS  CENTERED  AT  '.F12  10) 

310  FORMAT! 'TAU  STARTS  AT  '.F12  9) 

315  FORMAT! 'TAU  STOPS  AT  '.F12  9) 

320  FORMAT!  'W  IS  CENETERED  AT  '.F12  2) 

325  FORMAT!  'W  NOW  BEING  EVALUATED  =',F12.  2) 

330  FORMAT!  FOR  W  =  ' ,  F  1 2.  2.  3X,  MAGNITUDES  HAVE  BEEN  CALCULATED') 

END 


( 


144 


o  n  n 


c «••*••••***««* •**«-«*• ********* 

C*  PROGRAM  IMAIN  * 

*  WRITTEN  DY:  CAPT  JOHN  REED  * 

*  15  APRIL  1982  * 

**•****«****«*#***•»•###*****#**#**♦*»»*## 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


'IMAIN'  IS  A  FORTRAN  V  PROGRAM  WHICH  CALCULTES  THE  S I GNAL-TO-CLUTTER 
RATIO  FOR  A  SPECIFIED  RADAR  WAVEFORM,  CLUTTER  DISTRIBUTION,  AND 
DESIRED  TARGET.  APPROXIMATIONS  ARE  MADE  IN  CALCULATING  THE  VOLUME  OF 
THE  PRODUCT  OF  THE  AMBIGUITY  FCN  MAGNITUDE  AND  THE  CLUTTER  PROBABILITY 
DENSITY  FCN  IN  THIS  PROGRAM  THE  STANDARD  PROBABILTI Y  DENSITY  FCNS  FOR 
CLUTTER  ARE  UNIFORM  AND  GAUSSIAN.  THE  PROGRAM  ALLOWS  THE  USER  TO 
INPUT  ANY  OTHER  NON-STANDARD  DENSITY  FOR  CLUTTER. 

THE  FOLLOWING  ARE  USER  SELECTED  PARAMETERS  WHICH  THE  PROGRAM 
(OR  THE  SUBROUTINES  IT  CALLS)  WILL  REQUEST: 

1  PULSE  WIDTH. 

2.  WHETHER  OR  NOT  THE  PULSE  IS  REPEATED. 

A.  WHETHER  PULSE  REPETITION  IS  STAGGERED  OR  THE  SAME. 

B  PULSE  REPETITON  TIME(S)  <PRT>. 

C.  NUMBER  OF  REPETITIONS. 

4  WHETHER  OR  NOT  THE  PULSEIS)  IS  PHASE-SHIFT  CODED 

A  NUMBER  OF  ELEMENTS  IN  THE  PULSE  CODE 
B  PHASE  OF  EACH  ELEMENT. 

5  WHETHER  OR  NOT  THE  PULSE(S)  IS  FREQUENCY  MODULATED. 

A  BANDWIDTH  TO  BE  SWEPT  BY  THE  FM 

B.  IF  LINEAR  FM  OR  VFM. 

A  WHETHER  PULSE<S)  IS  RECTANGULAR  OR  GAUSSIAN  SHAPED 
7  STARTING  DOPPLER  FREQUENCY  FOR  CLUTTER  S.  AMBIGUITY  FCN 
a  FINAL  DOPPLER  FREQUENCY  FOR  CLUTTER  &  AMBIGUITY  FCN 

9  INTERVAL  BETWEEN  DOPPLER  FREQUENCIES. 

10  PROBABILITY  DENSITY  FCN  FOR  THE  CLUTTER 
11.  BOUNDS  ON  DELAY  < TAU )  OF  CLUTTER 

12  IF  CLUTTER  IS  NOT  A  STANDARD  PROB  DENSITY  FCN,  THEN 

A.  *  OF  AREAS  OF  CLUTTER. 

B  BOUNDS  OF  EACH  CLUTTER  AREA  IN  DELAY ( TAU  >  S,  DOPPLER. 

13  IF  CLUTTER  IS  A  STANDARD  PROB  DENSITY  FCN,  THEN 

A.  NECESSARY  STATISTICS  OF  DOPPLER  TO  DEFINE  DENSITY  FCN 
B  NECESSARY  STATISTICS  OF  DELAY (TAU )  TO  DEFINE  DENSITY  FCN 

IMAIN'  MAY  MAKE  CALLS  TO  THE  FOLLOWING  SUBROUTINES: 

INPUT, FHASE- FM, GPULS, RECPL, INTEGIWHICH  CALLS  SEVERAL  SUBROUTINES) 

IMAIN'  IS  COMPILED  AND  LOADED  ON  THE  ECLIPSE  USING  THE  FOLLOWING 
COMMANDS. 

FORTRAN  IMAIN 

RLDR/N  IMAIN  INPUT  PHASE  FM  GPULS  RECPL  INTEO  IDC  CPDID  SETUP  PROB 
EDFT  LB  9FLIB9 

THIS  PROGRAM  IS  STORED  IN  FILE  'IMAIN FR  '  IN  DIRECTORY  DP4  REED. 

the  FOLLOWING  VARIABLES  ARE  USED  IN  'IMAIN' 

PERIOD  -  ARRAY  CONTAINING  *  OF  SAMPLES  MADE  OF  EACH  PRT 
PRT  -  ARRAY  CONTAINING  THE  PRT  FOR  EACH  REPETITON 
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C  PW  -  PULSE  WIDTH 

C  NPRT  -  «  OF  REPETITIONS  OF  PULSE  TO  BE  MADE 

C  TOTAL  -  TOTAL  TIME  DURATION  OF  WAVEFORM 

C  NUM  -  «  OF  SAMPLES  TO  BE  MADE  OF  WAVEFORM 

C  SAMPLE-  #  OF  SAMPLES  MADE  OF  WAVEFORM  PLUS  ZERO-FILL  SAMPLES  NEEDED 

C  TO  PERFORM  LINEAR  CONVOLUTION  (=2*NUM> 

C  M  -  #  OF  SAMPLES  MADE  OF  WAVEFORM  AFTER  WEIGHTING 

C  N  -  #  OF  SAMPLES  MADE  OF  PULSE  AFTER  WEIGHTING 

C  FREQ  -  ARRAY  CONTAINING  THE  FM  OF  EACH  SAMPLE  OF  WAVEFORM 

C  Z  -  ARRAY  CONTAINING  THE  AMPLITUDE  OF  EACH  SAMPLE  OF  WAVEFORM 

C  CFACTOR-SCALINC  VALUE  OF  MAGNITUDE 

C  SCALE  -  SCALING  FACTOR  OF  MAGNITUDE 

C  VALUE  -  IMAGINARY  VALUE  «  W 

C  TRCS  -  RADAR  CROSS  SECTION  (RCS)  OF  DESIRED  TARGET 

C  CRCS  -  AVERAGE  RCS  OF  THE  CLUTTER 

C  AVGVOL-  AVERAGE  VOLUME  OF  THE  PRODUCT  OF  THE  AMB  FCN  MAG  AND 

C  CLUTTER  P ROB  DENSITY  FCN 

C  NS  -  TOTAL  *  OF  SCATTERERS 

C  SCR  -  SIGNAL-TO-CLUTTER  RATIO 

COMPLEX  Z  < 1 024 ) 

INTEGER  SAMPLE. PERIOD* 10) 

REAL  PRT<10), FREQ< 1024) 

TYPE  '  WELCOME  TO  IMAlN' 

TYPE  'THIS  PROGRAM  IS  DESIGNED  TO  COMPUTE  THE  SIGNAL-TO-CLUTTER ' 

TYPE  RATIO  FOR  A  SPECIFIED  RADAR  WAVEFORM,  CLUTTER  DISTRIBUTION,  ' 
TYPE  AND  TARGET.  YOU  WILL  BE  REQUIRED  TO  ENTER  PARAMETERS  AT' 

TYPE  'VARIOUS  TIMES  DURING  THE  PROGRAM.  WHEN  THE  PROGRAM  REQUESTS' 
TYPE  'AN  ENTRY.  SIMPLY  TYPE  YOUR  REPLY,  AND  IF  CORRECT,  HIT  THE  ' 
TYPE  CR  KEY.  USE  THE  <OEL>  KEY  TO  DELETE  CARACTERS  AND  THE  <\>' 
TYPE  'TO  DELETE  THE  ENTIRE  ENTRY.  USE  <~A>  TO  STOP  THE  PROGRAM.  ' 
TYPE  USE  <~S>  TO  STOP  THE  SCREEN  DISPLAY,  AND  <~Q>  TO  RESUME- ' 

PAUSE  'ENTER  ANY  CHARACTER  TO  BEGIN' 

CALL  INPUT<PERIOD,  PRT , PW  NPRT.  TOTAL,  SAMPLE, NUM, M, N) 

ACCEPT  'ENTER  0  IF  PULSE(S)  IS  PHASE-SHIFT  CODED,  1  IF  NOT  ',J 
IF  <  J  EQ.  1 )  CO  TO  3 

CALL  PH ASE(FREQ,  PRT,  PERIOD,  Z,  M.  N,  NPRT,  GFACTOR,  PW,  SAMPLE, TOTAL) 

GO  TO  20 

5  TYPE  PULSE < S )  IS  NOT  PHASE-SHIFT  CODED' 

ACCEPT  'ENTER  0  IF  PULSE(S)  IS  FREQUENCY  MODULATED,  1  IF  NOT  \J 
CALL  FM ( FREQ. PERIOD, NPRT, PW, N, J, SAMPLE) 

10  ACCEPT  'ENTER  0  FOR  RECTANGULAR-SHAPED  PULSE(S),  1  FOR  GAUSSIAN  ',0 

IF  <U  EQ  0)  GO  TO  15 

CALL  GPULS<PERIOD,  Z,  PW, GFACTOR,  NPRT, N) 

GO  TO  20 

15  CALL  RECPL (PERIOD, Z.  GFACTOR,  NPRT,  N) 

20  SCALE=GF ACTOR/ <  N*NPRT ) 

CALL  INTEG ( Z, FREQ,  TOTAL. SAMPLE.  M,  AVGVOL,  SCALE) 

ACCEPT  'ENTER  THE  RADAR  CROSS  SECTION  (RCS)  OF  DESIRED  TARGET  ',TRCS 
ACCEPT  'ENTER  THE  AVERAGE  RCS  OF  THE  CLUTTER  '.CRCS 
ACCEPT  'ENTER  THE  TOTAL  NUMBER  OF  SCATTERERS  ' , NS 
SCR*TRCS/ (NS*CRCS*AVGVOL) 

WRITEt 10, 300)  TRCS. CRCS 
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WR I  TE  (  1 0,  305 )  NS.  AVCVOL 
UIRITE<  10.  010)  SCR 

300  FORMAT!  'TARCET  RCS«'.F10  b.  5X,  'AVC  RCS  OF  CLUTTER- F 1 0  6) 

305  FORMAT (  TOTAL  #  OF  SCATTERERS- ' .  15.  5X,  'AVC  VOL  OF  1 NTEGR AL= ' . F 10  6) 

310  FORMAT!  SISNAL-TO-Ct UTTER  RATIO- '.FlO  6) 

END 
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SUBROUTINE  jnteg 
WRITTEN  BY  CAPT  JOHN  REED 
15  APRIL  1982 


C  * 
c  * 


C  '  I NTEG '  IS  A  FORTRAN  V  SUBROUTINE  WHICH  TAKES  A  RADAR  WAVEFORM, 

c  takes  the  clutter  distribution,  and  calculates  the  volume  of 
C  the  PRODUCT  OF  THE  MAGNITUDE  OF  THE  AMBIGUT1Y  FCN  AND  THE  CLUTTER 
C  PROBABILITY  DENSITY  FCN  IT  CALCULATES  THE  MAGNITUDES  OF  THE  AMB 
C  FCN  FOR  THE  CUTS  OF  DOPPLER  FREQUENCIES  SPECIFIED  BY  THE  USER.  THE 
C  VOLUME  IS  APPROXIMATED  BY  SUMMING  THE  PRODUCTS  OF  INCREMENTAL  AREAS  OF 
C  THE  DOPPLER  AND  DELAY  BOUNDS  WITH  THE  CLUTTER  DENSITY  FOR  THOSE  AREAS. 

C  THE  ACTUAL  ALGORITHM  IS  TO  TAKE  TWO  ADJACENT  SAMPLES  IN  TIME  OF 
C  THE  MAGNITUDES  OF  TWO  ADJACENT  AMBIGUITY  FUNCTION  CUTS.  THE  VALUE 
C  OF  THE  PROBABILITY  DENSITY  FUNCTION  AT  EACH  OF  THESE  FOUR  POINTS  IS 
C  DETERMINED  AND  MULTIPLIED  WITH  THE  CORRESPONDING  VALUES  OF  THE 
C  AMBIGUITY  FUNCTION  MAGNITUDES  THE  MAXIMUM  VALUE,  MINIMUM  VALUE, 

C  AND  AVERAGE  UALUE  OF  THESE  FOUR  PRODUCTS  IS  DETERMINED  AND  MULTIPLIED 
C  BY  THE  INCREMENTAL  AREA  THIS  RESULTS  IN  AN  INCREMENTAL  VALUE  OF 
C  MAXIMUM  VOLUME.  MINIMUM  VOLUME.  AND  AVERAGE  VOLUME  THE  SAMPLES  OF 
C  THE  MAGNITUDE  OF  THE  AMBIGUITY  FUNCTION  CUTS  ARE  SHIFTED  BY  ONE  AND 
C  THE  PROCESS  IS  REPEATED  UNTIL  ALL  SAMPLES  HAVE  BEEN  CONSIDERED  THEN 
C  THE  AMBIGUTIY  FUNCTION  CUTS  ARE  SHIFTED  BY  ONE  AND  THE  PROCESS  OF 
C  SAMPLING  THE  MAGNITUDES  IS  REPEATED.  AFTER  ALL  DOPPLER  CUTS  HAVE 
C  BEEN  CONSIDERED,  THE  SUM  OF  THE  INCREMENTAL  VOLUMES  GIVES  THE  OVERALL 
C  VOLUME  OF  THE  PRODUCT 
C 

C  '  INTEG  '  MAY  MAKE  CALLS  TO  THE  FOLLOWING  SUBROUTINES 
C  1DC, CPDID, SETUP, PROB.  S,  THE  FFT  LIBRARY  SUBROUTINE  EDFT  LB<DFT5> 

C  'INTEG'  RECEIVES  THE  FOLLOWING  VARIABLES  FROM  THE  PROGRAM  'IMAIN' 

C  Z.  FREQ,  TOTAL.  SAMPLE.  M,  SCALE. 

C  'INTEG'  RETURNS  THE  VARIABLE  'AVGVOL '  TO  'IMAIN'. 

C 

C  'INTEG'  IS  COMPILED  ON  THE  ECLIPSE  USING  THE  COMMAND 
C  FORTRAN  INTEG 

C 

C  THIS  SUBROUTINE  IS  STORED  IN  FILE  'INTEG  FR  '  IN  DIRECTORY  DP A  REED 
C 

C  THE  FOLLOWING  VARIABLES  ARE  USED  IN  'INTEG' 

C  Z  -  ARRAY  CONTAINING  THE  AMPLITUDE  OF  EACH  SAMPLE  OF  THE  CUT 

C  FREQ  -  ARRAY  CONTAINING  THE  FM  OF  EACH  SAMPLE  OF  A  CUT 

C  TOTAL  -  TOTAL  TIME  DURATION  OF  THE  WAVEFORM 

C  SAMPLE-  #  OF  SAMPLES  MADE  OF  WAVEFORM  PLUS  THE  ZERO-FILL 

C  SAMPLES  NEEDED  TO  PERFORM  THE  LINEAR  CONVOLUTION 

C  M  -  »  OF  SAMPLES  OF  WAVEFORM  AFTER  WEIGHTING 

C  SCALE  -  SCALING  FACTOR  OF  MAGNITUDE 

C  WS  -  STARTING  VALUE  OF  DOPPLER  FREQUENCIES  FOR  CLUTTER  t.  AMB  FCN 

C  WF  -  FINAL  VALUE  OF  DOPPLER  FREQS  FOR  CLUTTER  S>  AMB  FCN 

C  WIN  -  INCREMENT  OF  DOPPLER  FREQ  BETWEEN  WS  *  WF 

C  NUMW  -  *  OF  DOPPLER  FREQS  TO  BE  CONSIDERED 

C  APEA  -  INCREMENTAL  AREA  OF  DOPPLER  AND  DELAY ;  <  T I ME  (TAU)  BETWEEN 

C  SAMPLES  OF  AMB  FCN  MAG >*< INTERVAL  BETWEEN  DOPPLER  FREQS) 

C  MAXVOL-  MAXIMUM  VALUE  OF  THE  VOLUME  OF  THE  PRODUCT  OF  THE 
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C  AMG  FCN  MAC  6  CLUTTER  DENSITY 

C  MINVOL  -MINIMUM  VALUE  OF  THE  VOLUME 

C  AVGVOL-  AVERAGE  VALUE  OF  THE  VOLUME 

C  IP  -  USED  TO  I  DENT  I TFY  NON-STANDARD  CLUTTER  PPOB  DENSITY  FCN 

C  CTS  -  LOWER  BOUND  OF  TAU  FOR  CLUTTER 

C  CTF  -  UPPER  BOUND  OF  TAU  FOR  CLUTTER 

C  DCWF  -  ARRAY  CONTAINING  FINAL  VALUE  OFDOPPLER  FOR  EACH  CLUTTER  AREA 

C  DCWS  -  ARRAY  CONTAING  STARTING  DOPPLER  VALUE  FOR  EACH  CLUTTER  AREA 

C  DCTS  -  ARRAY  CONTAINING  STARTING  TAU  VALUE  FOR  EACH  CLUTTER  AREA 

C  DCTF  -  ARRAY  CONTAINING  FINAL  TAU  VALUE  FOR  EACH  CLUTTER  AREA 

I  C  WMEAN  -  STATISTIC  USED  TO  DESCRIBE  STANDARD  PROB  DENSITY  FCN 

i  C  WSTDV-  STATISTIC  OF  U  USED  TO  DESCRIBE  A  STANDARD  PROB  DENSITY  FCN 

C  TMEAN  -  STATISTIC  OF  TAU  USED  TO  DESCRIBE  A  STANDARD  PROB  DENSITY  FCN 

C  TSTDV  -  STATISTIC  OF  TAU  USED  TO  DESCRIBE  A  STANDARD  PROB  DENSITY  FCN 

C  W  -  DOPPLER  FREQ  OF  CUT  TO  BE  MADE 

C  VALUE  -  IMAGINARY  VALUE  *  W 

C  XVALS  -  ARRAY  CONTAINING  THE  TIME-REVERSED.  CONJUGATE  SAMPLES  OF  THE 

C  COMPLEX  ENVELOPE  OF  WAVEFORM;  LATER  CONTAINS  SAMPLES 

C  OF  THE  CUT  OF  THE  AMB  FCN 

C  YVALS  -  ARRAY  CONTAINING  THE  DOPPLER  AND  TIME-DELAYED  SAMPLES  OF  THE 

,  C  COMPLEX  ENVELOPE  OF  WAVEFORM 

C  MAG  -  ARRAY  CONTAING  THE  MAGNITUDE  SAMPLES  OF  A  CUT  OF  AMB  FCN 

C  TIME  -  VALUE  OF  T I ME  <  T AU  > 

C  DENSITY-VALUE  OF  THE  CLUTTER  PROB  DENSITY  FCN  FOR  A  SPECIFIED  W  Sr  TAU 

C  CSUM  -  TOTAL  VALUE  OF  THE  PROB  DENSITY  FCN  FOR  NON-STANDARD 

C  CLUTTER  DISTRIBUTION 

<i'  C  PROD  -  ARRAY  CONTAINING  THE  PRODUCT  OF  THE  PRODUCT  OF  THE 

C  AMB  FCN  MAC  V  CLUTTER  DENSITY 

c  MIN  -  USED  TO  STORE  MINIMUM  VALUE  OF  THE  FOUR  SAMPLES 

C  MAX  -  USED  TO  STORE  MAXIMUM  VALUE  OF  THE  FOUR  SAMPLES 

■I  C  I CDY  -  IDENTIFIES  DELAY  PROB  DENSITY  FCN  FOR  CLUTTER 

C  ICDP  -  IDENTIFIES  DOPPLER  PROB  DENSITY  FCN  FOR  CLUTTER 

SUBROUTINE  INTEG<2.  FREQ.  TOTAL.  SAMPLE, M,  AVGVOL, SCALE) 

COMPLEX  XVALS < 1024),  YVALS < 1024),  Z( 1024),  VALUE 
INTEGER  SAMPLE 

REAL  MAX.  MIN.  MAX VOL,  MINVOL,  MAG ( 1024).  PROD < 1024, 2) , FREQ < 1024) 

*  . DCWSUO),  DCWF ( 10).  DCTS ( 10).  DCTF ( 10) 

.  ACCEPT  "ENTER  STARTING  VALUE  OF  W  FOR  CLUTTER  !<  AMBIGUITY  FCN  ",WS 

ACCEPT  'ENTER  FINAL  VALUE  OF  W  FOR  CLUTTER  Si  AMBIGUITY  FCN  ',WF 
ACCEPT  "ENTER  INTERVAL  BETWEEN  W'S  ".WIN 
NUMW=ANINT( (WF-WS)/WIN) 

AREA»WIN*<  TOTAL/M) 

WRITE*  10.  300)  WS.  WF 
WRITEUO.  305)  WIN,  NUMW 
WRITEUO,  310)  AREA 
MAXV0L»0  0 
MINVOL-O  0 
AVGVOL =0  0 

TYPE  'IF  YOUR  CLUTTER  CAN  NOT  BE  MODELED  BY  UNIFORM  OR  GAUSSIAN' 

TYPE  DISTRIBUTIONS  IN  DELAY  S<  DOPPLER.  THEN  YOU  MUST  ENTER  THE' 

TYPE  DISTRIBUTION  POINT-BY-POINT  ' 


( 


ACCEPT  'ENTER  0  IF  YOU  INFuT  CLUTTER  DIST  POI NT-B Y-POI NT , ELSE  l'.IP 
ACCEPT  ENTER  LOWER  BOUND  ON  TAU  FOR  CLUTTER  '.CTS 
ACCEPT  'ENTER  UPPER  BOUND  ON  TAU  FOR  CLUTTER  ' . CTF 
IF  ( IP  EO  1)  GO  TO  10 

S  CALL  IDC  (  TOTAL,  WS,  WF,  N,  DCWF,  DCWS.  DCTS,  DCTF) 

CSUM=0  0 
GO  TO  15 

10  TYPE  ENTER  0  IF  CLUTTER  UNIFORMLY  DISTRIBUTED  IN  DELAY,  1  IF' 

ACCEPT  'DISTRIBUTED  AS  GAUSSIAN ', ICDY 

TYPE  ENTER  0  IF  CLUTTER  IS  DISTRIBUTED  UNIFORMLY  IN  DOPPLER, ' 

ACCEPT  '  1  IF  DISTRIBUTED  AS  GAUSSIAN ',  ICDP 
CALL  CPDIDI ICDY,  ICDP, WMEAN, WSTDV,  TMEAN, TSTDV) 

15  DO  90  IA-l.NUMW 

W=WS-*-WIN*(  IA-1  ) 

WRITE* 10, 315)  W 
VALUE=CMPLX(0.  0,  W) 

CALL  SETUP (Z, FREQ,  XVALS, YVALS,  VALUE,  TOTAL,  SAMPLE, M) 

CALL  DFT5<  XVALS-  SAMPLE,  0>- 
CALL  DFT5< YVALS,  SAMPLE- 0) 

DO  20  1=1, SAMPLE 

XVALS ( I )=CONJG<  XVALS! I ) ) 

YVALS ( I ) =CON JG ( YVALS  < I ) ) 

XVALS* I)*XVALS< I)*YVALS( I ) /SAMPLE 
20  CONTINUE 

CALL  DFT5< XVALS,  SAMPLE,  1 ) 

DO  25  1=1, SAMPLE 
MAG < I ) =C ABS  < XVALS  < I ) > #SC ALE 
25  CONTINUE 

DO  65  1*1, SAMPLE 
IF  < I  GE.  M+l )  GO  TO  30 
T I  ME=*  TOTAL/M**  1-1  )  ) -TOTAL 
GO  TO  35 

30  TIME=TOTAL/M*< I-M-l > 

35  IF  < IP  NE.  0)  GO  TO  55 

DO  50  K* 1 , N 

IF  <(W  GE  DCWS(K)>  AND  <W  LE  DCWF(K)))  GO  TO  40 
DENSI TY=0  0 
GO  TO  50 

40  IF  ((TIME  GE  DCTS  <  K  >  > .  AND  (TIME  LE  DCTF(K)))  GO  TO  45 

DENSITY=0  0 
GO  TO  50 

45  WRITE* 10, 340)  W, TIME 

ACCEPT  'ENTER  VALUE  OF  CLUTTER  DENSITY  AT  THIS  POINT  '.DENSITY 
C  SUM=CSUM+DENS I TY 
IF  (CSUM.  LE.  1.  0)  CO  TO  50 

TYPE  'EXCEEDED  LIMIT  FOR  CLUTTER  DISTRIBUTION.  .  ABORT  AND  START  AGAIN' 
50  CONTINUE 

GO  TO  60 

55  CALL  PROB  (  DENSITY,  TIME.  W.  WS,  WF,  CTS,  CTF.  WMEAN,  WSTDV,  TMEAN,  TSTDV.  ICDY, 

*  ICDP) 

60  PROD* I. 1 >=MAC( I >*MAC( I )*DENSITY 

65  CONTINUE 

IF  (  I A  EQ.  1  )  GO  TO  83 
DO  SO  1=1, SAMPLE- 1 
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SUM=0  0 
MAX=PROD  < I,  1 ) 

MIN=PROD< I,  1 > 

DO  75  K* 1 1  2 
DO  70  L=I. I  +  l 

IF  <PROD)L,K>  LT  MIN)  M1N=PR0D(L.  K) 

IF  (PROD<L.  K)  GT  MAX)  MAX  =  PROD)L,  K> 

SLM=SUM-t-PROD)L,  K) 

70  CONTINUE 

75  CONTINUE 

AVG—SUM/4.  0 
MAX VOL=MAX VOL*MAX  *AREA 
MINVOL=MINVOL-*-MIN*AREA 
AVGVOL=AVGVOL«-AVG*AREA 
80  CONTINUE 

WRITE) 10,  320)  IA 

WRITE) 10.  400)  I  A, MAX VOL.  MINVOL,  AVGVOL 
83  DO  85  1=1. SAMPLE 

PROD) I.  2)=PR0D( I,  I ) 

85  CONTINUE 

<?0  CONTINUE 

WRITE) 10. 325)  MINVOL 
WRITE) 10.  330)  MAX VOL 
WRITE) 10, 335)  AVGVOL 
RETURN 

300  FORMAT) -STARTING  VALUE  OF  W*',F12  2, 3X.  'FINAL  VALUE  OF  W--.F12  2) 

305  format) -interval  between  w  is-'.fiz  2,  3x,  -t»  of  w  to  be  used-  -,i5> 

310  FORMAT) 'AREA- F14  8) 

315  FORMAT)  'W  NOW  BEING  CONSIDERED- '.  FI 2.  2 ) 

320  FORMAT)  'NOW  COMPLETED  CALCULATIONS  FOR  W  INTERVAL- ',  1 5 ) 

325  FORMAT) 'THE  MINIMUM  VOLUME  OF  THE  INTEGRATED  PRODUCT  IS  '.ElO  4) 

330  FORMAT) 'THE  MAXIMUM  VOLUME  OF  THE  INTEGRATED  PRODUCT  IS  ',E10.  4) 

335  FORMAT)  'THE  AVERAGE  VOLUME  OF  THE  INTEGRATED  PRODUCT  IS  ',E10.  4) 

340  FORMAT)  'W»',F12  2,  5X,  'TAU= ' , FI  2.  2 ) 

400  FORMAT)  '  IA=  '.  14,  ZX,  MAX=  '  ,  ElO  4.  2X.  'MIN-  ',  ElO  4.  2X.  'AVG»  ',  ElO  4) 

END 
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*  SUBROUTINE  CPDID  * 

*  WRITTEN  BY  CAPT  JOHN  REED  * 

*  15  APRIL  1 982  * 


C  'CPDID'  IS  A  FORTRAN  V  SUBROUTINE  WHICH  DETERMINES  THE  STATISTICS 
C  NEEDED  TO  DEFINE  A  STANDARD  PROBABILITY  DENSITY  FUNCTION  OF  CLUTTER 

c 

C  'CPDID'  MAY  BE  CALLED  BY  SUBROUTINE  'INTEG' 

C  CPDID'  RECEIVES  THE  FOLLOWING  VARIABLES  FROM  'INTEG'  IP. ICDY. 1CDP 
C  CPDID'  RETUP "J  THE  FOLLOWING  VARIABLES  TO  'INTEG' 

C  WMEAN,  WS  DV,  TMEAN,  TSTDV 

C 

C  CPDID'  IS  STORED  IN  FILE  'CPDID  FR '  IN  DIRECTORY  DP4  REED 

C 

C  'CPDID'  IS  COMPILED  ON  THE  ECLIPSE  USING  THE  COMMAND: 

C  FORTRAN  CPDID 

C 

C  THE  FOLLOWING  VARIABLES  ARE  USED  IN  'CPDID': 

C  IP  -  IDENTIFIES  NON-STANDARD  PROB  DENSITY  FCN 

C  WMEAN  -  STATISTIC  OF  DOPPLER  USED  TO  DEFINE  DENSITY  FCN 

C  WSTDV  -  STATISTIC  OF  DOPPLER  USED  TO  DEFINE  DENSITY  FCN 

C  TMEAN  -  STATISTIC  OF  TAU  USED  TO  DEFINE  DENSITY  FCN 

C  TSTDV  -  STATISTIC  OF  TAU  USED  TO  DEFINE  DENSITY  FCN 

C  ICDY  -  IDENTIFIES  DELAY  PROB  DENSITY  FCN  FOR  CLUTTER 

C  I CDP  -  IDENTIFIES  DOPPLER  PROB  DENSITY  FCN  FOR  CLUTTER 

SUBROUTINE  CPDID< ICDY,  ICDP,  WMEAN,  WSTDV, TMEAN,  TSTDV) 

IF ( ICDY  EG  1 )  GO  TO  5 

TYPE  CLUTTER  IS  DISTRIBUTED  UNIFORMLY  IN  DELAY' 

GO  TO  10 

5  TYPE  CLUTTER  IS  DISTRIBUTED  GAUSSIANLY  IN  DELAY' 

ACCEPT  'ENTER  THE  MEAN  OF  DELAY ', TMEAN 
ACCEPT  'ENTER  THE  STANDARD  DEVIATION  OF  DELAY  '.TSTDV 
10  IF  ( ICDP  EG.  1 )  CO  TO  15 

TYPE  'CLUTTER  IS  DISTRIBUTED  UNIFORMLY  IN  DOPPLER' 

CO  TO  20 

15  TYPE  'CLUTTER  IS  DISTRIBUTED  CAUSSAINLY  IN  DOPPLER' 

ACCEPT  'ENTER  THE  MEAN  OF  DOPPLER  '.WMEAN 
ACCEPT  'ENTER  THE  STANDARD  DEVIATION  OF  DOPPLER  '.WSTDV 
20  RETURN 

END 


C»»««»*«****»«««*t*«*4ft««**4«**««*«**ft«4* 

C*  SUBROUTINE  IDC  * 

C*  WRITTEN  BY  CART  JOHN  REED  * 

C*  15  APRIL  1982  * 

£*##****##*#****#**■********■*#*#*■*#■****#** 

C  'IDC'  IS  A  FORTRAN  V  SUBROUTINE  WHICH  DEFINES  THE  CLUTTER  AREAS 
C  (IN  TERMS  OF  DOPPLER ( W)  S<  DELAY (TAU)  FOR  NON-STANDARD  CLUTTER 
C  DISTRIBUTIONS. 

C 

C  IDC'  MAY  BE  CALLED  BY  SUBROUTINE  '  I NTEG ' 

C  'IDC'  RECEIVES  THE  FOLLOWING  VARIABLES  FROM  ' 1 NTEG ' 

C  TOTAL,  WS.  WF. 

C  'IDC'  RETURNS  THE  FOLLOWING  VARIABLES  TO  'INTEG'. 

C  N,  DCWF,  DC  WS,  DCTS,  DCTF. 

C 

C  'IDC'  IS  COMPILED  ON  THE  ECLIPSE  USING  THE  COMMAND 
C  FORTRAN  IDC 

C 

C  IDC'  IS  STORED  IN  FILE  'IDC  FR '  IN  DIRECTORY  DP4REED 

C 

C  THE  FOLLOWING  VARIABLES  ARE  USED  IN  'IDC': 

C  TOTAL  -  TOTAL  TIME  DURATION  OF  WAVEFORM 

C  WS  -  STARTING  VALUE  OF  DOPPLER  FREQ  FOR  CLUTTER  S.  AMB  FCN 

C  WF  -  FINAL  VALUE  OF  DOPPLER  FREQ  FOR  CLUTTER  S.  AMB  FCN 

C  DCWF  -  ARRAY  CONTAINING  FINAL  DOPPLER  FREQ  FOR  EACH  CLUTTER  AREA 

C  DCWS  -  ARRAY  CONTAINING  STARTING  DOPPLER  FREQ  FOR  EACH  CLUTTER  AREA 

C  DCTS  -  ARRAY  CONTAINING  STARTING  VALUE  OF  TAU  FOR  EACH  CLUTTER  AREA 

C  DCTF  -  ARRAY  CONTAINING  PINAL  TAU  FOR  EACH  CLUTTER  AREA 

C  N  -  *  OF  DIFFERENT  CLUTTER  AREAS 

SUBROUTINE  IDC  (TOTAL-  WS.  WF,  N,  DCWF,  DCWS,  DCTS,  DCTF) 

REAL  DCWS( 10), DCWF( 10), DCTS( 10), DCTF (10) 

ACCEPT  ENTER  #  OF  AREAS  OF  CLUTTER  ',N 
DO  05  1=1, N 

ACCEPT  'ENTER  STARTING  VALUE  OF  W  FOR  NEXT  CLUTTER  AREA  ',DCWS(D 
ACCEPT  -ENTER  FINAL  VALUE  OF  W  FOR  THIS  CLUTTER  AREA  ',DCWF(I> 

ACCEPT  'ENTER  STARTING  VALUE  OF  TAU  FOR  THIS  CLUTTER  AREA  ',DCTS(I> 
ACCEPT  ENTER  FINAL  VALUE  OF  TAU  FDR  THIS  CLUTTER  AREA  ',DCTF(I> 

IF  (DCWSU)  LT  WS)  CO  TO  5 
GO  TO  10 

5  WRITE! 10. 300)  DCWS(I).WS 

DCWS( 1 >=WS 

IF  <DCUS(I>  GT  DCWF ( I ) )  DCWS U > »DCWF U > 

10  IF  ( DCWF  U  >  GT  WF)  CO  TO  15 

GO  TO  20 

15  WRITE* 10. 305)  DCWF(I),WF 

DCWF ( I )=WF 

IF  (DCWSU)  GT  DC  WF  U  >  )  DCWS  (  I  )  =99999 
20  IF(DCTSU)  LT  -TOTAL)  CO  TO  25 

GO  TO  30 

25  WRITE! 10. 310)  DC TS (I ) . -TOTAL 
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DOTS  < I ) “-TOTAL 

IF  ! DCTS! I )  GT.  DCTF < I ) )  DCTS ! I ) =DCTF ! I) 
IF  1DCTF!I)LT  TOTAL)  GO  TO  35 
WR1TE< 10. 315)  DCTF < I ). TOTAL 
DCTF( I ) -TOTAL 

IF  (DCTS(I)  GT  DCTF(I) )  DCTS( 1 1 =99999. 

CONTINUE 

RETURN 


FORMAT!  'DCWS!  I )»'.  F12.  2.  5X.  'LESS  THAN  US.  THUS  MAKING  DCWS 1 1 ) ■ ' , Fl 2  2 ) 
FORMAT!  'DCWF!  I  >»',  F12.  2,  5X,  'EXCEEDS  WF,  THUS  MAKING  DCWF! I >- F12  2) 
FORMAT!  'DCTS! I >»'. F12  10,  5X,  'LESS  THAN  -TOTAL.  MAKING  DCTS! I )“',  Fl2  10) 
FORMAT!  'DCTF!I)»',F12  10,  5X.  'EXCEEDS  TOTAL, MAKE  DCTF 1 1 ) » ' , F 12.  10) 

END 
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C ***************************************** 

C*  SUBROUT INE  PRGB  * 
C*  WRITTEN  BY  CAPT  JOHN  REED  * 
C*  15  APRIL  1983  * 
. . . . . . 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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■PROB'  IS  A  FORTRAN  V  S13ROUTINE  WHICH  DETERMINES  THE  VALUE  OF 
THE  PROBABILITY  DENSITY  FUNCTION  OF  CLUTTER  FOR  A  GIVEN  DELAY ( TAU > 
AND  DOPPLER  FREQUENCY  <W>.  THIS  SUBROUTINE  APPLIES  ONLY  FOR  A 
STANDARD  PROBABILITY  DENSITY  FUNCTION. 

PROB'  IS  CALLED  BY  THE  SUBROUTINE  INTEG' 

PROB'  RECEIVES  THE  FOLLOWING  VARIABLES  FROM  'INTEG' 

TIME,  W,  ICDY,  1C  DP.  WS,  WF.  CTS,  CTF,  WMEAN.  WSTDV.  TMEAN,  TSTDV. 

PROB'  RETURNS  THE  VARIABLE  DENSITY'  TO  'INTEG' 

PROB'  IS  COMPILED  ON  THE  ECLIPSE  USING  THE  COMMAND: 

FORTRAN  PROB 

PROB'  IS  STORED  IN  FILE  PROB  FR '  IN  DIRECTORY  DP4  REED. 

The  FOLLOWING  VARIABLES  ARE  USED  IN  'PROB' 

TIME  -  VALUE  OF  DELAY(TAU)  FOR  WHICH  DENSITY  IS  CALCULATED 
W  -  DOPPLER  FREQ  FOR  WHICH  DENSITY  IS  CALCULATED 
CTS  -  STARTING  VALUE  OF  TAU  FOR  CLUTTER 
CTF  -  FINAL  VALUE  OF  TAU  FOR  CLUTTER 
IP  -  IDENTIFIES  NON-STANDARD  PROB  DENSITY  FCN 
WF  -  FINAL  VALUE  OF  DOPPLER  FREQ  FOR  CLUTTER 
WS  -  STARTING  VALUE  OF  DOPPLER  FREQ  FOR  CLUTTER 
DENSITY-VALUE  OF  THE  PROB  DENSITY  FCN  FOR  GIVEN  W  &  TAU 
TMEAN  -  STATISTIC  OF  TAU  USED  TO  DESCRIBE  PROB  DENSITY  FCN 

TSTDV  -  STATISTIC  OF  TAU  USED  TO  DESCRIBE  PROB  DENSITY  FCN 

WMEAN  -  STATISTIC  OF  W  USED  TO  DESCRIBE  PROB  DENSITY  FCN 

WSTDV  -  STATISTIC  OF  W  USED  TO  DESCRIBE  PROB  DENSITY  FCN 

SCALE  -  SCALING  FACTOR  OF  DENSITY 

ICDY  -  IDENTIFIES  DELAY  PROB  DENSITY  FCN  FOR  CLUTTER 
I CDP  -  IDENTIFIES  DOPPLER  PROB  DENSITY  FCN  FOR  CLUTTER 
DELAYDEN-  VALUE  OF  DELAY  PROB  DENSITY  FCN  FOR  GIVEN  TAU 
DOPDEN-  VALUE  OF  DOPPLER  PROB  DENSITY  FCN  FOR  GIVEN  W 


SUBROUTINE  PROB ( DENSI TY.  T I  ME, W, WS, WF, CTS, CTF,  WMEAN, WSTDV. 
#  TMEAN, TSTDV, ICDY, ICDP) 


IF  ((TIME  LT  CTS)  OR  (TIME  GT  CTF) >  GO  TO  35 
IF  ( ICDY  EQ  1)  GO  TO  5 
DELAYDEN-1  OZ(CTF-CTS) 

GO  TO  10 

5  SCAi.E-1  O/SQRT  <  3  0*3  1 4 1  6*TSTDV*TSTDV ) 

AMT  1  *  < T I ME- TMEAN J / TSTDV 
DELAYDEN»SCALE*EXP ( -0  5*AMT 1 *AMT 1 > 

10  IF  (ICDP  EQ  1)  GO  TO  15 

DOPDEN» 1  OJ(WF-WS) 

CO  TO  30 
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15  SCALE* 1  0/ SORT ( 2  0*3.  M16*WSTDV*WSTDV) 

AMTS* ( W-WMEAN ) / WSTDV 
DOPDEN*SCALE*EXP  < -0  5*AMT2*AMT2 ) 

20  DENS I T Y*DELA YDEN*DQP  DEN 

00  TO  30 

25  DENSI TY=0  0 

30  RETURN 

END 
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£  ***#♦*####*##******-►*#**♦*#**********'•*# 

C*  PROGRAM  CLUTR  * 

C*  PRODUCED  DY  CAPT  JOHN  REED  * 

C*  21  MAY  1982  • 

C **************************************** 

C  'CLUTR'  IS  A  FORTRAN  V  PROGRAM  WHICH  DETERMINES  THE  VALUES  OF 
C  the  PROBABILITY  DENSITY  FUNCTION  OF  CLUTTER  AS  A  FUNCTION  OF  DELAY 
C  AND  DOPPLER  THE  USER  SPECIFIES  THE  DENSITY  FUNCTION  (CAN  BE 
C  UNIFORM  OR  GAUSSIAN)  FOR  BOTH  DELAY  AND  DOPPLER.  SINCE  THE  VARIABLES 
C  ARE  CONSIDERED  INDEPENDENT,  THE  CLUTTER  PROB  DENSITY  FCN  WILL  BE 
C  THE  PRODUCT  OF  THE  DELAY  AND  DOPPLER  DENSITY  FUNCTIONS  IT  STORES 
C  QUANT I  ZED  SAMPLES  OF  THE  PROB  DENSITY  FCN  IN  FILE  'CLUTTER'  WHICH 
C  CAN  BE  GRAPHED  USING  PROGRAM  'PLTTRNS 
C 

C  CLUTR'  MAY  CALL  THE  SUBROUTINE  'COMP'.  IF  SO,  IT  PASSES  THE 
C  VARIABLES  TAU, TMEAN.  TSTDV  (OR  W, WMEAN, WSTDV)  TO  'COMP'  AND 
C  RECEIVES  THE  VARIABLE  DEL < OR  DOP )  FROM  'COMP'. 

C 

C  CLUTR'  IS  COMPILED  ON  THE  ECLIPSE  USING  THE  COMMAND: 

C  FORTRAN  CLUTR 

C 

C  'CLUTR'  IS  STORED  IN  FILE  'CLUTR  FR  '  IN  DIRECTORY  DP4  REED 
C 

C  THE  FOLLOWING  VARIABLES  ARE  USED  IN  'CLUTR': 

C  TAU  -  DELAY  VALUE 

C  IDT  -  IDENTIFIES  TYPE  OF  DELAY  DISTRIBUTION 

C  IDD  -  IDENTIFIES  TYPE  OF  DOPPLER  DISTRIBUTION 

c  TMEAN  -  MEAN  VALUE  OF  DELAY  IF  GAUSSIAN  DISTRIBUTED 

C  TSTDV  -  STANDARD  DEVIATION  OF  DELAY  IF  GAUSSIAN  DISTRIBUTED 

C  WMEAN  -  MEAN  VALUE  OF  DOPPLER  IF  GAUSSIAN  DISTRIBUTED 

C  WSTDV  -  STANDARD  DEVIATION  OF  DOPPLER  IF  GAUSSIAN  DISTRIBUTED 

C  W  -  DOPPLER  VALUE 

C  DELDEN-  VALUE  OF  DELAY  DENSITY  FCN  FOR  A  GIVEN  DELAY 

C  DOPDEN-  VALUE  OF  DOPPLER  DENSITY  FCN  FOR  A  GIVEN  DOPPLER 

C  DENSITY-VALUE  OF  CLUTTER  PROB  DENSITY  FCN  FOR  A  GIVEN  DELAY  &  DOPPLER 

C  GRAPH  -  ARRAY  OF  QUANTIZED  SAMPLES  OF  THE  CLUTTER  PROB  DENSITY  FCN 

C  STORE  -  COMPLEX  ARRAY  OF  GRAPH  FOR  A  GIVEN  DOPPLER 


REAL  DENSITY<64,  64), MAX 
COMPLEX  STORE (64) 

INTEGER  GRAPH < 64, 64 ) 

CALL  DFILW( "CLUTTER", IER) 

TYPE  FOR  DFILW,  IER“'.IER 
CALL  CFILW< "CLUTTER”, 2, IER) 

TYPE  FOR  CFILW,  IER=', IER 
CALL  OPEN! 6,  "CLUTTER", 2,  IER) 

TYPE  'POR  OPEN.  IER='.IER 

ACCEPT  'ENTER  STARTING  VALUE  OF  DELAY  ON  GRAPH  ' , GST  ART 
ACCEPT  ENTER  FINAL  VALUE  OF  DELAY  ON  GRAPH  '.GSTOP 
ACCEPT  ENTER  STARTING  VALUE  OF  DELAY  ' , T ST  ART 
ACCEPT  'ENTER  FINAL  VALUE  OF  DELAY  '.TSTOP 
ACCEPT  'ENTER  STARTING  VALUE  OF  DOPPLER  ' , WSTART 
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ACCEPT  ENTER  FINAL  VALUE  OF  DCPPLER  ' , WSTQP 
T I NTERVAL3 < GSTQP-G START )  /64  0 
W I NTERVAL3 ' WS  TOP-WSTART )  /64  0 

ACCEPT  'ENTER  0  IF  DELAY  DISTRIBUTED  UNIFORMLY.  1  IF  GAUSSIAN  IDT 

IF  ( IDT  EQ  1)  GO  TO  5 

TYPE  DELAY  DISTRIBUTED  UNIFORMLY' 

CO  TO  10 

5  TYPE  DELAY  DISTRIBUTED  AS  GAUSSIAN' 

ACCEPT  'ENTER  MEAN  VALUE  OF  DELAY  VTMEAN 
ACCEPT  'ENTER  STANDARD  DEVIATION  OF  DELAY  ' .  TSTDV 
10  ACCEPT  ENTER  0  IF  DOPPLER  DISTRIBUTED  UNIFORMLY.  1  IF  GAUSSIAN  ',  IDD 

IF  < IDD  EQ  1 )  GO  TO  15 
TYPE  DOPPLER  DISTRIBUTED  UNIFORMLY' 

CO  TO  20 

15  TYPE  DOPPLER  DISTRIBUTED  AS  GAUSSIAN' 

ACCEPT  'ENTER  MEAN  VALUE  OF  DOPPLER  ' ,  WMEAN 
ACCEPT  'ENTER  STANDARD  DEVIATION  OF  DOPPLER  '.WSTDV 

20  MAX=0  0 

DO  50  I31,64 

W=WSTART ♦ < I -1 ) *W I NTERVAL 
DO  45  J=l,  64 

TAU*GSTAHT+< J-l ) *T I NTERVAL 

IF  (TAU  LT  TSTART  OR.  TAU  GT  TSTOP )  CO  TO  29 
IF  < IDT  EQ  1 )  GO  TO  25 
DELDEN=*1  0/ <  TSTOP -TST ART ) 

CO  TO  30 

25  CALL  COMP (TAU.  TMEAN,  TSTDV.  DEL) 

DELDEN=DEL 
CO  TO  30 

29  DELDEN-O.  0 

30  IF  < IDD  EQ  1 )  GO  TO  35 

D0PDEN=1.  0/ (WSTOP-WSTART) 

GO  TO  40 

35  CALL  COMP (W.  WMEAN.  WSTDV.  DOP ) 

DOPDEN=DOP 

40  DENSITV< j,  I )=DELDEN*DOPDEN 

IF  ( DENSITY ( U.  I )  CT  MAX)  MAX=DENS1 TY < J,  I ) 

45  CONTINUE 

50  CONTINUE 

DO  65  1-1, 64 
DO  60  J» 1 , 64 
K=*0 

VALUE  =  DENSITY(  J.  D/MAX 
DO  55  L»l. 200 
R-0. 005*L 

IF  (VALUE  GE  R)  K«L 
55  CONTINUE 

GRAPH!  J.  I  >»K 
60  CONTINUE 

65  CONTINUE 

DO  75  1-1, 64 
DO  70  J-l.  64 
STORE ( U ) 3GRAPH ( J. I ) 

70  CONTINUE 


C 
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CALL  WSQL* <6, I ■ 5TCPF- 1 ■ IEP i 
CONTINUE 

CALL  CLOSE (6.  3ER) 

TYPE  FOR  CLOSE,  IER=',IER 
WRITE! 10,  300)  TST  ART,  TSTOP 
WR I TE i 1 0,  305 )  WSTART,  WSTOP 
WRITE! 10,  310)  MAX 

TYPE  "TO  SET  DISPLAY,  RUN  PLTTRNS  CLUTTER/I  64/N'  ON  TEXTRONI X ” 
300  FORMAT!  'TAU  STARTS  AT  '.F12  9,  5X,  'TAU  STOPS  AT  ',F12.  9) 

305  FORMAT  1  'DOPPLER  STARTS  AT  ',F12.  3,  5X,  'DOPPLER  STOPS  AT  '.F12  3) 

310  FORMAT!  MAXIMUM  VALUE  OF  CLUTTER  DI STR I  BUT  I  ON- ' ,  F 1 2  6 ) 

END 
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